This is an open assessment looking at potential health effects of a national fish promotion program in Finland. The details of the assessment are described at Opasnet. This file contains the R code to run the assessment model.

Knit to html for best performance.

Calculation is based on BONUS GOHERR project and its fish health benefit-risk assessment.

What has been done for PFAS assessment?

  1. Amount should be gender and age-specific, because we need to worry about young mothers. Solution: take KKE amounts of fish, and scale those with the Goherr subgroup-specific proportions. DONE. Not age-specific atm.
  2. Amounts should reflect the actual fish consumption in Porvoo. Solution: postpone and use national statistics. DONE
  3. Infant’s dioxin concentration module must be added. Solution: Use Goherr model for expo_indir DONE
  4. The module must be updated to match PFAS as well. Solution: update the module to contain column Exposure_agent. Change body fat parameter to volume of distribution. DONE but not checked.
  5. Find why expo_indir is so much higher than EFSA toxicokinetic assessment. Make an alternative model. Solution: problem was with too high absorption fractions, probably. DONE
  6. Add ERF for PFAS (sum of PFOS, PFHxS, PFOA, PFNA). Solution: Make a new page for ERF of PFAS and add that to adjusted ERFs. First immunology; postpone cholesterol and low birth weight and liver toxicity. DONE
  7. Add case burdens for PFAS outcomes. Solution: give rough estimates for immunology to get started. DONE
  8. Add PFAS concentrations to data. Solution: Look at Porvoo concentrations first and make conc_pfas; combine that with conc_vit. DONE
  9. Add dioxin and MeHg concentrations. DONE
  10. Now that the model runs technically, look through each part to check that it makes sense. DONE
  11. If PAF contains some inputs with Age and some without Age (implying that they apply to any age group), only the explicit age groups remain in the end. Solution: Add RRorig to compensate for this. DONE
  12. TEQ and MeHg are missing from dose. Check why and correct. DONE
  13. The proportions of All-cause mortality and CHD2 mortality are 3:1 in Kotimaisen kalan edistämisohjelma but 3:2 in this analysis. Why? Probably omega3-concentration is higher here. Solution: Consider double counting, should the overlap be handled in some way? Also, CHD2 mortality ERF is non-linear, and here these is no background intake. Solution: add background. DONE
  14. In exposure, ALA + DHA < Omega3. Why? Because the values come directly from Fineli and probably omega3 contains other fatty acids as well. DONE
  15. BoDattr: age group 45 - 50 is way smaller than any other group. Why? Solution: it’s the misplaces 5-9-year age group. DONE
  16. conc_mehg does not contain Iter although it is probabilistic. Solution: Update formula. DONE
  17. TEQ response for cancer morbidity is way too large. Why? Is there a unit problem? Let’s look at age group 45-49. BoD 3800, BoDattr 970, dose (Female) 188, CSF 0.0005, incidence 0.00657, population (Female) 154000. This makes case burden = BoD / I / P = 3800/0.00657/154000 = 3.7 DALY/case which is ten times higher than Goherr case burden 0.39 DALY/case. The reason is that the former is about life-time exposure while the latter is per one-year exposure, which is the correct way. The current PAF implies that 21 and 31 % of cancers in females and males, respectively are caused by TEQ in all age groups. That’s an artifact also from the lifetime vs yearly mismatch. DONE
  18. Related to cancer yearly-lifetime mismatch, also case burden must be corrected. No need, because data comes from InpBoD. DONE
  19. MeHg exposure is ca. 2 ug/d when you would assume 1 mg/kg * 0.03 kg/d = 30 ug/d. Where is the difference? Because average fish only 0.1 mg/kg. In addition, salmon is missing. DONE
  20. Sperm concentration and dental defect show fairly large BoD. This is because the ERFs have exposure units pg/g, while the exposure has pg/d. Solution: implement the expo_indir submodule. DONE
  21. Why is average fish so much higher with PFAS concentration than other fish species and the expectation? Because individual samples are accidentally summed, not bootstrapped. Solution: update formula for conc_pfas. DONE
  22. Incorrect levels(Age). DONE
  23. expo_indir for TEQ is still unrealistically high. Solution: adjust parameters lower. DONE
  24. RR has column Exposure: Direct. This is redundant? Well, it comes from dose (and from exposure). It makes no harm. DONE
  25. PAF of TEQ TWI is -7 and TWI 2018 is 2. How come? ERF, incidence, and case_burden seem to be OK. Solution: There are parentheses missing from formula of PAF. Corrected. DONE
  26. Take EFSA food database and derive Finnish fish consumption distributions. DONE
  27. Age group undefined should go to the respective ages. Solution: Exposure should be age-specific. Then “undefined” could be removed. Also, correct incidence not to contain “undefined”. DONE
  28. PFAS BoD (BoDattr?) TWI is 0. Why? Some people should exceed it. Or is it due to lack of variability? Solution: Not any more. I don’t know what changed. DONE
  29. TEQ sperm concentration effect is surprisingly large. Why? Solution: Probably was due to age group mismatch. Now seens OK. DONE
  30. Exposure graph should be age-specific. DONE
  31. Fish consumption graph should be age-specific. DONE
  32. Many graphs get numbers wrong when openv$N>1. Solution: use trim. DONE
  33. Breast cancer, Vitamin D recommendation, and depression go to positive side although they are benefits. Why? It is a plotly thing, as the bars continue from the previous bar and therefore stack in a strange way. Solution: use ggplot for publishable graphs. DONE

What still needs to be done?

  1. MeHg IQ effect is surprisingly large. Why? The exposure increased as salmon and rainbow trout concentrations became available, but still it is bigger than expected. PAF is 110 % so this cannot be. The reason is that the response is based on mother’s intake, and not it is allocated to the child. Solution: adjust exposure to match mother’s exposures.
  2. MeHg concentration is lognormal and logSD may be > 1.7, i.e. there may easily be >1000-fold concentration differences. Solution: there should be an upper limit for concentrations.
  3. All-cause and CHD mortality are double counting. Solution: subtract CHD from all-cause.
  4. Remove PFAS and dioxin TWI from the final graphs.
  5. Vitamin D goes to positive side although it is a benefit.
  6. Fish-specific concentrations are constants. Why?
objects.latest("Op_en2261", code_name="RR2") # RR on page [[Health impact assessment]]
## Loading objects:
##   RR
RRorig <- RR@formula

RR@formula <- function(...) {
  out <- RRorig()
  out <- out * Ovariable(
    output=data.frame(Age = ages, Result=1),
    marginal=c(TRUE,FALSE)
  )
  return(out)
}

objects.latest("Op_en1836", code_name="EFSA_amount") #  [[Domestic fish consumption of the general population in Finland]]
## Loading objects:
##   amount
amountOrig <- amount@formula
amount@formula <- function(...) {
  out <- amountOrig()
  out <- out * fish_proportion
  return(out)
}
# This code was forked from https://github.com/jtuomist/fishhealth/blob/master/fishhealth.Rmd
# This code was previously forked from code Op_fi5923/model on page [[Kotimaisen kalan edistämisohjelma]]
# The code was even more previously forked from Op_fi5889/model on page [[Ruori]] and Op_en7748/model on page [[Goherr assessment]]

dat <- opbase.data("Op_fi5932", subset="Malliparametrit")[-1] # [[PFAS-yhdeisteiden tautitaakka]]
dec <- opbase.data("Op_fi5932", subset="Decisions")[-1]
DecisionTableParser(dec)

CTable <- opbase.data("Op_fi5932",subset="CollapseMarginals")
#for(i in 1:ncol(CTable)) {CTable[[i]] <- as.character(CTable[[i]])} # The default is currently character, not factor
CollapseTableParser(CTable)

cat("Laskennassa käytetty data.\n")
## Laskennassa käytetty data.
dat
cat("Tarkastellut päätökset.\n")
## Tarkastellut päätökset.
dec
cat("Aggregoidut marginaalit.\n")
## Aggregoidut marginaalit.
CTable
dummy <- Ovariable("dummy",data=data.frame(Age="dummy", Fish="dummy", Compound="dummy", Area="dummy", Result=1)) # Keep these columns marginals

fish_proportion <- Ovariable( # How population subgroups eat fish differently
  "fish_proportion",
  data = data.frame(prepare(dat,"fish proportion",c("Type","Exposure_agent","Response","Unit"))),
  unit="proportion of the mean")

total_amount <- Ovariable(
  "total_amount",
  data=prepare(dat, "amount", c("Type","Response","Exposure_agent","Unit")),
  unit="M kg/a")

amount_statistics <- Ovariable( # This comes from the Finnish fish market statistics
  "amount",
  dependencies = data.frame(Name=c("total_amount", "fish_proportion")),
  formula = function(...) {
    amount <- total_amount
      # Filleted weight, i.e. no loss.
    amount <- amount * 1000 / 5.52 /365.25 
      # M kg/a per 5.52M population --> g/d per average person.
    amount <- amount * fish_proportion
      # fish_proportion tells the relative amount in each subgroup
  
    # Match KKE-classification in amount with Fineli classification
    tmp <- Ovariable(
      output = data.frame(
        Kala = c("Kasvatettu", "Kaupallinen", "Kirjolohi", "Silakka", "Vapaa-ajan", "Muu tuonti", "Tuontikirjolohi", "Tuontilohi"),
        Fish = c("Whitefish", "Average fish","Rainbow trout", "Herring", "Average fish", "Average fish", "Rainbow trout", "Salmon"),
        Result = 1
      ),
      marginal = c(TRUE, TRUE, FALSE)
    )
    
    amount <- amount * tmp
    
    return(amount)
  },
  unit="g/d"
)
# Exposure:To child and To eater not needed, because dioxins are not (yet) included

#amount <- EvalOutput(amount)
#amount_statistics <- EvalOutput(amount_statistics)

#ggplot(amount@output, aes(x=amountResult+0.1, colour=Group))+stat_ecdf()+
#  facet_wrap(~Fish)+coord_cartesian(xlim=c(1,300))+
#  scale_x_log10()

population <- Ovariable(
  "population", 
  data=prepare(dat, "population", c("Type", "Exposure_agent", "Response","Unit")),
  unit = "#"
)
population@data$Age <- factor(population@data$Age, levels=levels(ages))

incidence <- Ovariable(
  "incidence",
  data = prepare(dat,"incidence",c("Type","Exposure_agent","Unit")),
  unit="1/person-year")
incidence@data$Age <- factor(incidence@data$Age, levels=levels(ages))

case_burden <- Ovariable(
  "case_burden",
  data = prepare(dat,"case burden",c("Type", "Exposure_agent","Unit")),
  unit="DALY/case")

ERFchoice <- Ovariable(
  "ERFchoice",
  data = 
    prepare(dat, "ERFchoice", c("Unit", "Type"))
)

#### InpBoD

tmp <- prepare(dat, "BoD", c("Type","Exposure_agent","Unit"))
tmp[tmp=="All causes"] <- "All-cause mortality"
tmp[tmp=="Depressive disorders"] <- "Depression"
tmp[tmp=="Neoplasms"] <- "Cancer morbidity yearly"
tmp[tmp=="Respiratory infections and tuberculosis"] <- "Immunosuppression" # Infections of 0-9-year-olds are assumed to represent the background BoD of immunosuppressive effect of PFAS
tmp[tmp=="Cardiovascular diseases"] <- "CHD2 mortality"
tmp$Age <- factor(tmp$Age, levels=levels(ages))
InpBoD <- EvalOutput(Ovariable( # Evaluated because is not a dependency but an Input
  "InpBoD",
  data = tmp,
  unit="DALY/a"
))
conc_vit <- Ovariable(
  "conc_vit",
  ddata = "Op_en1838", # [[Concentrations of beneficial nutrients in fish]]
  subset = "Fineli data for common fish species",
  unit="ALA mg/g, DHA mg/g, Fish g/g, Omega3 mg/g, Vitamin D µg/g f.w. after adjustment"
)
  df = conc_vit@data
  df$Nutrient[df$Nutrient=="D-vitamiini (µg)"] <- "Vitamin D"
  df$Nutrient[df$Nutrient=="rasvahapot n-3 moni-tyydyttymättömät (g)"] <- "Omega3"
  df$Nutrient[df$Nutrient=="rasvahappo 18:3 n-3 (alfalinoleenihappo) (mg)"] <- "ALA"
  df$Nutrient[df$Nutrient=="rasvahappo 22:6 n-3 (DHA) (mg)"] <- "DHA"
  df$Nutrient[df$Nutrient=="proteiini (g)"] <- "Fish"
  df$conc_vitResult[df$Nutrient=="Fish"] <- "1"
  df <- dropall(df[df$Nutrient %in% c("Vitamin D", "Omega3", "ALA", "DHA", "Fish") , ])
conc_vit@data <- df

######## Concentration of PFAS

# Data from EU-kalat3 (Finland excl Vanhankaupunginlahti): # pg/g fresh weight
#       POP     mean       sd      min   Q0.025   median   Q0.975      max
# 2.5% PFOS 2055.757 1404.045 305.0399 330.1365 1533.269 5029.697 5814.935

# Data from EU-kalat3 (Vanhankaupunginlahti, Helsinki) # ng/g f.w.
#      POP   mean       sd      min   Q0.025   median   Q0.975      max
#2.5% PFOS 14.428 11.94542 1.499441 1.607789 15.64988 35.32517 38.91994

conc_eukalat <- EvalOutput(Ovariable(
  "conc_eukalat",
  data = data.frame(
    Area = c("Suomi","Helsinki"),
    Compound="PFOS",
    Result=c("2.056 (3.301 - 5.030)", "14.428 (1.499 - 35.325)")),
  unit="ng/g fresh weight"
))

conc_pfas_raw <- Ovariable(
  "conc_pfas_raw",
  ddata="Op_fi5932",subset="PFAS concentrations",
  unit="ng/g f.w.")

conc_pfas <- Ovariable(
  "conc_pfas",
  dependencies = data.frame(Name="conc_pfas_raw"),
  formula = function(...) {
    out <- conc_pfas_raw
    out$Source <- NULL
    out$Area <- "Porvoo"
    out <- expand_index(out, list(Fish=list(Perch=c(
      "Average fish", "Pike","Rainbow trout","Roach", "Salmon", "Vendace", "Whitefish"))))
    out <- out[sample(nrow(out@output)),]
    out <- out[!duplicated(out@output[intersect(colnames(out@output),c("Fish","Compound","Iter"))]),]
    out@marginal <- c(TRUE, TRUE, FALSE, TRUE, TRUE)
    return(out)
  },
  unit="ng/g fresh weight"
)

#sum_pfas <- oapply(conc_pfas, cols=c("Kala","Compound"), FUN=sum)
#tmp <- conc_pfas / sum_pfas
#summary(tmp, marginals="Compound")
#
## This tells that PFOS consists of 71 - 97 % of the four key PFAS, while PFOA, PFNA, and PFHxS consist of 
# 0 - 10 %, 2 - 18 %, and 0 - 9 %, respectively.
# Even if we included the next most abundant congeners, i.e. PFDA and PFUnA, the overall picture would not change.

conc <- Ovariable(
  "conc",
  dependencies = data.frame(
    Name=c("conc_vit", "conc_pfas","conc_pcddf","conc_mehg"),
    Ident = c(NA,NA,"Op_en3104/pcddf_mean","Op_en4004/conc_mehg_allfish")),
  formula = function(...){
    conc_vit <- oapply(conc_vit, cols=c("Kala", "Adjust"),FUN=mean)
    colnames(conc_vit@output)[colnames(conc_vit@output)=="Nutrient"] <- "Compound"

#    conc_pfas <- oapply(conc_pfas, cols=c("Obs","Area"), FUN=mean) # Redundant because handled already
    conc_pfas$Compound[conc_pfas$Compound %in% c("PFOA","PFNA","PFHxS","PFOS")] <- "PFAS"
    conc_pfas <- oapply(conc_pfas, cols="Area", FUN=sum) # Works because only one area
    
    conc_pcddf$Compound <- "TEQ"
    conc_pcddf$Fish[conc_pcddf$Fish=="Baltic herring"] <- "Herring"
    
    conc_mehg <- oapply(conc_mehg[conc_mehg$Area=="Gulf of Finland",], NULL, mean, c("Area","Kala"))
    conc_mehg$Compound <- "MeHg"
    
    out <- OpasnetUtils::combine(conc_vit, conc_pfas, conc_pcddf, conc_mehg)
    return(out)
  }
)

conc_mehg <- Ovariable(
  "conc_mehg",
  dependencies = data.frame(Name="Hg", Ident="Op_en4004/Hg2"),
  formula = function(...) {

    out <- data.frame()
    for(i in 1:nrow(Hg)) {
      if(openv$N>1) {
        tmp <- data.frame(
          Iter=1:openv$N,
          Result=10^rnorm(openv$N, Hg$MeanLog[i], Hg$SDLog[i])
        )
      } else {
        tmp <- data.frame(Result = 10^out$Hg$MeanLog[i])
      }
      out <- rbind(out,
                   data.frame(Area=as.character(Hg$Area[i]),
                              Kala=as.character(Hg$Kala[i]),
                              tmp
                   )
      )
    }
    out <- Ovariable(output=out, marginal=colnames(out) != "Result") * 
      Ovariable(
        output = data.frame(
          Kala = c("ahven","hauki","kampela","kilkki","kirjolohi","lahna","lohi","made","silakka","simput (Cottus sp.)",
                   "särki","turska","vaellussiika","härkäsimppu","kilohaili","kuha","pasuri"),
          Fish = c("Perch","Pike","Flounder","Chaetiliidae","Rainbow trout","Bream","Salmon","Burbot","Herring","Cottidae",
                   "Roach","Average fish","Whitefish","Fourhorn sculpin","Sprat","Pike-perch","White bream"),
          Result = 1
        ),
        marginal=c(TRUE,TRUE,FALSE)
      )
    return(out)
  },
  unit = "ug/g fresh weight"
)

#conc <- EvalOutput(conc, verbose=TRUE)
#View(conc@output)
#View(conc_pcddf@output)
#View(conc_mehg@output)
###################################################################
# Code copied from http://en.opasnet.org/w/Goherr_assessment#

mc2dparam<- list(
  N2 = 10, # Number of iterations in the new Iter
  strength = 50, # Sample size to which the fun is to be applied. Resembles number of observations
  run2d = FALSE, # Should the mc2d function be used or not?
  info = 1, # Ovariable that contains additional indices, e.g. newmarginals.
  newmarginals = c("Exposure"), # Names of columns that are non-marginals but should be sampled enough to become marginals
  method = "bootstrap", # which method to use for 2D Monte Carlo? Currently bootsrap is the only option.
  fun = mean # Function for aggregating the first Iter dimension.
)

expo_background <- Ovariable(
  "expo_background",
  data = prepare(dat, "background exposure", c("Type","Response","Unit")),
  unit = "Vitamin D: µg /day; DHA, EPA, Omega3: mg /day"
)

expo_dir <- Ovariable(
  "expo_dir",
  dependencies=data.frame(Name=c("amount","conc","expo_background")),
  formula = function(...) {
    amount$Fish[amount$Fish %in% c("Imported fish","Fish product","Tuna")] <- "Average fish"
    amount <- amount * Ovariable(output=data.frame(
      Age = ages[-19],
      Group = c("Toddlers","Other children",rep("Adolescents",2), rep("Adults",9), rep("Elderly",5)),
      Result = 1
    ),
    marginal=c(TRUE,TRUE,FALSE))
    out <- amount * conc
    colnames(out@output)[colnames(out@output)=="Compound"] <- "Exposure_agent"
    out <- OpasnetUtils::combine(out, expo_background)
#    out <- out * Ovariable(output = data.frame(Age=ages[-19], Result=1), marginal=c(TRUE,FALSE))
    return(out)
  },
  unit = "PCDDF, PCB, TEQ: pg /day; Vitamin D, MeHg: µg /day; DHA, EPA, Omega3: mg /day"
)
#expo_dir <- EvalOutput(expo_dir, verbose=TRUE)

exposure <- Ovariable(
  "exposure",
  dependencies = data.frame(
    Name = c(
      "expo_dir", # direct exposure, i.e. the person eats or breaths the exposure agent themself
      "expo_indir", # indirect exposure, i.e. the person (typically fetus or infant) is exposed via someone else (mother)
      "mc2d" # 2D Monte Carlo function
    ),
    Ident = c(
      NA,
      "Op_en7797/expo_indir2", # [[Infant's dioxin exposure]] # expo_indir
      "Op_en7805/mc2d") # [[Two-dimensional Monte Carlo]]
  ),
  formula = function(...) {
    expo_indir <- expo_indir[expo_indir$Gender =="Female",] # Age must be selected too, if available, and then collapsed if >1 location.
    expo_indir$Gender <- NULL
    expo_indir <- expo_indir[expo_indir$Age %in% c("20 - 24","25 - 29", "30 - 34", "35 - 39", "40 - 44"),]
    expo_indir <- CollapseMarginal(expo_indir, "Age", "mean") # "sample" is not possible if openv$N==0
    expo_indir$Age <- "0 - 4"
    
    expo_dir$Exposure <- "Direct"
    
    out <- OpasnetUtils::combine(expo_dir, expo_indir)
    out <- mc2d(out)
    return(out)
  },
  unit = "PCDDF, PCB, TEQ: (To eater: pg /day; to child: pg /g fat); Vitamin D, MeHg: µg /day; DHA, EPA, Omega3: mg /day"
)

#oempty(c("expo_background", "expo_dir","expo_indir","exposure"))
#exposure <- EvalOutput(exposure, verbose=TRUE)
#View(exposure@output)
objects.latest("Op_en2261",code_name="BoDattr2") # [[Health impact assessment]]
## Loading objects:
##   BoDattr
tryCatch(BoDattr <- EvalOutput(BoDattr, verbose=TRUE))
##  Evaluating BoDattr ...
## 
##  - BoD fetched successfully!
## 
##  - PAF fetched successfully!
## - Evaluating BoD ...
## - - Evaluating incidence ...
## 
##  done(0.01 secs)!
## - - Checking incidence marginals ... Response, Age, incidenceSource recognized as marginal(s).
## - - Processing incidence decisions ... done!
## - - Evaluating case_burden ...
## 
##  done(0.01 secs)!
## - - Checking case_burden marginals ... Response, Iter, case_burdenSource recognized as marginal(s).
## - - Processing case_burden marginal collapses ... done!
## - - Evaluating population ...
## 
##  done(0 secs)!
## - - Checking population marginals ... Gender, Age, populationSource recognized as marginal(s).
## 
## - done(0.76 secs)!
## - Checking BoD marginals ... Response, Age, incidenceSource, Adjust, Gender, populationSource, Iter, BoDSource recognized as marginal(s).
## - Processing BoD inputs ... done!
## - Processing BoD marginal collapses ...
## Warning in oapply(variable, FUN = fun[[i]], cols = cols[[i]], na.rm = TRUE):
## While oapplying BoD, found NAs in indices: Adjust, InpBoDSource. They were
## automatically filled using fillna, which may result in a multiplied population.
## Please check your ovariable before using oapply.
##  done!
## - Processing BoD decisions ... done!
## - Evaluating PAF ...
## 
##  - - dose fetched successfully!
## 
##  - - ERF fetched successfully!
## 
##  - - frexposed fetched successfully!
## 
##  - - P_illness fetched successfully!
## 
##  - - sumExposcen fetched successfully!
## 
##  - - mc2d fetched successfully!
## - - Evaluating dose ...
## 
##  - - - BW fetched successfully!
## - - - Evaluating exposure ...
## 
##  - - - - expo_indir fetched successfully!
## - - - - Evaluating expo_dir ...
## - - - - - Evaluating amount ...
## - - - - - - Evaluating dummy ...
## 
##  done(0 secs)!
## - - - - - - Checking dummy marginals ... Age, Fish, Compound, Area, dummySource recognized as marginal(s).
## 
## ----- done(1.71 secs)!
## - - - - - Checking amount marginals ... Food, Iter, Group, Fish, Gender, fish_proportionSource, amountSource recognized as marginal(s).
## - - - - - Processing amount marginal collapses ... done!
## - - - - - Evaluating conc ...
## 
##  - - - - - - conc_pcddf fetched successfully!
## - - - - - - Evaluating conc_vit ...
## 
##  done(0.1 secs)!
## - - - - - - Checking conc_vit marginals ... Kala, Fish, Nutrient, Iter, conc_vitSource recognized as marginal(s).
## - - - - - - Processing conc_vit decisions ... done!
## - - - - - - Evaluating conc_pfas ...
## - - - - - - - Evaluating conc_pfas_raw ...
## 
##  done(0.07 secs)!
## - - - - - - - Checking conc_pfas_raw marginals ... Fish, Compound, Iter, conc_pfas_rawSource recognized as marginal(s).
## 
## ------ done(2.22 secs)!
## - - - - - - Checking conc_pfas marginals ... Fish, Compound, Iter, conc_pfas_rawSource, Area, conc_pfasSource recognized as marginal(s).
## - - - - - - Processing conc_pfas marginal collapses ... done!
## - - - - - - Evaluating conc_pcddf ...
## 
##  - - - - - - - TEF fetched successfully!
## - - - - - - - Evaluating TEF ...
## 
##  - - - - - - - - TEFversion fetched successfully!
## - - - - - - - - Evaluating TEFraw ...
## 
##  done(0.01 secs)!
## - - - - - - - - Checking TEFraw marginals ... Group, Compound, TEFversion, TEFrawSource recognized as marginal(s).
## 
## ------- done(23.95 secs)!
## - - - - - - - Checking TEF marginals ... Compound, Group, TEFversion, TEFrawSource, TEFSource recognized as marginal(s).
## 
## ------ done(1.05 mins)!
## - - - - - - Checking conc_pcddf marginals ... Fish, conc_pcddfSource recognized as marginal(s).
## - - - - - - Evaluating conc_mehg ...
## 
##  - - - - - - - Hg fetched successfully!
## Warning in rnorm(openv$N, Hg$MeanLog[i], Hg$SDLog[i]): NAs produced
## Warning in rnorm(openv$N, Hg$MeanLog[i], Hg$SDLog[i]): NAs produced
## ------ done(27.69 secs)!
## - - - - - - Checking conc_mehg marginals ... Area, Kala, Iter, Fish, conc_mehgSource recognized as marginal(s).
## 
## ----- done(2.07 mins)!
## - - - - - Checking conc marginals ... Fish, Compound, Iter, concSource recognized as marginal(s).
## - - - - - Evaluating expo_background ...
## 
##  done(0.01 secs)!
## - - - - - Checking expo_background marginals ... Exposure_agent, Gender, Iter, expo_backgroundSource recognized as marginal(s).
## 
## ---- done(2.39 mins)!
## - - - - Checking expo_dir marginals ... Iter, Group, Fish, Gender, fish_proportionSource, amountSource, Age, Exposure_agent, concSource, expo_dirSource recognized as marginal(s).
## - - - - Processing expo_dir marginal collapses ...
## Warning in oapply(variable, FUN = fun[[i]], cols = cols[[i]], na.rm = TRUE):
## While oapplying expo_dir, found NAs in indices: Age. They were automatically
## filled using fillna, which may result in a multiplied population. Please check
## your ovariable before using oapply.
##  done!
## - - - - Evaluating expo_indir ...
## - - - - - Processing expo_dir marginal collapses ... done!
## - - - - - Evaluating t0.5 ...
## 
##  done(0.01 secs)!
## - - - - - Checking t0.5 marginals ... Exposure_agent, Iter, t0.5Source recognized as marginal(s).
## - - - - - Evaluating f_ing ...
## 
##  done(0 secs)!
## - - - - - Checking f_ing marginals ... Exposure_agent, Iter, f_ingSource recognized as marginal(s).
## - - - - - Evaluating f_mtoc ...
## 
##  done(0 secs)!
## - - - - - Checking f_mtoc marginals ... Exposure_agent, Iter, f_mtocSource recognized as marginal(s).
## - - - - - Evaluating BF ...
## 
##  done(0 secs)!
## - - - - - Checking BF marginals ... Exposure_agent, BFSource recognized as marginal(s).
## 
## ---- done(4.87 secs)!
## - - - - Checking expo_indir marginals ... Iter, Gender, Age, Exposure_agent, expo_dirSource, f_ingSource, t0.5Source, f_mtocSource, BFSource, Exposure, expo_indirSource recognized as marginal(s).
## - - - - Processing expo_indir marginal collapses ... done!
## 
## --- done(3.23 mins)!
## - - - Checking exposure marginals ... Iter, Gender, Age, Exposure_agent, Exposure, expo_dirSource, exposureSource recognized as marginal(s).
## - - - Processing exposure marginal collapses ...
## Warning in oapply(variable, FUN = fun[[i]], cols = cols[[i]], na.rm = TRUE):
## While oapplying exposure, found NAs in indices: Gender. They were automatically
## filled using fillna, which may result in a multiplied population. Please check
## your ovariable before using oapply.
##  done!
## - - - Evaluating BW ...
## 
##  done(0 secs)!
## - - - Checking BW marginals ... BWSource recognized as marginal(s).
## 
## -- done(3.7 mins)!
## - - Checking dose marginals ... Iter, Gender, Age, Exposure_agent, Exposure, Scaling, exposureSource, BWSource, doseSource recognized as marginal(s).
## - - Processing dose marginal collapses ... done!
## - - Evaluating ERF ...
## 
##  - - - ERF_env fetched successfully!
## 
##  - - - ERF_omega3 fetched successfully!
## 
##  - - - ERF_mehg fetched successfully!
## 
##  - - - ERF_diox fetched successfully!
## 
##  - - - ERF_vit fetched successfully!
## 
##  - - - ERF_micr fetched successfully!
## 
##  - - - ERF_pfas fetched successfully!
## - - - Evaluating ERF_env ...
## 
##  done(0.08 secs)!
## - - - Checking ERF_env marginals ... Exposure_agent, Response, Subgroup, Exposure, ER_function, Scaling, Exposure_unit, Observation, Iter, ERF_envSource recognized as marginal(s).
## - - - Evaluating ERF_omega3 ...
## 
##  done(0.03 secs)!
## - - - Checking ERF_omega3 marginals ... Exposure_agent, Response, Exposure, Exposure_unit, ER_function, Scaling, Observation, Iter, ERF_omega3Source recognized as marginal(s).
## - - - Evaluating ERF_mehg ...
## 
##  done(0.01 secs)!
## - - - Checking ERF_mehg marginals ... Exposure_agent, Response, Exposure, Exposure_unit, ER_function, Scaling, Observation, Iter, ERF_mehgSource recognized as marginal(s).
## - - - Evaluating ERF_diox ...
## 
##  done(0.02 secs)!
## - - - Checking ERF_diox marginals ... Exposure_agent, Response, Exposure, Exposure_unit, ER_function, Scaling, Observation, Iter, ERF_dioxSource recognized as marginal(s).
## - - - Evaluating ERF_vit ...
## 
##  done(0 secs)!
## - - - Checking ERF_vit marginals ... Exposure_agent, Response, Exposure, Exposure_unit, ER_function, Scaling, Observation, ERF_vitSource recognized as marginal(s).
## - - - Evaluating ERF_micr ...
## 
##  done(0 secs)!
## - - - Checking ERF_micr marginals ... Exposure_agent, Response, Exposure, Exposure_unit, ER_function, Scaling, Observation, ERF_micrSource recognized as marginal(s).
## - - - Evaluating ERF_pfas ...
## 
##  done(0 secs)!
## - - - Checking ERF_pfas marginals ... Exposure_agent, Response, Exposure, Exposure_unit, ER_function, Scaling, Observation, ERF_pfasSource recognized as marginal(s).
## - - - Evaluating ERFchoice ...
## 
##  done(0 secs)!
## - - - Checking ERFchoice marginals ... Exposure_agent, Response, Scaling, Exposure, ER_function, ERFchoiceSource recognized as marginal(s).
## 
## -- done(3.09 mins)!
## - - Checking ERF marginals ... Exposure_agent, Response, Exposure, ER_function, Scaling, Observation, Iter, ERFSource recognized as marginal(s).
## - - Processing ERF marginal collapses ... done!
## - - Evaluating RR ...
## - - - Processing dose marginal collapses ... done!
## - - - Processing ERF marginal collapses ... done!
## 
## -- done(52.66 secs)!
## - - Checking RR marginals ... Exposure_agent, Response, ER_function, Scaling, Iter, ERFSource, Gender, Age, Exposure, doseSource, RRSource recognized as marginal(s).
## - - Evaluating frexposed ...
## 
##  done(0 secs)!
## - - Checking frexposed marginals ... frexposedSource recognized as marginal(s).
## - - Evaluating P_illness ...
## 
##  done(0 secs)!
## - - Checking P_illness marginals ... Response, Illness, Age, P_illnessSource recognized as marginal(s).
## 
## - done(11.25 mins)!
## - Checking PAF marginals ... Exposure_agent, Response, ER_function, Scaling, Iter, ERFSource, Gender, Age, Exposure, doseSource, frexposedSource, incidenceSource, Adjust, RRSource, PAFSource recognized as marginal(s).
## - Processing PAF marginal collapses ...
## Warning in oapply(variable, FUN = fun[[i]], cols = cols[[i]], na.rm = TRUE):
## While oapplying PAF, found NAs in indices: Adjust. They were automatically
## filled using fillna, which may result in a multiplied population. Please check
## your ovariable before using oapply.
##  done!
## 
##  done(12.62 mins)!
##  Checking BoDattr marginals ... Adjust, Response, Age, Gender, Iter, InpBoDSource, Exposure_agent, PAFSource, BoDattrSource recognized as marginal(s).
BoDattr$Age <- factor(BoDattr$Age, levels=levels(ages))

tst <- BoDattr[BoDattr$Response=="CHD2 mortality",] * -1
tst$Response <- "All-cause mortality"
tst <- OpasnetUtils::combine(BoDattr, tst, name = "BoDattr")
oprint(summary(amount,"mean"))
##             Group          Fish Gender        mean
## 1         Elderly         Bream Female  0.07310000
## 2  Other children         Bream Female  0.02125000
## 3     Adolescents  Fish product Female 10.03239193
## 4          Adults  Fish product Female  9.68891209
## 5         Elderly  Fish product Female 13.54361847
## 6         Infants  Fish product Female  0.32917881
## 7  Other children  Fish product Female 12.01013046
## 8        Toddlers  Fish product Female  2.91654229
## 9     Adolescents      Flounder Female  0.04675000
## 10         Adults      Flounder Female  0.15312521
## 11        Infants      Flounder Female  0.06344377
## 12 Other children      Flounder Female  0.03953492
## 13       Toddlers      Flounder Female  0.03331799
## 14    Adolescents       Herring Female  0.30063704
## 15         Adults       Herring Female  0.82153806
## 16        Elderly       Herring Female  2.53177575
## 17        Infants       Herring Female  0.00425000
## 18 Other children       Herring Female  0.23244724
## 19       Toddlers       Herring Female  0.06290565
## 20    Adolescents Imported fish Female  0.02915649
## 21         Adults Imported fish Female  0.87277057
## 22        Elderly Imported fish Female  0.98544293
## 23 Other children Imported fish Female  0.83097009
## 24       Toddlers Imported fish Female  0.31018711
## 25    Adolescents         Perch Female  0.07735000
## 26         Adults         Perch Female  0.28340779
## 27        Elderly         Perch Female  0.70761299
## 28        Infants         Perch Female  0.00340000
## 29 Other children         Perch Female  0.20077449
## 30       Toddlers         Perch Female  0.03904384
## 31    Adolescents        Salmon Female  6.55818983
## 32         Adults        Salmon Female 10.30522875
## 33        Elderly        Salmon Female 14.76496224
## 34        Infants        Salmon Female  0.24475105
## 35 Other children        Salmon Female  3.73026532
## 36       Toddlers        Salmon Female  2.38253281
## 37    Adolescents          Tuna Female  0.53789148
## 38         Adults          Tuna Female  2.79320966
## 39        Elderly          Tuna Female  1.29878515
## 40 Other children          Tuna Female  0.60739863
## 41       Toddlers          Tuna Female  0.04377803
## 42    Adolescents     Whitefish Female  0.15587904
## 43         Adults     Whitefish Female  3.18852976
## 44        Elderly     Whitefish Female  2.57402181
## 45 Other children     Whitefish Female  0.13742324
## 46       Toddlers     Whitefish Female  0.04623433
## 47        Elderly         Bream   Male  0.09890000
## 48 Other children         Bream   Male  0.02875000
## 49    Adolescents  Fish product   Male 13.57323614
## 50         Adults  Fish product   Male 13.10852812
## 51        Elderly  Fish product   Male 18.32371911
## 52        Infants  Fish product   Male  0.44535957
## 53 Other children  Fish product   Male 16.24900004
## 54       Toddlers  Fish product   Male  3.94591016
## 55    Adolescents      Flounder   Male  0.06325000
## 56         Adults      Flounder   Male  0.20716941
## 57        Infants      Flounder   Male  0.08583568
## 58 Other children      Flounder   Male  0.05348842
## 59       Toddlers      Flounder   Male  0.04507728
## 60    Adolescents       Herring   Male  0.40674423
## 61         Adults       Herring   Male  1.11149267
## 62        Elderly       Herring   Male  3.42534367
## 63        Infants       Herring   Male  0.00575000
## 64 Other children       Herring   Male  0.31448744
## 65       Toddlers       Herring   Male  0.08510765
## 66    Adolescents Imported fish   Male  0.03944701
## 67         Adults Imported fish   Male  1.18080724
## 68        Elderly Imported fish   Male  1.33324631
## 69 Other children Imported fish   Male  1.12425366
## 70       Toddlers Imported fish   Male  0.41966491
## 71    Adolescents         Perch   Male  0.10465000
## 72         Adults         Perch   Male  0.38343407
## 73        Elderly         Perch   Male  0.95735875
## 74        Infants         Perch   Male  0.00460000
## 75 Other children         Perch   Male  0.27163608
## 76       Toddlers         Perch   Male  0.05282402
## 77    Adolescents        Salmon   Male  8.87284506
## 78         Adults        Salmon   Male 13.94236831
## 79        Elderly        Salmon   Male 19.97612539
## 80        Infants        Salmon   Male  0.33113377
## 81 Other children        Salmon   Male  5.04682956
## 82       Toddlers        Salmon   Male  3.22342675
## 83    Adolescents          Tuna   Male  0.72773554
## 84         Adults          Tuna   Male  3.77904837
## 85        Elderly          Tuna   Male  1.75717991
## 86 Other children          Tuna   Male  0.82177461
## 87       Toddlers          Tuna   Male  0.05922909
## 88    Adolescents     Whitefish   Male  0.21089518
## 89         Adults     Whitefish   Male  4.31389320
## 90        Elderly     Whitefish   Male  3.48250009
## 91 Other children     Whitefish   Male  0.18592557
## 92       Toddlers     Whitefish   Male  0.06255233
oprint(summary(BoD,marginals=c("Age","Response"),"mean"))
##           Age                                          Response         mean
## 1       0 - 4                           Cancer morbidity yearly 3.102639e+02
## 2     10 - 14                           Cancer morbidity yearly 3.471589e+02
## 3     15 - 19                           Cancer morbidity yearly 4.155185e+02
## 4     20 - 24                           Cancer morbidity yearly 5.688098e+02
## 5     25 - 29                           Cancer morbidity yearly 7.967250e+02
## 6     30 - 34                           Cancer morbidity yearly 1.146693e+03
## 7     35 - 39                           Cancer morbidity yearly 1.678742e+03
## 8     40 - 44                           Cancer morbidity yearly 2.450202e+03
## 9     45 - 49                           Cancer morbidity yearly 3.846950e+03
## 10      5 - 9                           Cancer morbidity yearly 3.441288e+02
## 11    50 - 54                           Cancer morbidity yearly 6.981211e+03
## 12    55 - 59                           Cancer morbidity yearly 1.124565e+04
## 13    60 - 64                           Cancer morbidity yearly 1.634208e+04
## 14    65 - 69                           Cancer morbidity yearly 2.187666e+04
## 15    70 - 74                           Cancer morbidity yearly 2.555144e+04
## 16    75 - 79                           Cancer morbidity yearly 1.800670e+04
## 17    80 - 84                           Cancer morbidity yearly 1.465322e+04
## 18        85+                           Cancer morbidity yearly 1.196704e+04
## 19      0 - 4      Dioxin recommendation tolerable daily intake 3.088021e+01
## 20    10 - 14      Dioxin recommendation tolerable daily intake 3.726333e+01
## 21    15 - 19      Dioxin recommendation tolerable daily intake 3.579681e+01
## 22    20 - 24      Dioxin recommendation tolerable daily intake 3.788100e+01
## 23    25 - 29      Dioxin recommendation tolerable daily intake 4.291754e+01
## 24    30 - 34      Dioxin recommendation tolerable daily intake 4.232425e+01
## 25    35 - 39      Dioxin recommendation tolerable daily intake 4.353255e+01
## 26    40 - 44      Dioxin recommendation tolerable daily intake 4.180625e+01
## 27    45 - 49      Dioxin recommendation tolerable daily intake 3.793277e+01
## 28      5 - 9      Dioxin recommendation tolerable daily intake 3.695710e+01
## 29    50 - 54      Dioxin recommendation tolerable daily intake 4.293069e+01
## 30    55 - 59      Dioxin recommendation tolerable daily intake 4.450858e+01
## 31    60 - 64      Dioxin recommendation tolerable daily intake 4.339222e+01
## 32    65 - 69      Dioxin recommendation tolerable daily intake 4.307138e+01
## 33    70 - 74      Dioxin recommendation tolerable daily intake 4.219912e+01
## 34    75 - 79      Dioxin recommendation tolerable daily intake 2.560163e+01
## 35    80 - 84      Dioxin recommendation tolerable daily intake 1.947431e+01
## 36        85+      Dioxin recommendation tolerable daily intake 1.822112e+01
## 37  Undefined      Dioxin recommendation tolerable daily intake 6.666909e+02
## 38      0 - 4 Dioxin recommendation tolerable daily intake 2018 8.414927e+01
## 39    10 - 14 Dioxin recommendation tolerable daily intake 2018 1.015434e+02
## 40    15 - 19 Dioxin recommendation tolerable daily intake 2018 9.754712e+01
## 41    20 - 24 Dioxin recommendation tolerable daily intake 2018 1.032266e+02
## 42    25 - 29 Dioxin recommendation tolerable daily intake 2018 1.169513e+02
## 43    30 - 34 Dioxin recommendation tolerable daily intake 2018 1.153345e+02
## 44    35 - 39 Dioxin recommendation tolerable daily intake 2018 1.186272e+02
## 45    40 - 44 Dioxin recommendation tolerable daily intake 2018 1.139230e+02
## 46    45 - 49 Dioxin recommendation tolerable daily intake 2018 1.033676e+02
## 47      5 - 9 Dioxin recommendation tolerable daily intake 2018 1.007089e+02
## 48    50 - 54 Dioxin recommendation tolerable daily intake 2018 1.169871e+02
## 49    55 - 59 Dioxin recommendation tolerable daily intake 2018 1.212869e+02
## 50    60 - 64 Dioxin recommendation tolerable daily intake 2018 1.182448e+02
## 51    65 - 69 Dioxin recommendation tolerable daily intake 2018 1.173705e+02
## 52    70 - 74 Dioxin recommendation tolerable daily intake 2018 1.149936e+02
## 53    75 - 79 Dioxin recommendation tolerable daily intake 2018 6.976501e+01
## 54    80 - 84 Dioxin recommendation tolerable daily intake 2018 5.306793e+01
## 55        85+ Dioxin recommendation tolerable daily intake 2018 4.965295e+01
## 56  Undefined Dioxin recommendation tolerable daily intake 2018 1.816748e+03
## 57      0 - 4                                 Immunosuppression 3.024113e+02
## 58      5 - 9                                 Immunosuppression 2.715484e+02
## 59      0 - 4                         Loss in child's IQ points 1.694833e+04
## 60      0 - 4                                          PFAS TWI 2.827461e+02
## 61    10 - 14                                          PFAS TWI 3.411914e+02
## 62    15 - 19                                          PFAS TWI 3.277636e+02
## 63    20 - 24                                          PFAS TWI 3.468469e+02
## 64    25 - 29                                          PFAS TWI 3.929626e+02
## 65    30 - 34                                          PFAS TWI 3.875302e+02
## 66    35 - 39                                          PFAS TWI 3.985938e+02
## 67    40 - 44                                          PFAS TWI 3.827873e+02
## 68    45 - 49                                          PFAS TWI 3.473209e+02
## 69      5 - 9                                          PFAS TWI 3.383874e+02
## 70    50 - 54                                          PFAS TWI 3.930830e+02
## 71    55 - 59                                          PFAS TWI 4.075305e+02
## 72    60 - 64                                          PFAS TWI 3.973089e+02
## 73    65 - 69                                          PFAS TWI 3.943712e+02
## 74    70 - 74                                          PFAS TWI 3.863846e+02
## 75    75 - 79                                          PFAS TWI 2.344142e+02
## 76    80 - 84                                          PFAS TWI 1.783111e+02
## 77        85+                                          PFAS TWI 1.668366e+02
## 78  Undefined                                          PFAS TWI 6.104370e+03
## 79      0 - 4                               Sperm concentration 4.428657e+03
## 80      0 - 4                          Vitamin D recommendation 6.203907e+01
## 81    10 - 14                          Vitamin D recommendation 7.486292e+01
## 82    15 - 19                          Vitamin D recommendation 7.191664e+01
## 83    20 - 24                          Vitamin D recommendation 7.610382e+01
## 84    25 - 29                          Vitamin D recommendation 8.622235e+01
## 85    30 - 34                          Vitamin D recommendation 8.503041e+01
## 86    35 - 39                          Vitamin D recommendation 8.745792e+01
## 87    40 - 44                          Vitamin D recommendation 8.398973e+01
## 88    45 - 49                          Vitamin D recommendation 7.620782e+01
## 89      5 - 9                          Vitamin D recommendation 7.424767e+01
## 90    50 - 54                          Vitamin D recommendation 8.624877e+01
## 91    55 - 59                          Vitamin D recommendation 8.941879e+01
## 92    60 - 64                          Vitamin D recommendation 8.717600e+01
## 93    65 - 69                          Vitamin D recommendation 8.653142e+01
## 94    70 - 74                          Vitamin D recommendation 8.477903e+01
## 95    75 - 79                          Vitamin D recommendation 5.143427e+01
## 96    80 - 84                          Vitamin D recommendation 3.912434e+01
## 97        85+                          Vitamin D recommendation 3.660665e+01
## 98  Undefined                          Vitamin D recommendation 1.339398e+03
## 99      0 - 4                           Yes or no dental defect 3.445681e+02
## 100     0 - 4                               All-cause mortality 4.882678e+03
## 101   10 - 14                               All-cause mortality 1.086325e+03
## 102   15 - 19                               All-cause mortality 2.966719e+03
## 103   20 - 24                               All-cause mortality 5.275181e+03
## 104   25 - 29                               All-cause mortality 6.264269e+03
## 105   30 - 34                               All-cause mortality 6.710701e+03
## 106   35 - 39                               All-cause mortality 8.275831e+03
## 107   40 - 44                               All-cause mortality 1.022931e+04
## 108   45 - 49                               All-cause mortality 1.369141e+04
## 109     5 - 9                               All-cause mortality 9.310611e+02
## 110   50 - 54                               All-cause mortality 2.132050e+04
## 111   55 - 59                               All-cause mortality 2.983913e+04
## 112   60 - 64                               All-cause mortality 3.921433e+04
## 113   65 - 69                               All-cause mortality 4.929973e+04
## 114   70 - 74                               All-cause mortality 5.979016e+04
## 115   75 - 79                               All-cause mortality 4.944779e+04
## 116   80 - 84                               All-cause mortality 5.195164e+04
## 117       85+                               All-cause mortality 6.952864e+04
## 118     0 - 4                                        Depression 4.476772e-01
## 119   10 - 14                                        Depression 5.846342e+02
## 120   15 - 19                                        Depression 1.307743e+03
## 121   20 - 24                                        Depression 1.825022e+03
## 122   25 - 29                                        Depression 1.812272e+03
## 123   30 - 34                                        Depression 1.603272e+03
## 124   35 - 39                                        Depression 1.717558e+03
## 125   40 - 44                                        Depression 1.701713e+03
## 126   45 - 49                                        Depression 1.536503e+03
## 127     5 - 9                                        Depression 7.944539e+01
## 128   50 - 54                                        Depression 1.722057e+03
## 129   55 - 59                                        Depression 1.762131e+03
## 130   60 - 64                                        Depression 1.694530e+03
## 131   65 - 69                                        Depression 1.551205e+03
## 132   70 - 74                                        Depression 1.424708e+03
## 133   75 - 79                                        Depression 7.936395e+02
## 134   80 - 84                                        Depression 5.813108e+02
## 135       85+                                        Depression 5.511287e+02
## 136     0 - 4                                    CHD2 mortality 5.330556e+01
## 137   10 - 14                                    CHD2 mortality 7.548896e+01
## 138   15 - 19                                    CHD2 mortality 1.461496e+02
## 139   20 - 24                                    CHD2 mortality 2.541203e+02
## 140   25 - 29                                    CHD2 mortality 4.384905e+02
## 141   30 - 34                                    CHD2 mortality 6.382844e+02
## 142   35 - 39                                    CHD2 mortality 1.168969e+03
## 143   40 - 44                                    CHD2 mortality 1.977620e+03
## 144   45 - 49                                    CHD2 mortality 3.114102e+03
## 145     5 - 9                                    CHD2 mortality 4.872613e+01
## 146   50 - 54                                    CHD2 mortality 5.473365e+03
## 147   55 - 59                                    CHD2 mortality 8.741711e+03
## 148   60 - 64                                    CHD2 mortality 1.306399e+04
## 149   65 - 69                                    CHD2 mortality 1.811567e+04
## 150   70 - 74                                    CHD2 mortality 2.453937e+04
## 151   75 - 79                                    CHD2 mortality 2.340159e+04
## 152   80 - 84                                    CHD2 mortality 2.792111e+04
## 153       85+                                    CHD2 mortality 4.182868e+04
## 154   15 - 19                                     Breast cancer 2.061577e+00
## 155   20 - 24                                     Breast cancer 7.838659e+00
## 156   25 - 29                                     Breast cancer 3.660033e+01
## 157   30 - 34                                     Breast cancer 1.455554e+02
## 158   35 - 39                                     Breast cancer 2.649923e+02
## 159   40 - 44                                     Breast cancer 4.313909e+02
## 160   45 - 49                                     Breast cancer 7.064921e+02
## 161   50 - 54                                     Breast cancer 1.043878e+03
## 162   55 - 59                                     Breast cancer 1.342853e+03
## 163   60 - 64                                     Breast cancer 1.523321e+03
## 164   65 - 69                                     Breast cancer 1.653962e+03
## 165   70 - 74                                     Breast cancer 1.673425e+03
## 166   75 - 79                                     Breast cancer 1.235337e+03
## 167   80 - 84                                     Breast cancer 9.479436e+02
## 168       85+                                     Breast cancer 8.755721e+02
oprint(summary(BoDattr,marginals=c("Age","Response"),"mean"))
##         Age                                          Response          mean
## 1     0 - 4                           Cancer morbidity yearly  6.894085e-01
## 2     5 - 9                           Cancer morbidity yearly  8.169153e-01
## 3   10 - 14                           Cancer morbidity yearly  1.445137e+00
## 4   15 - 19                           Cancer morbidity yearly  1.754437e+00
## 5   20 - 24                           Cancer morbidity yearly  4.663839e+00
## 6   25 - 29                           Cancer morbidity yearly  6.506471e+00
## 7   30 - 34                           Cancer morbidity yearly  9.186651e+00
## 8   35 - 39                           Cancer morbidity yearly  1.353225e+01
## 9   40 - 44                           Cancer morbidity yearly  1.957787e+01
## 10  45 - 49                           Cancer morbidity yearly  3.046821e+01
## 11  50 - 54                           Cancer morbidity yearly  5.636304e+01
## 12  55 - 59                           Cancer morbidity yearly  9.222244e+01
## 13  60 - 64                           Cancer morbidity yearly  1.342997e+02
## 14  65 - 69                           Cancer morbidity yearly  3.203521e+02
## 15  70 - 74                           Cancer morbidity yearly  3.754101e+02
## 16  75 - 79                           Cancer morbidity yearly  2.612220e+02
## 17  80 - 84                           Cancer morbidity yearly  2.111356e+02
## 18      85+                           Cancer morbidity yearly  1.694984e+02
## 19    0 - 4 Dioxin recommendation tolerable daily intake 2018  1.307113e+02
## 20    5 - 9 Dioxin recommendation tolerable daily intake 2018  7.752637e+01
## 21  10 - 14 Dioxin recommendation tolerable daily intake 2018  8.403418e+01
## 22  15 - 19 Dioxin recommendation tolerable daily intake 2018  8.072697e+01
## 23  20 - 24 Dioxin recommendation tolerable daily intake 2018  7.946513e+01
## 24  25 - 29 Dioxin recommendation tolerable daily intake 2018  9.003056e+01
## 25  30 - 34 Dioxin recommendation tolerable daily intake 2018  8.878598e+01
## 26  35 - 39 Dioxin recommendation tolerable daily intake 2018  9.132071e+01
## 27  40 - 44 Dioxin recommendation tolerable daily intake 2018  8.769934e+01
## 28  45 - 49 Dioxin recommendation tolerable daily intake 2018  7.957372e+01
## 29  50 - 54 Dioxin recommendation tolerable daily intake 2018  9.005815e+01
## 30  55 - 59 Dioxin recommendation tolerable daily intake 2018  9.336819e+01
## 31  60 - 64 Dioxin recommendation tolerable daily intake 2018  9.102633e+01
## 32  65 - 69 Dioxin recommendation tolerable daily intake 2018  1.034493e+02
## 33  70 - 74 Dioxin recommendation tolerable daily intake 2018  1.013543e+02
## 34  75 - 79 Dioxin recommendation tolerable daily intake 2018  6.149027e+01
## 35  80 - 84 Dioxin recommendation tolerable daily intake 2018  4.677361e+01
## 36      85+ Dioxin recommendation tolerable daily intake 2018  4.376368e+01
## 37    0 - 4                                 Immunosuppression  1.409758e+00
## 38    5 - 9                                 Immunosuppression  3.950077e+00
## 39    0 - 4                         Loss in child's IQ points  1.169670e+01
## 40    0 - 4                                          PFAS TWI  8.396932e+00
## 41    5 - 9                                          PFAS TWI  6.689257e+01
## 42  10 - 14                                          PFAS TWI  6.477286e+01
## 43  15 - 19                                          PFAS TWI  6.229605e+01
## 44  20 - 24                                          PFAS TWI  9.059494e+01
## 45  25 - 29                                          PFAS TWI  1.026361e+02
## 46  30 - 34                                          PFAS TWI  1.012443e+02
## 47  35 - 39                                          PFAS TWI  1.041267e+02
## 48  40 - 44                                          PFAS TWI  9.995775e+01
## 49  45 - 49                                          PFAS TWI  9.062490e+01
## 50  50 - 54                                          PFAS TWI  1.024735e+02
## 51  55 - 59                                          PFAS TWI  1.061396e+02
## 52  60 - 64                                          PFAS TWI  1.033391e+02
## 53  65 - 69                                          PFAS TWI  1.090532e+02
## 54  70 - 74                                          PFAS TWI  1.066389e+02
## 55  75 - 79                                          PFAS TWI  6.442962e+01
## 56  80 - 84                                          PFAS TWI  4.869127e+01
## 57      85+                                          PFAS TWI  4.488841e+01
## 58    0 - 4                               Sperm concentration  3.673540e+02
## 59    0 - 4                          Vitamin D recommendation  2.789953e+02
## 60    5 - 9                          Vitamin D recommendation  2.954665e+02
## 61  10 - 14                          Vitamin D recommendation  3.091448e+02
## 62  15 - 19                          Vitamin D recommendation  2.968834e+02
## 63  20 - 24                          Vitamin D recommendation  2.836553e+02
## 64  25 - 29                          Vitamin D recommendation  3.213781e+02
## 65  30 - 34                          Vitamin D recommendation  3.168763e+02
## 66  35 - 39                          Vitamin D recommendation  3.259402e+02
## 67  40 - 44                          Vitamin D recommendation  3.131014e+02
## 68  45 - 49                          Vitamin D recommendation  2.842471e+02
## 69  50 - 54                          Vitamin D recommendation  3.218989e+02
## 70  55 - 59                          Vitamin D recommendation  3.339483e+02
## 71  60 - 64                          Vitamin D recommendation  3.258734e+02
## 72  65 - 69                          Vitamin D recommendation  3.037856e+02
## 73  70 - 74                          Vitamin D recommendation  2.979933e+02
## 74  75 - 79                          Vitamin D recommendation  1.812547e+02
## 75  80 - 84                          Vitamin D recommendation  1.384307e+02
## 76      85+                          Vitamin D recommendation  1.306929e+02
## 77    0 - 4                           Yes or no dental defect  2.255433e+02
## 78    0 - 4                               All-cause mortality -6.997371e+01
## 79    5 - 9                               All-cause mortality -3.839834e+01
## 80  10 - 14                               All-cause mortality -4.622956e+01
## 81  15 - 19                               All-cause mortality -1.309934e+02
## 82  20 - 24                               All-cause mortality -3.452314e+02
## 83  25 - 29                               All-cause mortality -4.119663e+02
## 84  30 - 34                               All-cause mortality -4.379799e+02
## 85  35 - 39                               All-cause mortality -5.407760e+02
## 86  40 - 44                               All-cause mortality -6.646000e+02
## 87  45 - 49                               All-cause mortality -8.853071e+02
## 88  50 - 54                               All-cause mortality -1.371145e+03
## 89  55 - 59                               All-cause mortality -1.921390e+03
## 90  60 - 64                               All-cause mortality -2.522767e+03
## 91  65 - 69                               All-cause mortality -3.921957e+03
## 92  70 - 74                               All-cause mortality -4.724217e+03
## 93  75 - 79                               All-cause mortality -3.859363e+03
## 94  80 - 84                               All-cause mortality -3.975712e+03
## 95      85+                               All-cause mortality -5.154125e+03
## 96    0 - 4                                        Depression -1.494882e-02
## 97    5 - 9                                        Depression -7.312900e+00
## 98  10 - 14                                        Depression -5.496555e+01
## 99  15 - 19                                        Depression -1.242428e+02
## 100 20 - 24                                        Depression -2.309439e+02
## 101 25 - 29                                        Depression -2.308706e+02
## 102 30 - 34                                        Depression -2.025407e+02
## 103 35 - 39                                        Depression -2.161037e+02
## 104 40 - 44                                        Depression -2.149180e+02
## 105 45 - 49                                        Depression -1.955520e+02
## 106 50 - 54                                        Depression -2.149183e+02
## 107 55 - 59                                        Depression -2.217746e+02
## 108 60 - 64                                        Depression -2.124055e+02
## 109 65 - 69                                        Depression -2.425608e+02
## 110 70 - 74                                        Depression -2.197377e+02
## 111 75 - 79                                        Depression -1.236256e+02
## 112 80 - 84                                        Depression -9.091906e+01
## 113     85+                                        Depression -8.467332e+01
## 114   0 - 4                                    CHD2 mortality -6.018588e+00
## 115   5 - 9                                    CHD2 mortality -5.962606e+00
## 116 10 - 14                                    CHD2 mortality -9.536885e+00
## 117 15 - 19                                    CHD2 mortality -1.849165e+01
## 118 20 - 24                                    CHD2 mortality -3.079494e+01
## 119 25 - 29                                    CHD2 mortality -5.345880e+01
## 120 30 - 34                                    CHD2 mortality -7.779232e+01
## 121 35 - 39                                    CHD2 mortality -1.432076e+02
## 122 40 - 44                                    CHD2 mortality -2.428114e+02
## 123 45 - 49                                    CHD2 mortality -3.823758e+02
## 124 50 - 54                                    CHD2 mortality -6.728882e+02
## 125 55 - 59                                    CHD2 mortality -1.075147e+03
## 126 60 - 64                                    CHD2 mortality -1.604622e+03
## 127 65 - 69                                    CHD2 mortality -2.257671e+03
## 128 70 - 74                                    CHD2 mortality -3.050607e+03
## 129 75 - 79                                    CHD2 mortality -2.895608e+03
## 130 80 - 84                                    CHD2 mortality -3.438948e+03
## 131     85+                                    CHD2 mortality -5.113855e+03
## 132 15 - 19                                     Breast cancer -2.863128e-01
## 133 20 - 24                                     Breast cancer -1.245421e+00
## 134 25 - 29                                     Breast cancer -5.731800e+00
## 135 30 - 34                                     Breast cancer -2.292193e+01
## 136 35 - 39                                     Breast cancer -4.144614e+01
## 137 40 - 44                                     Breast cancer -6.708729e+01
## 138 45 - 49                                     Breast cancer -1.095386e+02
## 139 50 - 54                                     Breast cancer -1.638916e+02
## 140 55 - 59                                     Breast cancer -2.103540e+02
## 141 60 - 64                                     Breast cancer -2.387390e+02
## 142 65 - 69                                     Breast cancer -3.208953e+02
## 143 70 - 74                                     Breast cancer -3.233861e+02
## 144 75 - 79                                     Breast cancer -2.393205e+02
## 145 80 - 84                                     Breast cancer -1.836715e+02
## 146     85+                                     Breast cancer -1.697011e+02
oprint(summary(BoDattr,marginals=c("Exposure_agent","Response"),"mean"))
##    Exposure_agent                                          Response
## 1             TEQ                           Cancer morbidity yearly
## 2             TEQ Dioxin recommendation tolerable daily intake 2018
## 3            PFAS                                 Immunosuppression
## 4             DHA                         Loss in child's IQ points
## 5            MeHg                         Loss in child's IQ points
## 6            PFAS                                          PFAS TWI
## 7             TEQ                               Sperm concentration
## 8       Vitamin D                          Vitamin D recommendation
## 9             TEQ                           Yes or no dental defect
## 10           Fish                               All-cause mortality
## 11           Fish                                        Depression
## 12         Omega3                                    CHD2 mortality
## 13         Omega3                                     Breast cancer
##            mean
## 1     94.952474
## 2     84.508784
## 3      2.679917
## 4   -248.500988
## 5    271.894397
## 6     82.066488
## 7    367.354018
## 8    281.087002
## 9    225.543298
## 10 -1723.451728
## 11  -160.448890
## 12 -1171.099826
## 13  -139.881112
oprint(summary(BoDattr,marginals=c("Gender","Response"),"mean"))
##    Gender                                          Response         mean
## 1  Female                           Cancer morbidity yearly    72.216623
## 2    Male                           Cancer morbidity yearly   117.688326
## 3  Female Dioxin recommendation tolerable daily intake 2018    85.655597
## 4    Male Dioxin recommendation tolerable daily intake 2018    83.361970
## 5  Female                                 Immunosuppression     2.174741
## 6    Male                                 Immunosuppression     3.185094
## 7  Female                         Loss in child's IQ points    10.895974
## 8    Male                         Loss in child's IQ points    12.497435
## 9  Female                                          PFAS TWI    75.169428
## 10   Male                                          PFAS TWI    88.963549
## 11 Female                               Sperm concentration   327.716920
## 12   Male                               Sperm concentration   406.991116
## 13 Female                          Vitamin D recommendation   299.605848
## 14   Male                          Vitamin D recommendation   262.568157
## 15 Female                           Yes or no dental defect   199.063593
## 16   Male                           Yes or no dental defect   252.023003
## 17 Female                               All-cause mortality -1218.527576
## 18   Male                               All-cause mortality -2228.375881
## 19 Female                                        Depression  -177.905392
## 20   Male                                        Depression  -142.992388
## 21 Female                                    CHD2 mortality  -958.507116
## 22   Male                                    CHD2 mortality -1383.692537
## 23 Female                                     Breast cancer  -277.673762
## 24   Male                                     Breast cancer    -2.088462
oprint(summary(case_burden,"mean"))
##                                            Response        mean
## 1      Dioxin recommendation tolerable daily intake 0.002193848
## 2 Dioxin recommendation tolerable daily intake 2018 0.002055036
## 3                         Loss in child's IQ points 0.111114237
## 4                                          PFAS TWI 0.002209610
## 5                               Sperm concentration 2.472082715
## 6                          Vitamin D recommendation 0.002203746
## 7                           Yes or no dental defect 0.060105744
oprint(summary(conc,"mean"))
##                Fish  Compound        mean
## 1      Average fish       ALA  0.69000000
## 2             Bream       ALA  0.22000000
## 3           Herring       ALA  1.74000000
## 4              Pike       ALA  0.08000000
## 5     Rainbow trout       ALA  4.81000000
## 6             Roach       ALA  0.10000000
## 7            Salmon       ALA  7.96000000
## 8           Vendace       ALA  1.35000000
## 9         Whitefish       ALA  2.22000000
## 10     Average fish       DHA  2.54000000
## 11            Bream       DHA  2.73000000
## 12          Herring       DHA  5.86000000
## 13             Pike       DHA  0.30000000
## 14    Rainbow trout       DHA  7.57000000
## 15            Roach       DHA  2.87000000
## 16           Salmon       DHA  6.69000000
## 17          Vendace       DHA  3.00000000
## 18        Whitefish       DHA  3.94000000
## 19     Average fish      Fish  1.00000000
## 20            Bream      Fish  1.00000000
## 21          Herring      Fish  1.00000000
## 22             Pike      Fish  1.00000000
## 23    Rainbow trout      Fish  1.00000000
## 24            Roach      Fish  1.00000000
## 25           Salmon      Fish  1.00000000
## 26          Vendace      Fish  1.00000000
## 27        Whitefish      Fish  1.00000000
## 28     Average fish      MeHg  0.04173655
## 29            Bream      MeHg  0.74219290
## 30           Burbot      MeHg  1.28856628
## 31         Flounder      MeHg  0.05654086
## 32 Fourhorn sculpin      MeHg  0.68213993
## 33          Herring      MeHg  0.03006897
## 34            Perch      MeHg  0.21755263
## 35             Pike      MeHg  1.39428894
## 36       Pike-perch      MeHg  0.35321728
## 37            Roach      MeHg  0.32359742
## 38           Salmon      MeHg  0.17732483
## 39            Sprat      MeHg  0.04861800
## 40      White bream      MeHg  1.37967806
## 41        Whitefish      MeHg  0.07528727
## 42     Average fish    Omega3  7.00000000
## 43            Bream    Omega3  6.00000000
## 44          Herring    Omega3 24.00000000
## 45             Pike    Omega3  0.48755001
## 46    Rainbow trout    Omega3 18.00000000
## 47            Roach    Omega3  5.00000000
## 48           Salmon    Omega3 23.00000000
## 49          Vendace    Omega3 10.00000000
## 50        Whitefish    Omega3 10.00000000
## 51     Average fish      PFAS  8.89745683
## 52            Bream      PFAS  4.94701243
## 53              Eel      PFAS  7.90154169
## 54          Herring      PFAS  1.65501151
## 55            Perch      PFAS  8.99982753
## 56             Pike      PFAS  8.74241994
## 57       Pike-perch      PFAS  2.59682481
## 58    Rainbow trout      PFAS  8.21572049
## 59            Roach      PFAS  9.21044812
## 60           Salmon      PFAS  8.63064676
## 61          Vendace      PFAS  8.68377750
## 62        Whitefish      PFAS  8.57080660
## 63            Bream       TEQ  6.20175433
## 64           Burbot       TEQ  0.82549312
## 65              Cod       TEQ  1.67367981
## 66         Flounder       TEQ 15.36558024
## 67          Herring       TEQ 90.08722156
## 68            Perch       TEQ  3.56897420
## 69             Pike       TEQ  2.14862962
## 70       Pike-perch       TEQ  2.70814568
## 71    River lamprey       TEQ  7.34301416
## 72            Roach       TEQ  0.72124877
## 73           Salmon       TEQ 20.00778296
## 74        Sea trout       TEQ  9.31778459
## 75            Sprat       TEQ  4.15456559
## 76          Vendace       TEQ  0.99253888
## 77        Whitefish       TEQ  5.74967872
## 78     Average fish Vitamin D  0.10500000
## 79            Bream Vitamin D  0.14000000
## 80          Herring Vitamin D  0.15600000
## 81             Pike Vitamin D  0.02100000
## 82    Rainbow trout Vitamin D  0.05100000
## 83            Roach Vitamin D  0.10000000
## 84           Salmon Vitamin D  0.06700000
## 85          Vendace Vitamin D  0.09400000
## 86        Whitefish Vitamin D  0.14400000
oprint(summary(dose,"mean"))
##     Gender     Age Exposure_agent Exposure Scaling         mean
## 1   Female   0 - 4            ALA   Direct      BW 3.061958e-01
## 2     Male   0 - 4            ALA   Direct      BW 4.142649e-01
## 3   Female   5 - 9            ALA   Direct      BW 5.669513e-01
## 4     Male   5 - 9            ALA   Direct      BW 7.670517e-01
## 5   Female 10 - 14            ALA   Direct      BW 8.626566e-01
## 6     Male 10 - 14            ALA   Direct      BW 1.167124e+00
## 7   Female 15 - 19            ALA   Direct      BW 8.626566e-01
## 8     Male 15 - 19            ALA   Direct      BW 1.167124e+00
## 9   Female 20 - 24            ALA   Direct      BW 1.425036e+00
## 10    Male 20 - 24            ALA   Direct      BW 1.927990e+00
## 11  Female 25 - 29            ALA   Direct      BW 1.425036e+00
## 12    Male 25 - 29            ALA   Direct      BW 1.927990e+00
## 13  Female 30 - 34            ALA   Direct      BW 1.425036e+00
## 14    Male 30 - 34            ALA   Direct      BW 1.927990e+00
## 15  Female 35 - 39            ALA   Direct      BW 1.425036e+00
## 16    Male 35 - 39            ALA   Direct      BW 1.927990e+00
## 17  Female 40 - 44            ALA   Direct      BW 1.425036e+00
## 18    Male 40 - 44            ALA   Direct      BW 1.927990e+00
## 19  Female 45 - 49            ALA   Direct      BW 1.425036e+00
## 20    Male 45 - 49            ALA   Direct      BW 1.927990e+00
## 21  Female 50 - 54            ALA   Direct      BW 1.425036e+00
## 22    Male 50 - 54            ALA   Direct      BW 1.927990e+00
## 23  Female 55 - 59            ALA   Direct      BW 1.425036e+00
## 24    Male 55 - 59            ALA   Direct      BW 1.927990e+00
## 25  Female 60 - 64            ALA   Direct      BW 1.425036e+00
## 26    Male 60 - 64            ALA   Direct      BW 1.927990e+00
## 27  Female 65 - 69            ALA   Direct      BW 1.979800e+00
## 28    Male 65 - 69            ALA   Direct      BW 2.678553e+00
## 29  Female 70 - 74            ALA   Direct      BW 1.979800e+00
## 30    Male 70 - 74            ALA   Direct      BW 2.678553e+00
## 31  Female 75 - 79            ALA   Direct      BW 1.979800e+00
## 32    Male 75 - 79            ALA   Direct      BW 2.678553e+00
## 33  Female 80 - 84            ALA   Direct      BW 1.979800e+00
## 34    Male 80 - 84            ALA   Direct      BW 2.678553e+00
## 35  Female     85+            ALA   Direct      BW 1.979800e+00
## 36    Male     85+            ALA   Direct      BW 2.678553e+00
## 37  Female   0 - 4            DHA   Direct      BW 7.933049e-01
## 38    Male   0 - 4            DHA   Direct      BW 1.093282e+00
## 39  Female   5 - 9            DHA   Direct      BW 1.311580e+00
## 40    Male   5 - 9            DHA   Direct      BW 1.794477e+00
## 41  Female 10 - 14            DHA   Direct      BW 1.484387e+00
## 42    Male 10 - 14            DHA   Direct      BW 2.028275e+00
## 43  Female 15 - 19            DHA   Direct      BW 1.484387e+00
## 44    Male 15 - 19            DHA   Direct      BW 2.028275e+00
## 45  Female 20 - 24            DHA   Direct      BW 2.156782e+00
## 46    Male 20 - 24            DHA   Direct      BW 2.937986e+00
## 47  Female 25 - 29            DHA   Direct      BW 2.156782e+00
## 48    Male 25 - 29            DHA   Direct      BW 2.937986e+00
## 49  Female 30 - 34            DHA   Direct      BW 2.156782e+00
## 50    Male 30 - 34            DHA   Direct      BW 2.937986e+00
## 51  Female 35 - 39            DHA   Direct      BW 2.156782e+00
## 52    Male 35 - 39            DHA   Direct      BW 2.937986e+00
## 53  Female 40 - 44            DHA   Direct      BW 2.156782e+00
## 54    Male 40 - 44            DHA   Direct      BW 2.937986e+00
## 55  Female 45 - 49            DHA   Direct      BW 2.156782e+00
## 56    Male 45 - 49            DHA   Direct      BW 2.937986e+00
## 57  Female 50 - 54            DHA   Direct      BW 2.156782e+00
## 58    Male 50 - 54            DHA   Direct      BW 2.937986e+00
## 59  Female 55 - 59            DHA   Direct      BW 2.156782e+00
## 60    Male 55 - 59            DHA   Direct      BW 2.937986e+00
## 61  Female 60 - 64            DHA   Direct      BW 2.156782e+00
## 62    Male 60 - 64            DHA   Direct      BW 2.937986e+00
## 63  Female 65 - 69            DHA   Direct      BW 2.784172e+00
## 64    Male 65 - 69            DHA   Direct      BW 3.786808e+00
## 65  Female 70 - 74            DHA   Direct      BW 2.784172e+00
## 66    Male 70 - 74            DHA   Direct      BW 3.786808e+00
## 67  Female 75 - 79            DHA   Direct      BW 2.784172e+00
## 68    Male 75 - 79            DHA   Direct      BW 3.786808e+00
## 69  Female 80 - 84            DHA   Direct      BW 2.784172e+00
## 70    Male 80 - 84            DHA   Direct      BW 3.786808e+00
## 71  Female     85+            DHA   Direct      BW 2.784172e+00
## 72    Male     85+            DHA   Direct      BW 3.786808e+00
## 73  Female   0 - 4            EPA   Direct      BW 4.244296e-01
## 74    Male   0 - 4            EPA   Direct      BW 6.210307e-01
## 75  Female   5 - 9            EPA   Direct      BW 4.244296e-01
## 76    Male   5 - 9            EPA   Direct      BW 6.210307e-01
## 77  Female 10 - 14            EPA   Direct      BW 4.244296e-01
## 78    Male 10 - 14            EPA   Direct      BW 6.210307e-01
## 79  Female 15 - 19            EPA   Direct      BW 4.244296e-01
## 80    Male 15 - 19            EPA   Direct      BW 6.210307e-01
## 81  Female 20 - 24            EPA   Direct      BW 4.244296e-01
## 82    Male 20 - 24            EPA   Direct      BW 6.210307e-01
## 83  Female 25 - 29            EPA   Direct      BW 4.244296e-01
## 84    Male 25 - 29            EPA   Direct      BW 6.210307e-01
## 85  Female 30 - 34            EPA   Direct      BW 4.244296e-01
## 86    Male 30 - 34            EPA   Direct      BW 6.210307e-01
## 87  Female 35 - 39            EPA   Direct      BW 4.244296e-01
## 88    Male 35 - 39            EPA   Direct      BW 6.210307e-01
## 89  Female 40 - 44            EPA   Direct      BW 4.244296e-01
## 90    Male 40 - 44            EPA   Direct      BW 6.210307e-01
## 91  Female 45 - 49            EPA   Direct      BW 4.244296e-01
## 92    Male 45 - 49            EPA   Direct      BW 6.210307e-01
## 93  Female 50 - 54            EPA   Direct      BW 4.244296e-01
## 94    Male 50 - 54            EPA   Direct      BW 6.210307e-01
## 95  Female 55 - 59            EPA   Direct      BW 4.244296e-01
## 96    Male 55 - 59            EPA   Direct      BW 6.210307e-01
## 97  Female 60 - 64            EPA   Direct      BW 4.244296e-01
## 98    Male 60 - 64            EPA   Direct      BW 6.210307e-01
## 99  Female 65 - 69            EPA   Direct      BW 4.244296e-01
## 100   Male 65 - 69            EPA   Direct      BW 6.210307e-01
## 101 Female 70 - 74            EPA   Direct      BW 4.244296e-01
## 102   Male 70 - 74            EPA   Direct      BW 6.210307e-01
## 103 Female 75 - 79            EPA   Direct      BW 4.244296e-01
## 104   Male 75 - 79            EPA   Direct      BW 6.210307e-01
## 105 Female 80 - 84            EPA   Direct      BW 4.244296e-01
## 106   Male 80 - 84            EPA   Direct      BW 6.210307e-01
## 107 Female     85+            EPA   Direct      BW 4.244296e-01
## 108   Male     85+            EPA   Direct      BW 6.210307e-01
## 109 Female   0 - 4           Fish   Direct      BW 8.231686e-02
## 110   Male   0 - 4           Fish   Direct      BW 1.113699e-01
## 111 Female   5 - 9           Fish   Direct      BW 2.509984e-01
## 112   Male   5 - 9           Fish   Direct      BW 3.395860e-01
## 113 Female 10 - 14           Fish   Direct      BW 2.516307e-01
## 114   Male 10 - 14           Fish   Direct      BW 3.404415e-01
## 115 Female 15 - 19           Fish   Direct      BW 2.516307e-01
## 116   Male 15 - 19           Fish   Direct      BW 3.404415e-01
## 117 Female 20 - 24           Fish   Direct      BW 3.952884e-01
## 118   Male 20 - 24           Fish   Direct      BW 5.348020e-01
## 119 Female 25 - 29           Fish   Direct      BW 3.952884e-01
## 120   Male 25 - 29           Fish   Direct      BW 5.348020e-01
## 121 Female 30 - 34           Fish   Direct      BW 3.952884e-01
## 122   Male 30 - 34           Fish   Direct      BW 5.348020e-01
## 123 Female 35 - 39           Fish   Direct      BW 3.952884e-01
## 124   Male 35 - 39           Fish   Direct      BW 5.348020e-01
## 125 Female 40 - 44           Fish   Direct      BW 3.952884e-01
## 126   Male 40 - 44           Fish   Direct      BW 5.348020e-01
## 127 Female 45 - 49           Fish   Direct      BW 3.952884e-01
## 128   Male 45 - 49           Fish   Direct      BW 5.348020e-01
## 129 Female 50 - 54           Fish   Direct      BW 3.952884e-01
## 130   Male 50 - 54           Fish   Direct      BW 5.348020e-01
## 131 Female 55 - 59           Fish   Direct      BW 3.952884e-01
## 132   Male 55 - 59           Fish   Direct      BW 5.348020e-01
## 133 Female 60 - 64           Fish   Direct      BW 3.952884e-01
## 134   Male 60 - 64           Fish   Direct      BW 5.348020e-01
## 135 Female 65 - 69           Fish   Direct      BW 5.110244e-01
## 136   Male 65 - 69           Fish   Direct      BW 6.913859e-01
## 137 Female 70 - 74           Fish   Direct      BW 5.110244e-01
## 138   Male 70 - 74           Fish   Direct      BW 6.913859e-01
## 139 Female 75 - 79           Fish   Direct      BW 5.110244e-01
## 140   Male 75 - 79           Fish   Direct      BW 6.913859e-01
## 141 Female 80 - 84           Fish   Direct      BW 5.110244e-01
## 142   Male 80 - 84           Fish   Direct      BW 6.913859e-01
## 143 Female     85+           Fish   Direct      BW 5.110244e-01
## 144   Male     85+           Fish   Direct      BW 6.913859e-01
## 145 Female   0 - 4           MeHg   Direct      BW 8.531094e-03
## 146   Male   0 - 4           MeHg   Direct      BW 1.154207e-02
## 147 Female   5 - 9           MeHg   Direct      BW 1.895836e-02
## 148   Male   5 - 9           MeHg   Direct      BW 2.564955e-02
## 149 Female 10 - 14           MeHg   Direct      BW 2.330731e-02
## 150   Male 10 - 14           MeHg   Direct      BW 3.153343e-02
## 151 Female 15 - 19           MeHg   Direct      BW 2.330731e-02
## 152   Male 15 - 19           MeHg   Direct      BW 3.153343e-02
## 153 Female 20 - 24           MeHg   Direct      BW 4.040592e-02
## 154   Male 20 - 24           MeHg   Direct      BW 5.466683e-02
## 155 Female 25 - 29           MeHg   Direct      BW 4.040592e-02
## 156   Male 25 - 29           MeHg   Direct      BW 5.466683e-02
## 157 Female 30 - 34           MeHg   Direct      BW 4.040592e-02
## 158   Male 30 - 34           MeHg   Direct      BW 5.466683e-02
## 159 Female 35 - 39           MeHg   Direct      BW 4.040592e-02
## 160   Male 35 - 39           MeHg   Direct      BW 5.466683e-02
## 161 Female 40 - 44           MeHg   Direct      BW 4.040592e-02
## 162   Male 40 - 44           MeHg   Direct      BW 5.466683e-02
## 163 Female 45 - 49           MeHg   Direct      BW 4.040592e-02
## 164   Male 45 - 49           MeHg   Direct      BW 5.466683e-02
## 165 Female 50 - 54           MeHg   Direct      BW 4.040592e-02
## 166   Male 50 - 54           MeHg   Direct      BW 5.466683e-02
## 167 Female 55 - 59           MeHg   Direct      BW 4.040592e-02
## 168   Male 55 - 59           MeHg   Direct      BW 5.466683e-02
## 169 Female 60 - 64           MeHg   Direct      BW 4.040592e-02
## 170   Male 60 - 64           MeHg   Direct      BW 5.466683e-02
## 171 Female 65 - 69           MeHg   Direct      BW 5.346493e-02
## 172   Male 65 - 69           MeHg   Direct      BW 7.233491e-02
## 173 Female 70 - 74           MeHg   Direct      BW 5.346493e-02
## 174   Male 70 - 74           MeHg   Direct      BW 7.233491e-02
## 175 Female 75 - 79           MeHg   Direct      BW 5.346493e-02
## 176   Male 75 - 79           MeHg   Direct      BW 7.233491e-02
## 177 Female 80 - 84           MeHg   Direct      BW 5.346493e-02
## 178   Male 80 - 84           MeHg   Direct      BW 7.233491e-02
## 179 Female     85+           MeHg   Direct      BW 5.346493e-02
## 180   Male     85+           MeHg   Direct      BW 7.233491e-02
## 181 Female   0 - 4         Omega3   Direct      BW 2.020964e+00
## 182   Male   0 - 4         Omega3   Direct      BW 2.753231e+00
## 183 Female   5 - 9         Omega3   Direct      BW 3.554566e+00
## 184   Male   5 - 9         Omega3   Direct      BW 4.828105e+00
## 185 Female 10 - 14         Omega3   Direct      BW 4.223030e+00
## 186   Male 10 - 14         Omega3   Direct      BW 5.732497e+00
## 187 Female 15 - 19         Omega3   Direct      BW 4.223030e+00
## 188   Male 15 - 19         Omega3   Direct      BW 5.732497e+00
## 189 Female 20 - 24         Omega3   Direct      BW 6.341576e+00
## 190   Male 20 - 24         Omega3   Direct      BW 8.598765e+00
## 191 Female 25 - 29         Omega3   Direct      BW 6.341576e+00
## 192   Male 25 - 29         Omega3   Direct      BW 8.598765e+00
## 193 Female 30 - 34         Omega3   Direct      BW 6.341576e+00
## 194   Male 30 - 34         Omega3   Direct      BW 8.598765e+00
## 195 Female 35 - 39         Omega3   Direct      BW 6.341576e+00
## 196   Male 35 - 39         Omega3   Direct      BW 8.598765e+00
## 197 Female 40 - 44         Omega3   Direct      BW 6.341576e+00
## 198   Male 40 - 44         Omega3   Direct      BW 8.598765e+00
## 199 Female 45 - 49         Omega3   Direct      BW 6.341576e+00
## 200   Male 45 - 49         Omega3   Direct      BW 8.598765e+00
## 201 Female 50 - 54         Omega3   Direct      BW 6.341576e+00
## 202   Male 50 - 54         Omega3   Direct      BW 8.598765e+00
## 203 Female 55 - 59         Omega3   Direct      BW 6.341576e+00
## 204   Male 55 - 59         Omega3   Direct      BW 8.598765e+00
## 205 Female 60 - 64         Omega3   Direct      BW 6.341576e+00
## 206   Male 60 - 64         Omega3   Direct      BW 8.598765e+00
## 207 Female 65 - 69         Omega3   Direct      BW 8.559058e+00
## 208   Male 65 - 69         Omega3   Direct      BW 1.159889e+01
## 209 Female 70 - 74         Omega3   Direct      BW 8.559058e+00
## 210   Male 70 - 74         Omega3   Direct      BW 1.159889e+01
## 211 Female 75 - 79         Omega3   Direct      BW 8.559058e+00
## 212   Male 75 - 79         Omega3   Direct      BW 1.159889e+01
## 213 Female 80 - 84         Omega3   Direct      BW 8.559058e+00
## 214   Male 80 - 84         Omega3   Direct      BW 1.159889e+01
## 215 Female     85+         Omega3   Direct      BW 8.559058e+00
## 216   Male     85+         Omega3   Direct      BW 1.159889e+01
## 217 Female   0 - 4           PFAS   Direct      BW 6.749939e-01
## 218   Male   0 - 4           PFAS   Direct      BW 9.132271e-01
## 219 Female   5 - 9           PFAS   Direct      BW 2.135159e+00
## 220   Male   5 - 9           PFAS   Direct      BW 2.888744e+00
## 221 Female 10 - 14           PFAS   Direct      BW 2.167332e+00
## 222   Male 10 - 14           PFAS   Direct      BW 2.932272e+00
## 223 Female 15 - 19           PFAS   Direct      BW 2.167332e+00
## 224   Male 15 - 19           PFAS   Direct      BW 2.932272e+00
## 225 Female 20 - 24           PFAS   Direct      BW 3.296255e+00
## 226   Male 20 - 24           PFAS   Direct      BW 4.459638e+00
## 227 Female 25 - 29           PFAS   Direct      BW 3.296255e+00
## 228   Male 25 - 29           PFAS   Direct      BW 4.459638e+00
## 229 Female 30 - 34           PFAS   Direct      BW 3.296255e+00
## 230   Male 30 - 34           PFAS   Direct      BW 4.459638e+00
## 231 Female 35 - 39           PFAS   Direct      BW 3.296255e+00
## 232   Male 35 - 39           PFAS   Direct      BW 4.459638e+00
## 233 Female 40 - 44           PFAS   Direct      BW 3.296255e+00
## 234   Male 40 - 44           PFAS   Direct      BW 4.459638e+00
## 235 Female 45 - 49           PFAS   Direct      BW 3.296255e+00
## 236   Male 45 - 49           PFAS   Direct      BW 4.459638e+00
## 237 Female 50 - 54           PFAS   Direct      BW 3.296255e+00
## 238   Male 50 - 54           PFAS   Direct      BW 4.459638e+00
## 239 Female 55 - 59           PFAS   Direct      BW 3.296255e+00
## 240   Male 55 - 59           PFAS   Direct      BW 4.459638e+00
## 241 Female 60 - 64           PFAS   Direct      BW 3.296255e+00
## 242   Male 60 - 64           PFAS   Direct      BW 4.459638e+00
## 243 Female 65 - 69           PFAS   Direct      BW 4.410885e+00
## 244   Male 65 - 69           PFAS   Direct      BW 5.967668e+00
## 245 Female 70 - 74           PFAS   Direct      BW 4.410885e+00
## 246   Male 70 - 74           PFAS   Direct      BW 5.967668e+00
## 247 Female 75 - 79           PFAS   Direct      BW 4.410885e+00
## 248   Male 75 - 79           PFAS   Direct      BW 5.967668e+00
## 249 Female 80 - 84           PFAS   Direct      BW 4.410885e+00
## 250   Male 80 - 84           PFAS   Direct      BW 5.967668e+00
## 251 Female     85+           PFAS   Direct      BW 4.410885e+00
## 252   Male     85+           PFAS   Direct      BW 5.967668e+00
## 253 Female   0 - 4            TEQ   Direct      BW 7.750475e-01
## 254   Male   0 - 4            TEQ   Direct      BW 1.048594e+00
## 255 Female   5 - 9            TEQ   Direct      BW 1.397440e+00
## 256   Male   5 - 9            TEQ   Direct      BW 1.890655e+00
## 257 Female 10 - 14            TEQ   Direct      BW 2.288415e+00
## 258   Male 10 - 14            TEQ   Direct      BW 3.096091e+00
## 259 Female 15 - 19            TEQ   Direct      BW 2.288415e+00
## 260   Male 15 - 19            TEQ   Direct      BW 3.096091e+00
## 261 Female 20 - 24            TEQ   Direct      BW 4.312746e+00
## 262   Male 20 - 24            TEQ   Direct      BW 5.834892e+00
## 263 Female 25 - 29            TEQ   Direct      BW 4.312746e+00
## 264   Male 25 - 29            TEQ   Direct      BW 5.834892e+00
## 265 Female 30 - 34            TEQ   Direct      BW 4.312746e+00
## 266   Male 30 - 34            TEQ   Direct      BW 5.834892e+00
## 267 Female 35 - 39            TEQ   Direct      BW 4.312746e+00
## 268   Male 35 - 39            TEQ   Direct      BW 5.834892e+00
## 269 Female 40 - 44            TEQ   Direct      BW 4.312746e+00
## 270   Male 40 - 44            TEQ   Direct      BW 5.834892e+00
## 271 Female 45 - 49            TEQ   Direct      BW 4.312746e+00
## 272   Male 45 - 49            TEQ   Direct      BW 5.834892e+00
## 273 Female 50 - 54            TEQ   Direct      BW 4.312746e+00
## 274   Male 50 - 54            TEQ   Direct      BW 5.834892e+00
## 275 Female 55 - 59            TEQ   Direct      BW 4.312746e+00
## 276   Male 55 - 59            TEQ   Direct      BW 5.834892e+00
## 277 Female 60 - 64            TEQ   Direct      BW 4.312746e+00
## 278   Male 60 - 64            TEQ   Direct      BW 5.834892e+00
## 279 Female 65 - 69            TEQ   Direct      BW 7.732477e+00
## 280   Male 65 - 69            TEQ   Direct      BW 1.046159e+01
## 281 Female 70 - 74            TEQ   Direct      BW 7.732477e+00
## 282   Male 70 - 74            TEQ   Direct      BW 1.046159e+01
## 283 Female 75 - 79            TEQ   Direct      BW 7.732477e+00
## 284   Male 75 - 79            TEQ   Direct      BW 1.046159e+01
## 285 Female 80 - 84            TEQ   Direct      BW 7.732477e+00
## 286   Male 80 - 84            TEQ   Direct      BW 1.046159e+01
## 287 Female     85+            TEQ   Direct      BW 7.732477e+00
## 288   Male     85+            TEQ   Direct      BW 1.046159e+01
## 289 Female   0 - 4      Vitamin D   Direct      BW 4.807978e-02
## 290   Male   0 - 4      Vitamin D   Direct      BW 7.419831e-02
## 291 Female   5 - 9      Vitamin D   Direct      BW 6.524467e-02
## 292   Male   5 - 9      Vitamin D   Direct      BW 9.742139e-02
## 293 Female 10 - 14      Vitamin D   Direct      BW 6.382524e-02
## 294   Male 10 - 14      Vitamin D   Direct      BW 9.550099e-02
## 295 Female 15 - 19      Vitamin D   Direct      BW 6.382524e-02
## 296   Male 15 - 19      Vitamin D   Direct      BW 9.550099e-02
## 297 Female 20 - 24      Vitamin D   Direct      BW 7.894433e-02
## 298   Male 20 - 24      Vitamin D   Direct      BW 1.159562e-01
## 299 Female 25 - 29      Vitamin D   Direct      BW 7.894433e-02
## 300   Male 25 - 29      Vitamin D   Direct      BW 1.159562e-01
## 301 Female 30 - 34      Vitamin D   Direct      BW 7.894433e-02
## 302   Male 30 - 34      Vitamin D   Direct      BW 1.159562e-01
## 303 Female 35 - 39      Vitamin D   Direct      BW 7.894433e-02
## 304   Male 35 - 39      Vitamin D   Direct      BW 1.159562e-01
## 305 Female 40 - 44      Vitamin D   Direct      BW 7.894433e-02
## 306   Male 40 - 44      Vitamin D   Direct      BW 1.159562e-01
## 307 Female 45 - 49      Vitamin D   Direct      BW 7.894433e-02
## 308   Male 45 - 49      Vitamin D   Direct      BW 1.159562e-01
## 309 Female 50 - 54      Vitamin D   Direct      BW 7.894433e-02
## 310   Male 50 - 54      Vitamin D   Direct      BW 1.159562e-01
## 311 Female 55 - 59      Vitamin D   Direct      BW 7.894433e-02
## 312   Male 55 - 59      Vitamin D   Direct      BW 1.159562e-01
## 313 Female 60 - 64      Vitamin D   Direct      BW 7.894433e-02
## 314   Male 60 - 64      Vitamin D   Direct      BW 1.159562e-01
## 315 Female 65 - 69      Vitamin D   Direct      BW 8.961582e-02
## 316   Male 65 - 69      Vitamin D   Direct      BW 1.303941e-01
## 317 Female 70 - 74      Vitamin D   Direct      BW 8.961582e-02
## 318   Male 70 - 74      Vitamin D   Direct      BW 1.303941e-01
## 319 Female 75 - 79      Vitamin D   Direct      BW 8.961582e-02
## 320   Male 75 - 79      Vitamin D   Direct      BW 1.303941e-01
## 321 Female 80 - 84      Vitamin D   Direct      BW 8.961582e-02
## 322   Male 80 - 84      Vitamin D   Direct      BW 1.303941e-01
## 323 Female     85+      Vitamin D   Direct      BW 8.961582e-02
## 324   Male     85+      Vitamin D   Direct      BW 1.303941e-01
## 325 Female   0 - 4           PFAS To child      BW 1.923055e-01
## 326   Male   0 - 4           PFAS To child      BW 1.923055e-01
## 327 Female   0 - 4            TEQ To child      BW 5.770853e-01
## 328   Male   0 - 4            TEQ To child      BW 5.770853e-01
## 329 Female   0 - 4            ALA   Direct   Log10         -Inf
## 330   Male   0 - 4            ALA   Direct   Log10         -Inf
## 331 Female   5 - 9            ALA   Direct   Log10         -Inf
## 332   Male   5 - 9            ALA   Direct   Log10         -Inf
## 333 Female 10 - 14            ALA   Direct   Log10         -Inf
## 334   Male 10 - 14            ALA   Direct   Log10         -Inf
## 335 Female 15 - 19            ALA   Direct   Log10         -Inf
## 336   Male 15 - 19            ALA   Direct   Log10         -Inf
## 337 Female 20 - 24            ALA   Direct   Log10         -Inf
## 338   Male 20 - 24            ALA   Direct   Log10         -Inf
## 339 Female 25 - 29            ALA   Direct   Log10         -Inf
## 340   Male 25 - 29            ALA   Direct   Log10         -Inf
## 341 Female 30 - 34            ALA   Direct   Log10         -Inf
## 342   Male 30 - 34            ALA   Direct   Log10         -Inf
## 343 Female 35 - 39            ALA   Direct   Log10         -Inf
## 344   Male 35 - 39            ALA   Direct   Log10         -Inf
## 345 Female 40 - 44            ALA   Direct   Log10         -Inf
## 346   Male 40 - 44            ALA   Direct   Log10         -Inf
## 347 Female 45 - 49            ALA   Direct   Log10         -Inf
## 348   Male 45 - 49            ALA   Direct   Log10         -Inf
## 349 Female 50 - 54            ALA   Direct   Log10         -Inf
## 350   Male 50 - 54            ALA   Direct   Log10         -Inf
## 351 Female 55 - 59            ALA   Direct   Log10         -Inf
## 352   Male 55 - 59            ALA   Direct   Log10         -Inf
## 353 Female 60 - 64            ALA   Direct   Log10         -Inf
## 354   Male 60 - 64            ALA   Direct   Log10         -Inf
## 355 Female 65 - 69            ALA   Direct   Log10         -Inf
## 356   Male 65 - 69            ALA   Direct   Log10         -Inf
## 357 Female 70 - 74            ALA   Direct   Log10         -Inf
## 358   Male 70 - 74            ALA   Direct   Log10         -Inf
## 359 Female 75 - 79            ALA   Direct   Log10         -Inf
## 360   Male 75 - 79            ALA   Direct   Log10         -Inf
## 361 Female 80 - 84            ALA   Direct   Log10         -Inf
## 362   Male 80 - 84            ALA   Direct   Log10         -Inf
## 363 Female     85+            ALA   Direct   Log10         -Inf
## 364   Male     85+            ALA   Direct   Log10         -Inf
## 365 Female   0 - 4            DHA   Direct   Log10 1.582628e+00
## 366   Male   0 - 4            DHA   Direct   Log10 1.722662e+00
## 367 Female   5 - 9            DHA   Direct   Log10 1.742764e+00
## 368   Male   5 - 9            DHA   Direct   Log10 1.896121e+00
## 369 Female 10 - 14            DHA   Direct   Log10 1.796361e+00
## 370   Male 10 - 14            DHA   Direct   Log10 1.935693e+00
## 371 Female 15 - 19            DHA   Direct   Log10 1.796361e+00
## 372   Male 15 - 19            DHA   Direct   Log10 1.935693e+00
## 373 Female 20 - 24            DHA   Direct   Log10 1.783013e+00
## 374   Male 20 - 24            DHA   Direct   Log10 1.936152e+00
## 375 Female 25 - 29            DHA   Direct   Log10 1.783013e+00
## 376   Male 25 - 29            DHA   Direct   Log10 1.936152e+00
## 377 Female 30 - 34            DHA   Direct   Log10 1.783013e+00
## 378   Male 30 - 34            DHA   Direct   Log10 1.936152e+00
## 379 Female 35 - 39            DHA   Direct   Log10 1.783013e+00
## 380   Male 35 - 39            DHA   Direct   Log10 1.936152e+00
## 381 Female 40 - 44            DHA   Direct   Log10 1.783013e+00
## 382   Male 40 - 44            DHA   Direct   Log10 1.936152e+00
## 383 Female 45 - 49            DHA   Direct   Log10 1.783013e+00
## 384   Male 45 - 49            DHA   Direct   Log10 1.936152e+00
## 385 Female 50 - 54            DHA   Direct   Log10 1.783013e+00
## 386   Male 50 - 54            DHA   Direct   Log10 1.936152e+00
## 387 Female 55 - 59            DHA   Direct   Log10 1.783013e+00
## 388   Male 55 - 59            DHA   Direct   Log10 1.936152e+00
## 389 Female 60 - 64            DHA   Direct   Log10 1.783013e+00
## 390   Male 60 - 64            DHA   Direct   Log10 1.936152e+00
## 391 Female 65 - 69            DHA   Direct   Log10 1.859742e+00
## 392   Male 65 - 69            DHA   Direct   Log10 1.993377e+00
## 393 Female 70 - 74            DHA   Direct   Log10 1.859742e+00
## 394   Male 70 - 74            DHA   Direct   Log10 1.993377e+00
## 395 Female 75 - 79            DHA   Direct   Log10 1.859742e+00
## 396   Male 75 - 79            DHA   Direct   Log10 1.993377e+00
## 397 Female 80 - 84            DHA   Direct   Log10 1.859742e+00
## 398   Male 80 - 84            DHA   Direct   Log10 1.993377e+00
## 399 Female     85+            DHA   Direct   Log10 1.859742e+00
## 400   Male     85+            DHA   Direct   Log10 1.993377e+00
## 401 Female   0 - 4            EPA   Direct   Log10 1.346994e+00
## 402   Male   0 - 4            EPA   Direct   Log10 1.505465e+00
## 403 Female   5 - 9            EPA   Direct   Log10 1.346994e+00
## 404   Male   5 - 9            EPA   Direct   Log10 1.505465e+00
## 405 Female 10 - 14            EPA   Direct   Log10 1.346994e+00
## 406   Male 10 - 14            EPA   Direct   Log10 1.505465e+00
## 407 Female 15 - 19            EPA   Direct   Log10 1.346994e+00
## 408   Male 15 - 19            EPA   Direct   Log10 1.505465e+00
## 409 Female 20 - 24            EPA   Direct   Log10 1.346994e+00
## 410   Male 20 - 24            EPA   Direct   Log10 1.505465e+00
## 411 Female 25 - 29            EPA   Direct   Log10 1.346994e+00
## 412   Male 25 - 29            EPA   Direct   Log10 1.505465e+00
## 413 Female 30 - 34            EPA   Direct   Log10 1.346994e+00
## 414   Male 30 - 34            EPA   Direct   Log10 1.505465e+00
## 415 Female 35 - 39            EPA   Direct   Log10 1.346994e+00
## 416   Male 35 - 39            EPA   Direct   Log10 1.505465e+00
## 417 Female 40 - 44            EPA   Direct   Log10 1.346994e+00
## 418   Male 40 - 44            EPA   Direct   Log10 1.505465e+00
## 419 Female 45 - 49            EPA   Direct   Log10 1.346994e+00
## 420   Male 45 - 49            EPA   Direct   Log10 1.505465e+00
## 421 Female 50 - 54            EPA   Direct   Log10 1.346994e+00
## 422   Male 50 - 54            EPA   Direct   Log10 1.505465e+00
## 423 Female 55 - 59            EPA   Direct   Log10 1.346994e+00
## 424   Male 55 - 59            EPA   Direct   Log10 1.505465e+00
## 425 Female 60 - 64            EPA   Direct   Log10 1.346994e+00
## 426   Male 60 - 64            EPA   Direct   Log10 1.505465e+00
## 427 Female 65 - 69            EPA   Direct   Log10 1.346994e+00
## 428   Male 65 - 69            EPA   Direct   Log10 1.505465e+00
## 429 Female 70 - 74            EPA   Direct   Log10 1.346994e+00
## 430   Male 70 - 74            EPA   Direct   Log10 1.505465e+00
## 431 Female 75 - 79            EPA   Direct   Log10 1.346994e+00
## 432   Male 75 - 79            EPA   Direct   Log10 1.505465e+00
## 433 Female 80 - 84            EPA   Direct   Log10 1.346994e+00
## 434   Male 80 - 84            EPA   Direct   Log10 1.505465e+00
## 435 Female     85+            EPA   Direct   Log10 1.346994e+00
## 436   Male     85+            EPA   Direct   Log10 1.505465e+00
## 437 Female   0 - 4           Fish   Direct   Log10         -Inf
## 438   Male   0 - 4           Fish   Direct   Log10         -Inf
## 439 Female   5 - 9           Fish   Direct   Log10         -Inf
## 440   Male   5 - 9           Fish   Direct   Log10         -Inf
## 441 Female 10 - 14           Fish   Direct   Log10         -Inf
## 442   Male 10 - 14           Fish   Direct   Log10         -Inf
## 443 Female 15 - 19           Fish   Direct   Log10         -Inf
## 444   Male 15 - 19           Fish   Direct   Log10         -Inf
## 445 Female 20 - 24           Fish   Direct   Log10         -Inf
## 446   Male 20 - 24           Fish   Direct   Log10         -Inf
## 447 Female 25 - 29           Fish   Direct   Log10         -Inf
## 448   Male 25 - 29           Fish   Direct   Log10         -Inf
## 449 Female 30 - 34           Fish   Direct   Log10         -Inf
## 450   Male 30 - 34           Fish   Direct   Log10         -Inf
## 451 Female 35 - 39           Fish   Direct   Log10         -Inf
## 452   Male 35 - 39           Fish   Direct   Log10         -Inf
## 453 Female 40 - 44           Fish   Direct   Log10         -Inf
## 454   Male 40 - 44           Fish   Direct   Log10         -Inf
## 455 Female 45 - 49           Fish   Direct   Log10         -Inf
## 456   Male 45 - 49           Fish   Direct   Log10         -Inf
## 457 Female 50 - 54           Fish   Direct   Log10         -Inf
## 458   Male 50 - 54           Fish   Direct   Log10         -Inf
## 459 Female 55 - 59           Fish   Direct   Log10         -Inf
## 460   Male 55 - 59           Fish   Direct   Log10         -Inf
## 461 Female 60 - 64           Fish   Direct   Log10         -Inf
## 462   Male 60 - 64           Fish   Direct   Log10         -Inf
## 463 Female 65 - 69           Fish   Direct   Log10         -Inf
## 464   Male 65 - 69           Fish   Direct   Log10         -Inf
## 465 Female 70 - 74           Fish   Direct   Log10         -Inf
## 466   Male 70 - 74           Fish   Direct   Log10         -Inf
## 467 Female 75 - 79           Fish   Direct   Log10         -Inf
## 468   Male 75 - 79           Fish   Direct   Log10         -Inf
## 469 Female 80 - 84           Fish   Direct   Log10         -Inf
## 470   Male 80 - 84           Fish   Direct   Log10         -Inf
## 471 Female     85+           Fish   Direct   Log10         -Inf
## 472   Male     85+           Fish   Direct   Log10         -Inf
## 473 Female   0 - 4           MeHg   Direct   Log10         -Inf
## 474   Male   0 - 4           MeHg   Direct   Log10         -Inf
## 475 Female   5 - 9           MeHg   Direct   Log10         -Inf
## 476   Male   5 - 9           MeHg   Direct   Log10         -Inf
## 477 Female 10 - 14           MeHg   Direct   Log10         -Inf
## 478   Male 10 - 14           MeHg   Direct   Log10         -Inf
## 479 Female 15 - 19           MeHg   Direct   Log10         -Inf
## 480   Male 15 - 19           MeHg   Direct   Log10         -Inf
## 481 Female 20 - 24           MeHg   Direct   Log10         -Inf
## 482   Male 20 - 24           MeHg   Direct   Log10         -Inf
## 483 Female 25 - 29           MeHg   Direct   Log10         -Inf
## 484   Male 25 - 29           MeHg   Direct   Log10         -Inf
## 485 Female 30 - 34           MeHg   Direct   Log10         -Inf
## 486   Male 30 - 34           MeHg   Direct   Log10         -Inf
## 487 Female 35 - 39           MeHg   Direct   Log10         -Inf
## 488   Male 35 - 39           MeHg   Direct   Log10         -Inf
## 489 Female 40 - 44           MeHg   Direct   Log10         -Inf
## 490   Male 40 - 44           MeHg   Direct   Log10         -Inf
## 491 Female 45 - 49           MeHg   Direct   Log10         -Inf
## 492   Male 45 - 49           MeHg   Direct   Log10         -Inf
## 493 Female 50 - 54           MeHg   Direct   Log10         -Inf
## 494   Male 50 - 54           MeHg   Direct   Log10         -Inf
## 495 Female 55 - 59           MeHg   Direct   Log10         -Inf
## 496   Male 55 - 59           MeHg   Direct   Log10         -Inf
## 497 Female 60 - 64           MeHg   Direct   Log10         -Inf
## 498   Male 60 - 64           MeHg   Direct   Log10         -Inf
## 499 Female 65 - 69           MeHg   Direct   Log10         -Inf
## 500   Male 65 - 69           MeHg   Direct   Log10         -Inf
## 501 Female 70 - 74           MeHg   Direct   Log10         -Inf
## 502   Male 70 - 74           MeHg   Direct   Log10         -Inf
## 503 Female 75 - 79           MeHg   Direct   Log10         -Inf
## 504   Male 75 - 79           MeHg   Direct   Log10         -Inf
## 505 Female 80 - 84           MeHg   Direct   Log10         -Inf
## 506   Male 80 - 84           MeHg   Direct   Log10         -Inf
## 507 Female     85+           MeHg   Direct   Log10         -Inf
## 508   Male     85+           MeHg   Direct   Log10         -Inf
## 509 Female   0 - 4         Omega3   Direct   Log10 1.941662e+00
## 510   Male   0 - 4         Omega3   Direct   Log10 2.070025e+00
## 511 Female   5 - 9         Omega3   Direct   Log10 2.123304e+00
## 512   Male   5 - 9         Omega3   Direct   Log10 2.254701e+00
## 513 Female 10 - 14         Omega3   Direct   Log10 2.196235e+00
## 514   Male 10 - 14         Omega3   Direct   Log10 2.311898e+00
## 515 Female 15 - 19         Omega3   Direct   Log10 2.196235e+00
## 516   Male 15 - 19         Omega3   Direct   Log10 2.311898e+00
## 517 Female 20 - 24         Omega3   Direct   Log10 2.166542e+00
## 518   Male 20 - 24         Omega3   Direct   Log10 2.300500e+00
## 519 Female 25 - 29         Omega3   Direct   Log10 2.166542e+00
## 520   Male 25 - 29         Omega3   Direct   Log10 2.300500e+00
## 521 Female 30 - 34         Omega3   Direct   Log10 2.166542e+00
## 522   Male 30 - 34         Omega3   Direct   Log10 2.300500e+00
## 523 Female 35 - 39         Omega3   Direct   Log10 2.166542e+00
## 524   Male 35 - 39         Omega3   Direct   Log10 2.300500e+00
## 525 Female 40 - 44         Omega3   Direct   Log10 2.166542e+00
## 526   Male 40 - 44         Omega3   Direct   Log10 2.300500e+00
## 527 Female 45 - 49         Omega3   Direct   Log10 2.166542e+00
## 528   Male 45 - 49         Omega3   Direct   Log10 2.300500e+00
## 529 Female 50 - 54         Omega3   Direct   Log10 2.166542e+00
## 530   Male 50 - 54         Omega3   Direct   Log10 2.300500e+00
## 531 Female 55 - 59         Omega3   Direct   Log10 2.166542e+00
## 532   Male 55 - 59         Omega3   Direct   Log10 2.300500e+00
## 533 Female 60 - 64         Omega3   Direct   Log10 2.166542e+00
## 534   Male 60 - 64         Omega3   Direct   Log10 2.300500e+00
## 535 Female 65 - 69         Omega3   Direct   Log10 2.250029e+00
## 536   Male 65 - 69         Omega3   Direct   Log10 2.360819e+00
## 537 Female 70 - 74         Omega3   Direct   Log10 2.250029e+00
## 538   Male 70 - 74         Omega3   Direct   Log10 2.360819e+00
## 539 Female 75 - 79         Omega3   Direct   Log10 2.250029e+00
## 540   Male 75 - 79         Omega3   Direct   Log10 2.360819e+00
## 541 Female 80 - 84         Omega3   Direct   Log10 2.250029e+00
## 542   Male 80 - 84         Omega3   Direct   Log10 2.360819e+00
## 543 Female     85+         Omega3   Direct   Log10 2.250029e+00
## 544   Male     85+         Omega3   Direct   Log10 2.360819e+00
## 545 Female   0 - 4           PFAS   Direct   Log10         -Inf
## 546   Male   0 - 4           PFAS   Direct   Log10         -Inf
## 547 Female   5 - 9           PFAS   Direct   Log10         -Inf
## 548   Male   5 - 9           PFAS   Direct   Log10         -Inf
## 549 Female 10 - 14           PFAS   Direct   Log10         -Inf
## 550   Male 10 - 14           PFAS   Direct   Log10         -Inf
## 551 Female 15 - 19           PFAS   Direct   Log10         -Inf
## 552   Male 15 - 19           PFAS   Direct   Log10         -Inf
## 553 Female 20 - 24           PFAS   Direct   Log10         -Inf
## 554   Male 20 - 24           PFAS   Direct   Log10         -Inf
## 555 Female 25 - 29           PFAS   Direct   Log10         -Inf
## 556   Male 25 - 29           PFAS   Direct   Log10         -Inf
## 557 Female 30 - 34           PFAS   Direct   Log10         -Inf
## 558   Male 30 - 34           PFAS   Direct   Log10         -Inf
## 559 Female 35 - 39           PFAS   Direct   Log10         -Inf
## 560   Male 35 - 39           PFAS   Direct   Log10         -Inf
## 561 Female 40 - 44           PFAS   Direct   Log10         -Inf
## 562   Male 40 - 44           PFAS   Direct   Log10         -Inf
## 563 Female 45 - 49           PFAS   Direct   Log10         -Inf
## 564   Male 45 - 49           PFAS   Direct   Log10         -Inf
## 565 Female 50 - 54           PFAS   Direct   Log10         -Inf
## 566   Male 50 - 54           PFAS   Direct   Log10         -Inf
## 567 Female 55 - 59           PFAS   Direct   Log10         -Inf
## 568   Male 55 - 59           PFAS   Direct   Log10         -Inf
## 569 Female 60 - 64           PFAS   Direct   Log10         -Inf
## 570   Male 60 - 64           PFAS   Direct   Log10         -Inf
## 571 Female 65 - 69           PFAS   Direct   Log10         -Inf
## 572   Male 65 - 69           PFAS   Direct   Log10         -Inf
## 573 Female 70 - 74           PFAS   Direct   Log10         -Inf
## 574   Male 70 - 74           PFAS   Direct   Log10         -Inf
## 575 Female 75 - 79           PFAS   Direct   Log10         -Inf
## 576   Male 75 - 79           PFAS   Direct   Log10         -Inf
## 577 Female 80 - 84           PFAS   Direct   Log10         -Inf
## 578   Male 80 - 84           PFAS   Direct   Log10         -Inf
## 579 Female     85+           PFAS   Direct   Log10         -Inf
## 580   Male     85+           PFAS   Direct   Log10         -Inf
## 581 Female   0 - 4            TEQ   Direct   Log10         -Inf
## 582   Male   0 - 4            TEQ   Direct   Log10         -Inf
## 583 Female   5 - 9            TEQ   Direct   Log10         -Inf
## 584   Male   5 - 9            TEQ   Direct   Log10         -Inf
## 585 Female 10 - 14            TEQ   Direct   Log10         -Inf
## 586   Male 10 - 14            TEQ   Direct   Log10         -Inf
## 587 Female 15 - 19            TEQ   Direct   Log10         -Inf
## 588   Male 15 - 19            TEQ   Direct   Log10         -Inf
## 589 Female 20 - 24            TEQ   Direct   Log10         -Inf
## 590   Male 20 - 24            TEQ   Direct   Log10         -Inf
## 591 Female 25 - 29            TEQ   Direct   Log10         -Inf
## 592   Male 25 - 29            TEQ   Direct   Log10         -Inf
## 593 Female 30 - 34            TEQ   Direct   Log10         -Inf
## 594   Male 30 - 34            TEQ   Direct   Log10         -Inf
## 595 Female 35 - 39            TEQ   Direct   Log10         -Inf
## 596   Male 35 - 39            TEQ   Direct   Log10         -Inf
## 597 Female 40 - 44            TEQ   Direct   Log10         -Inf
## 598   Male 40 - 44            TEQ   Direct   Log10         -Inf
## 599 Female 45 - 49            TEQ   Direct   Log10         -Inf
## 600   Male 45 - 49            TEQ   Direct   Log10         -Inf
## 601 Female 50 - 54            TEQ   Direct   Log10         -Inf
## 602   Male 50 - 54            TEQ   Direct   Log10         -Inf
## 603 Female 55 - 59            TEQ   Direct   Log10         -Inf
## 604   Male 55 - 59            TEQ   Direct   Log10         -Inf
## 605 Female 60 - 64            TEQ   Direct   Log10         -Inf
## 606   Male 60 - 64            TEQ   Direct   Log10         -Inf
## 607 Female 65 - 69            TEQ   Direct   Log10         -Inf
## 608   Male 65 - 69            TEQ   Direct   Log10         -Inf
## 609 Female 70 - 74            TEQ   Direct   Log10         -Inf
## 610   Male 70 - 74            TEQ   Direct   Log10         -Inf
## 611 Female 75 - 79            TEQ   Direct   Log10         -Inf
## 612   Male 75 - 79            TEQ   Direct   Log10         -Inf
## 613 Female 80 - 84            TEQ   Direct   Log10         -Inf
## 614   Male 80 - 84            TEQ   Direct   Log10         -Inf
## 615 Female     85+            TEQ   Direct   Log10         -Inf
## 616   Male     85+            TEQ   Direct   Log10         -Inf
## 617 Female   0 - 4      Vitamin D   Direct   Log10 4.158916e-01
## 618   Male   0 - 4      Vitamin D   Direct   Log10 6.077755e-01
## 619 Female   5 - 9      Vitamin D   Direct   Log10 5.385490e-01
## 620   Male   5 - 9      Vitamin D   Direct   Log10 7.129671e-01
## 621 Female 10 - 14      Vitamin D   Direct   Log10 5.509670e-01
## 622   Male 10 - 14      Vitamin D   Direct   Log10 7.269154e-01
## 623 Female 15 - 19      Vitamin D   Direct   Log10 5.509670e-01
## 624   Male 15 - 19      Vitamin D   Direct   Log10 7.269154e-01
## 625 Female 20 - 24      Vitamin D   Direct   Log10 5.526277e-01
## 626   Male 20 - 24      Vitamin D   Direct   Log10 7.410014e-01
## 627 Female 25 - 29      Vitamin D   Direct   Log10 5.526277e-01
## 628   Male 25 - 29      Vitamin D   Direct   Log10 7.410014e-01
## 629 Female 30 - 34      Vitamin D   Direct   Log10 5.526277e-01
## 630   Male 30 - 34      Vitamin D   Direct   Log10 7.410014e-01
## 631 Female 35 - 39      Vitamin D   Direct   Log10 5.526277e-01
## 632   Male 35 - 39      Vitamin D   Direct   Log10 7.410014e-01
## 633 Female 40 - 44      Vitamin D   Direct   Log10 5.526277e-01
## 634   Male 40 - 44      Vitamin D   Direct   Log10 7.410014e-01
## 635 Female 45 - 49      Vitamin D   Direct   Log10 5.526277e-01
## 636   Male 45 - 49      Vitamin D   Direct   Log10 7.410014e-01
## 637 Female 50 - 54      Vitamin D   Direct   Log10 5.526277e-01
## 638   Male 50 - 54      Vitamin D   Direct   Log10 7.410014e-01
## 639 Female 55 - 59      Vitamin D   Direct   Log10 5.526277e-01
## 640   Male 55 - 59      Vitamin D   Direct   Log10 7.410014e-01
## 641 Female 60 - 64      Vitamin D   Direct   Log10 5.526277e-01
## 642   Male 60 - 64      Vitamin D   Direct   Log10 7.410014e-01
## 643 Female 65 - 69      Vitamin D   Direct   Log10 5.958847e-01
## 644   Male 65 - 69      Vitamin D   Direct   Log10 7.826329e-01
## 645 Female 70 - 74      Vitamin D   Direct   Log10 5.958847e-01
## 646   Male 70 - 74      Vitamin D   Direct   Log10 7.826329e-01
## 647 Female 75 - 79      Vitamin D   Direct   Log10 5.958847e-01
## 648   Male 75 - 79      Vitamin D   Direct   Log10 7.826329e-01
## 649 Female 80 - 84      Vitamin D   Direct   Log10 5.958847e-01
## 650   Male 80 - 84      Vitamin D   Direct   Log10 7.826329e-01
## 651 Female     85+      Vitamin D   Direct   Log10 5.958847e-01
## 652   Male     85+      Vitamin D   Direct   Log10 7.826329e-01
## 653 Female   0 - 4           PFAS To child   Log10         -Inf
## 654   Male   0 - 4           PFAS To child   Log10         -Inf
## 655 Female   0 - 4            TEQ To child   Log10         -Inf
## 656   Male   0 - 4            TEQ To child   Log10         -Inf
## 657 Female   0 - 4            ALA   Direct    None 2.143371e+01
## 658   Male   0 - 4            ALA   Direct    None 2.899855e+01
## 659 Female   5 - 9            ALA   Direct    None 3.968659e+01
## 660   Male   5 - 9            ALA   Direct    None 5.369362e+01
## 661 Female 10 - 14            ALA   Direct    None 6.038596e+01
## 662   Male 10 - 14            ALA   Direct    None 8.169866e+01
## 663 Female 15 - 19            ALA   Direct    None 6.038596e+01
## 664   Male 15 - 19            ALA   Direct    None 8.169866e+01
## 665 Female 20 - 24            ALA   Direct    None 9.975251e+01
## 666   Male 20 - 24            ALA   Direct    None 1.349593e+02
## 667 Female 25 - 29            ALA   Direct    None 9.975251e+01
## 668   Male 25 - 29            ALA   Direct    None 1.349593e+02
## 669 Female 30 - 34            ALA   Direct    None 9.975251e+01
## 670   Male 30 - 34            ALA   Direct    None 1.349593e+02
## 671 Female 35 - 39            ALA   Direct    None 9.975251e+01
## 672   Male 35 - 39            ALA   Direct    None 1.349593e+02
## 673 Female 40 - 44            ALA   Direct    None 9.975251e+01
## 674   Male 40 - 44            ALA   Direct    None 1.349593e+02
## 675 Female 45 - 49            ALA   Direct    None 9.975251e+01
## 676   Male 45 - 49            ALA   Direct    None 1.349593e+02
## 677 Female 50 - 54            ALA   Direct    None 9.975251e+01
## 678   Male 50 - 54            ALA   Direct    None 1.349593e+02
## 679 Female 55 - 59            ALA   Direct    None 9.975251e+01
## 680   Male 55 - 59            ALA   Direct    None 1.349593e+02
## 681 Female 60 - 64            ALA   Direct    None 9.975251e+01
## 682   Male 60 - 64            ALA   Direct    None 1.349593e+02
## 683 Female 65 - 69            ALA   Direct    None 1.385860e+02
## 684   Male 65 - 69            ALA   Direct    None 1.874987e+02
## 685 Female 70 - 74            ALA   Direct    None 1.385860e+02
## 686   Male 70 - 74            ALA   Direct    None 1.874987e+02
## 687 Female 75 - 79            ALA   Direct    None 1.385860e+02
## 688   Male 75 - 79            ALA   Direct    None 1.874987e+02
## 689 Female 80 - 84            ALA   Direct    None 1.385860e+02
## 690   Male 80 - 84            ALA   Direct    None 1.874987e+02
## 691 Female     85+            ALA   Direct    None 1.385860e+02
## 692   Male     85+            ALA   Direct    None 1.874987e+02
## 693 Female   0 - 4            DHA   Direct    None 5.553134e+01
## 694   Male   0 - 4            DHA   Direct    None 7.652973e+01
## 695 Female   5 - 9            DHA   Direct    None 9.181058e+01
## 696   Male   5 - 9            DHA   Direct    None 1.256134e+02
## 697 Female 10 - 14            DHA   Direct    None 1.039071e+02
## 698   Male 10 - 14            DHA   Direct    None 1.419793e+02
## 699 Female 15 - 19            DHA   Direct    None 1.039071e+02
## 700   Male 15 - 19            DHA   Direct    None 1.419793e+02
## 701 Female 20 - 24            DHA   Direct    None 1.509747e+02
## 702   Male 20 - 24            DHA   Direct    None 2.056590e+02
## 703 Female 25 - 29            DHA   Direct    None 1.509747e+02
## 704   Male 25 - 29            DHA   Direct    None 2.056590e+02
## 705 Female 30 - 34            DHA   Direct    None 1.509747e+02
## 706   Male 30 - 34            DHA   Direct    None 2.056590e+02
## 707 Female 35 - 39            DHA   Direct    None 1.509747e+02
## 708   Male 35 - 39            DHA   Direct    None 2.056590e+02
## 709 Female 40 - 44            DHA   Direct    None 1.509747e+02
## 710   Male 40 - 44            DHA   Direct    None 2.056590e+02
## 711 Female 45 - 49            DHA   Direct    None 1.509747e+02
## 712   Male 45 - 49            DHA   Direct    None 2.056590e+02
## 713 Female 50 - 54            DHA   Direct    None 1.509747e+02
## 714   Male 50 - 54            DHA   Direct    None 2.056590e+02
## 715 Female 55 - 59            DHA   Direct    None 1.509747e+02
## 716   Male 55 - 59            DHA   Direct    None 2.056590e+02
## 717 Female 60 - 64            DHA   Direct    None 1.509747e+02
## 718   Male 60 - 64            DHA   Direct    None 2.056590e+02
## 719 Female 65 - 69            DHA   Direct    None 1.948921e+02
## 720   Male 65 - 69            DHA   Direct    None 2.650766e+02
## 721 Female 70 - 74            DHA   Direct    None 1.948921e+02
## 722   Male 70 - 74            DHA   Direct    None 2.650766e+02
## 723 Female 75 - 79            DHA   Direct    None 1.948921e+02
## 724   Male 75 - 79            DHA   Direct    None 2.650766e+02
## 725 Female 80 - 84            DHA   Direct    None 1.948921e+02
## 726   Male 80 - 84            DHA   Direct    None 2.650766e+02
## 727 Female     85+            DHA   Direct    None 1.948921e+02
## 728   Male     85+            DHA   Direct    None 2.650766e+02
## 729 Female   0 - 4            EPA   Direct    None 2.971007e+01
## 730   Male   0 - 4            EPA   Direct    None 4.347215e+01
## 731 Female   5 - 9            EPA   Direct    None 2.971007e+01
## 732   Male   5 - 9            EPA   Direct    None 4.347215e+01
## 733 Female 10 - 14            EPA   Direct    None 2.971007e+01
## 734   Male 10 - 14            EPA   Direct    None 4.347215e+01
## 735 Female 15 - 19            EPA   Direct    None 2.971007e+01
## 736   Male 15 - 19            EPA   Direct    None 4.347215e+01
## 737 Female 20 - 24            EPA   Direct    None 2.971007e+01
## 738   Male 20 - 24            EPA   Direct    None 4.347215e+01
## 739 Female 25 - 29            EPA   Direct    None 2.971007e+01
## 740   Male 25 - 29            EPA   Direct    None 4.347215e+01
## 741 Female 30 - 34            EPA   Direct    None 2.971007e+01
## 742   Male 30 - 34            EPA   Direct    None 4.347215e+01
## 743 Female 35 - 39            EPA   Direct    None 2.971007e+01
## 744   Male 35 - 39            EPA   Direct    None 4.347215e+01
## 745 Female 40 - 44            EPA   Direct    None 2.971007e+01
## 746   Male 40 - 44            EPA   Direct    None 4.347215e+01
## 747 Female 45 - 49            EPA   Direct    None 2.971007e+01
## 748   Male 45 - 49            EPA   Direct    None 4.347215e+01
## 749 Female 50 - 54            EPA   Direct    None 2.971007e+01
## 750   Male 50 - 54            EPA   Direct    None 4.347215e+01
## 751 Female 55 - 59            EPA   Direct    None 2.971007e+01
## 752   Male 55 - 59            EPA   Direct    None 4.347215e+01
## 753 Female 60 - 64            EPA   Direct    None 2.971007e+01
## 754   Male 60 - 64            EPA   Direct    None 4.347215e+01
## 755 Female 65 - 69            EPA   Direct    None 2.971007e+01
## 756   Male 65 - 69            EPA   Direct    None 4.347215e+01
## 757 Female 70 - 74            EPA   Direct    None 2.971007e+01
## 758   Male 70 - 74            EPA   Direct    None 4.347215e+01
## 759 Female 75 - 79            EPA   Direct    None 2.971007e+01
## 760   Male 75 - 79            EPA   Direct    None 4.347215e+01
## 761 Female 80 - 84            EPA   Direct    None 2.971007e+01
## 762   Male 80 - 84            EPA   Direct    None 4.347215e+01
## 763 Female     85+            EPA   Direct    None 2.971007e+01
## 764   Male     85+            EPA   Direct    None 4.347215e+01
## 765 Female   0 - 4           Fish   Direct    None 5.762180e+00
## 766   Male   0 - 4           Fish   Direct    None 7.795891e+00
## 767 Female   5 - 9           Fish   Direct    None 1.756988e+01
## 768   Male   5 - 9           Fish   Direct    None 2.377102e+01
## 769 Female 10 - 14           Fish   Direct    None 1.761415e+01
## 770   Male 10 - 14           Fish   Direct    None 2.383090e+01
## 771 Female 15 - 19           Fish   Direct    None 1.761415e+01
## 772   Male 15 - 19           Fish   Direct    None 2.383090e+01
## 773 Female 20 - 24           Fish   Direct    None 2.767019e+01
## 774   Male 20 - 24           Fish   Direct    None 3.743614e+01
## 775 Female 25 - 29           Fish   Direct    None 2.767019e+01
## 776   Male 25 - 29           Fish   Direct    None 3.743614e+01
## 777 Female 30 - 34           Fish   Direct    None 2.767019e+01
## 778   Male 30 - 34           Fish   Direct    None 3.743614e+01
## 779 Female 35 - 39           Fish   Direct    None 2.767019e+01
## 780   Male 35 - 39           Fish   Direct    None 3.743614e+01
## 781 Female 40 - 44           Fish   Direct    None 2.767019e+01
## 782   Male 40 - 44           Fish   Direct    None 3.743614e+01
## 783 Female 45 - 49           Fish   Direct    None 2.767019e+01
## 784   Male 45 - 49           Fish   Direct    None 3.743614e+01
## 785 Female 50 - 54           Fish   Direct    None 2.767019e+01
## 786   Male 50 - 54           Fish   Direct    None 3.743614e+01
## 787 Female 55 - 59           Fish   Direct    None 2.767019e+01
## 788   Male 55 - 59           Fish   Direct    None 3.743614e+01
## 789 Female 60 - 64           Fish   Direct    None 2.767019e+01
## 790   Male 60 - 64           Fish   Direct    None 3.743614e+01
## 791 Female 65 - 69           Fish   Direct    None 3.577171e+01
## 792   Male 65 - 69           Fish   Direct    None 4.839701e+01
## 793 Female 70 - 74           Fish   Direct    None 3.577171e+01
## 794   Male 70 - 74           Fish   Direct    None 4.839701e+01
## 795 Female 75 - 79           Fish   Direct    None 3.577171e+01
## 796   Male 75 - 79           Fish   Direct    None 4.839701e+01
## 797 Female 80 - 84           Fish   Direct    None 3.577171e+01
## 798   Male 80 - 84           Fish   Direct    None 4.839701e+01
## 799 Female     85+           Fish   Direct    None 3.577171e+01
## 800   Male     85+           Fish   Direct    None 4.839701e+01
## 801 Female   0 - 4           MeHg   Direct    None 5.971765e-01
## 802   Male   0 - 4           MeHg   Direct    None 8.079447e-01
## 803 Female   5 - 9           MeHg   Direct    None 1.327085e+00
## 804   Male   5 - 9           MeHg   Direct    None 1.795468e+00
## 805 Female 10 - 14           MeHg   Direct    None 1.631512e+00
## 806   Male 10 - 14           MeHg   Direct    None 2.207340e+00
## 807 Female 15 - 19           MeHg   Direct    None 1.631512e+00
## 808   Male 15 - 19           MeHg   Direct    None 2.207340e+00
## 809 Female 20 - 24           MeHg   Direct    None 2.828414e+00
## 810   Male 20 - 24           MeHg   Direct    None 3.826678e+00
## 811 Female 25 - 29           MeHg   Direct    None 2.828414e+00
## 812   Male 25 - 29           MeHg   Direct    None 3.826678e+00
## 813 Female 30 - 34           MeHg   Direct    None 2.828414e+00
## 814   Male 30 - 34           MeHg   Direct    None 3.826678e+00
## 815 Female 35 - 39           MeHg   Direct    None 2.828414e+00
## 816   Male 35 - 39           MeHg   Direct    None 3.826678e+00
## 817 Female 40 - 44           MeHg   Direct    None 2.828414e+00
## 818   Male 40 - 44           MeHg   Direct    None 3.826678e+00
## 819 Female 45 - 49           MeHg   Direct    None 2.828414e+00
## 820   Male 45 - 49           MeHg   Direct    None 3.826678e+00
## 821 Female 50 - 54           MeHg   Direct    None 2.828414e+00
## 822   Male 50 - 54           MeHg   Direct    None 3.826678e+00
## 823 Female 55 - 59           MeHg   Direct    None 2.828414e+00
## 824   Male 55 - 59           MeHg   Direct    None 3.826678e+00
## 825 Female 60 - 64           MeHg   Direct    None 2.828414e+00
## 826   Male 60 - 64           MeHg   Direct    None 3.826678e+00
## 827 Female 65 - 69           MeHg   Direct    None 3.742545e+00
## 828   Male 65 - 69           MeHg   Direct    None 5.063444e+00
## 829 Female 70 - 74           MeHg   Direct    None 3.742545e+00
## 830   Male 70 - 74           MeHg   Direct    None 5.063444e+00
## 831 Female 75 - 79           MeHg   Direct    None 3.742545e+00
## 832   Male 75 - 79           MeHg   Direct    None 5.063444e+00
## 833 Female 80 - 84           MeHg   Direct    None 3.742545e+00
## 834   Male 80 - 84           MeHg   Direct    None 5.063444e+00
## 835 Female     85+           MeHg   Direct    None 3.742545e+00
## 836   Male     85+           MeHg   Direct    None 5.063444e+00
## 837 Female   0 - 4         Omega3   Direct    None 1.414675e+02
## 838   Male   0 - 4         Omega3   Direct    None 1.927262e+02
## 839 Female   5 - 9         Omega3   Direct    None 2.488197e+02
## 840   Male   5 - 9         Omega3   Direct    None 3.379674e+02
## 841 Female 10 - 14         Omega3   Direct    None 2.956121e+02
## 842   Male 10 - 14         Omega3   Direct    None 4.012748e+02
## 843 Female 15 - 19         Omega3   Direct    None 2.956121e+02
## 844   Male 15 - 19         Omega3   Direct    None 4.012748e+02
## 845 Female 20 - 24         Omega3   Direct    None 4.439103e+02
## 846   Male 20 - 24         Omega3   Direct    None 6.019136e+02
## 847 Female 25 - 29         Omega3   Direct    None 4.439103e+02
## 848   Male 25 - 29         Omega3   Direct    None 6.019136e+02
## 849 Female 30 - 34         Omega3   Direct    None 4.439103e+02
## 850   Male 30 - 34         Omega3   Direct    None 6.019136e+02
## 851 Female 35 - 39         Omega3   Direct    None 4.439103e+02
## 852   Male 35 - 39         Omega3   Direct    None 6.019136e+02
## 853 Female 40 - 44         Omega3   Direct    None 4.439103e+02
## 854   Male 40 - 44         Omega3   Direct    None 6.019136e+02
## 855 Female 45 - 49         Omega3   Direct    None 4.439103e+02
## 856   Male 45 - 49         Omega3   Direct    None 6.019136e+02
## 857 Female 50 - 54         Omega3   Direct    None 4.439103e+02
## 858   Male 50 - 54         Omega3   Direct    None 6.019136e+02
## 859 Female 55 - 59         Omega3   Direct    None 4.439103e+02
## 860   Male 55 - 59         Omega3   Direct    None 6.019136e+02
## 861 Female 60 - 64         Omega3   Direct    None 4.439103e+02
## 862   Male 60 - 64         Omega3   Direct    None 6.019136e+02
## 863 Female 65 - 69         Omega3   Direct    None 5.991341e+02
## 864   Male 65 - 69         Omega3   Direct    None 8.119222e+02
## 865 Female 70 - 74         Omega3   Direct    None 5.991341e+02
## 866   Male 70 - 74         Omega3   Direct    None 8.119222e+02
## 867 Female 75 - 79         Omega3   Direct    None 5.991341e+02
## 868   Male 75 - 79         Omega3   Direct    None 8.119222e+02
## 869 Female 80 - 84         Omega3   Direct    None 5.991341e+02
## 870   Male 80 - 84         Omega3   Direct    None 8.119222e+02
## 871 Female     85+         Omega3   Direct    None 5.991341e+02
## 872   Male     85+         Omega3   Direct    None 8.119222e+02
## 873 Female   0 - 4           PFAS   Direct    None 4.724958e+01
## 874   Male   0 - 4           PFAS   Direct    None 6.392590e+01
## 875 Female   5 - 9           PFAS   Direct    None 1.494611e+02
## 876   Male   5 - 9           PFAS   Direct    None 2.022121e+02
## 877 Female 10 - 14           PFAS   Direct    None 1.517132e+02
## 878   Male 10 - 14           PFAS   Direct    None 2.052591e+02
## 879 Female 15 - 19           PFAS   Direct    None 1.517132e+02
## 880   Male 15 - 19           PFAS   Direct    None 2.052591e+02
## 881 Female 20 - 24           PFAS   Direct    None 2.307378e+02
## 882   Male 20 - 24           PFAS   Direct    None 3.121747e+02
## 883 Female 25 - 29           PFAS   Direct    None 2.307378e+02
## 884   Male 25 - 29           PFAS   Direct    None 3.121747e+02
## 885 Female 30 - 34           PFAS   Direct    None 2.307378e+02
## 886   Male 30 - 34           PFAS   Direct    None 3.121747e+02
## 887 Female 35 - 39           PFAS   Direct    None 2.307378e+02
## 888   Male 35 - 39           PFAS   Direct    None 3.121747e+02
## 889 Female 40 - 44           PFAS   Direct    None 2.307378e+02
## 890   Male 40 - 44           PFAS   Direct    None 3.121747e+02
## 891 Female 45 - 49           PFAS   Direct    None 2.307378e+02
## 892   Male 45 - 49           PFAS   Direct    None 3.121747e+02
## 893 Female 50 - 54           PFAS   Direct    None 2.307378e+02
## 894   Male 50 - 54           PFAS   Direct    None 3.121747e+02
## 895 Female 55 - 59           PFAS   Direct    None 2.307378e+02
## 896   Male 55 - 59           PFAS   Direct    None 3.121747e+02
## 897 Female 60 - 64           PFAS   Direct    None 2.307378e+02
## 898   Male 60 - 64           PFAS   Direct    None 3.121747e+02
## 899 Female 65 - 69           PFAS   Direct    None 3.087620e+02
## 900   Male 65 - 69           PFAS   Direct    None 4.177368e+02
## 901 Female 70 - 74           PFAS   Direct    None 3.087620e+02
## 902   Male 70 - 74           PFAS   Direct    None 4.177368e+02
## 903 Female 75 - 79           PFAS   Direct    None 3.087620e+02
## 904   Male 75 - 79           PFAS   Direct    None 4.177368e+02
## 905 Female 80 - 84           PFAS   Direct    None 3.087620e+02
## 906   Male 80 - 84           PFAS   Direct    None 4.177368e+02
## 907 Female     85+           PFAS   Direct    None 3.087620e+02
## 908   Male     85+           PFAS   Direct    None 4.177368e+02
## 909 Female   0 - 4            TEQ   Direct    None 5.425332e+01
## 910   Male   0 - 4            TEQ   Direct    None 7.340156e+01
## 911 Female   5 - 9            TEQ   Direct    None 9.782083e+01
## 912   Male   5 - 9            TEQ   Direct    None 1.323458e+02
## 913 Female 10 - 14            TEQ   Direct    None 1.601890e+02
## 914   Male 10 - 14            TEQ   Direct    None 2.167264e+02
## 915 Female 15 - 19            TEQ   Direct    None 1.601890e+02
## 916   Male 15 - 19            TEQ   Direct    None 2.167264e+02
## 917 Female 20 - 24            TEQ   Direct    None 3.018922e+02
## 918   Male 20 - 24            TEQ   Direct    None 4.084424e+02
## 919 Female 25 - 29            TEQ   Direct    None 3.018922e+02
## 920   Male 25 - 29            TEQ   Direct    None 4.084424e+02
## 921 Female 30 - 34            TEQ   Direct    None 3.018922e+02
## 922   Male 30 - 34            TEQ   Direct    None 4.084424e+02
## 923 Female 35 - 39            TEQ   Direct    None 3.018922e+02
## 924   Male 35 - 39            TEQ   Direct    None 4.084424e+02
## 925 Female 40 - 44            TEQ   Direct    None 3.018922e+02
## 926   Male 40 - 44            TEQ   Direct    None 4.084424e+02
## 927 Female 45 - 49            TEQ   Direct    None 3.018922e+02
## 928   Male 45 - 49            TEQ   Direct    None 4.084424e+02
## 929 Female 50 - 54            TEQ   Direct    None 3.018922e+02
## 930   Male 50 - 54            TEQ   Direct    None 4.084424e+02
## 931 Female 55 - 59            TEQ   Direct    None 3.018922e+02
## 932   Male 55 - 59            TEQ   Direct    None 4.084424e+02
## 933 Female 60 - 64            TEQ   Direct    None 3.018922e+02
## 934   Male 60 - 64            TEQ   Direct    None 4.084424e+02
## 935 Female 65 - 69            TEQ   Direct    None 5.412734e+02
## 936   Male 65 - 69            TEQ   Direct    None 7.323111e+02
## 937 Female 70 - 74            TEQ   Direct    None 5.412734e+02
## 938   Male 70 - 74            TEQ   Direct    None 7.323111e+02
## 939 Female 75 - 79            TEQ   Direct    None 5.412734e+02
## 940   Male 75 - 79            TEQ   Direct    None 7.323111e+02
## 941 Female 80 - 84            TEQ   Direct    None 5.412734e+02
## 942   Male 80 - 84            TEQ   Direct    None 7.323111e+02
## 943 Female     85+            TEQ   Direct    None 5.412734e+02
## 944   Male     85+            TEQ   Direct    None 7.323111e+02
## 945 Female   0 - 4      Vitamin D   Direct    None 3.365585e+00
## 946   Male   0 - 4      Vitamin D   Direct    None 5.193882e+00
## 947 Female   5 - 9      Vitamin D   Direct    None 4.567127e+00
## 948   Male   5 - 9      Vitamin D   Direct    None 6.819497e+00
## 949 Female 10 - 14      Vitamin D   Direct    None 4.467767e+00
## 950   Male 10 - 14      Vitamin D   Direct    None 6.685069e+00
## 951 Female 15 - 19      Vitamin D   Direct    None 4.467767e+00
## 952   Male 15 - 19      Vitamin D   Direct    None 6.685069e+00
## 953 Female 20 - 24      Vitamin D   Direct    None 5.526103e+00
## 954   Male 20 - 24      Vitamin D   Direct    None 8.116936e+00
## 955 Female 25 - 29      Vitamin D   Direct    None 5.526103e+00
## 956   Male 25 - 29      Vitamin D   Direct    None 8.116936e+00
## 957 Female 30 - 34      Vitamin D   Direct    None 5.526103e+00
## 958   Male 30 - 34      Vitamin D   Direct    None 8.116936e+00
## 959 Female 35 - 39      Vitamin D   Direct    None 5.526103e+00
## 960   Male 35 - 39      Vitamin D   Direct    None 8.116936e+00
## 961 Female 40 - 44      Vitamin D   Direct    None 5.526103e+00
## 962   Male 40 - 44      Vitamin D   Direct    None 8.116936e+00
## 963 Female 45 - 49      Vitamin D   Direct    None 5.526103e+00
## 964   Male 45 - 49      Vitamin D   Direct    None 8.116936e+00
## 965 Female 50 - 54      Vitamin D   Direct    None 5.526103e+00
## 966   Male 50 - 54      Vitamin D   Direct    None 8.116936e+00
## 967 Female 55 - 59      Vitamin D   Direct    None 5.526103e+00
## 968   Male 55 - 59      Vitamin D   Direct    None 8.116936e+00
## 969 Female 60 - 64      Vitamin D   Direct    None 5.526103e+00
## 970   Male 60 - 64      Vitamin D   Direct    None 8.116936e+00
## 971 Female 65 - 69      Vitamin D   Direct    None 6.273107e+00
## 972   Male 65 - 69      Vitamin D   Direct    None 9.127589e+00
## 973 Female 70 - 74      Vitamin D   Direct    None 6.273107e+00
## 974   Male 70 - 74      Vitamin D   Direct    None 9.127589e+00
## 975 Female 75 - 79      Vitamin D   Direct    None 6.273107e+00
## 976   Male 75 - 79      Vitamin D   Direct    None 9.127589e+00
## 977 Female 80 - 84      Vitamin D   Direct    None 6.273107e+00
## 978   Male 80 - 84      Vitamin D   Direct    None 9.127589e+00
## 979 Female     85+      Vitamin D   Direct    None 6.273107e+00
## 980   Male     85+      Vitamin D   Direct    None 9.127589e+00
## 981 Female   0 - 4           PFAS To child    None 1.346139e+01
## 982   Male   0 - 4           PFAS To child    None 1.346139e+01
## 983 Female   0 - 4            TEQ To child    None 4.039597e+01
## 984   Male   0 - 4            TEQ To child    None 4.039597e+01
oprint(summary(ERF,"mean"))
##    Exposure_agent                                          Response
## 1            Fish                               All-cause mortality
## 2          Omega3                                     Breast cancer
## 3            Fish                                        Depression
## 4             DHA                         Loss in child's IQ points
## 5             TEQ                               Sperm concentration
## 6             TEQ                           Yes or no dental defect
## 7          Omega3                                    CHD2 mortality
## 8       Vitamin D                          Vitamin D recommendation
## 9            MeHg                         Loss in child's IQ points
## 10           PFAS                                 Immunosuppression
## 11           PFAS                                          PFAS TWI
## 12            TEQ                           Cancer morbidity yearly
## 13            TEQ Dioxin recommendation tolerable daily intake 2018
## 14           Fish                               All-cause mortality
## 15         Omega3                                     Breast cancer
## 16           Fish                                        Depression
## 17            DHA                         Loss in child's IQ points
## 18            TEQ                               Sperm concentration
## 19            TEQ                           Yes or no dental defect
## 20         Omega3                                    CHD2 mortality
## 21      Vitamin D                          Vitamin D recommendation
## 22           MeHg                         Loss in child's IQ points
## 23           PFAS                                 Immunosuppression
## 24           PFAS                                          PFAS TWI
## 25            TEQ                           Cancer morbidity yearly
## 26            TEQ Dioxin recommendation tolerable daily intake 2018
##      ER_function Scaling Observation          mean
## 1             RR    None         ERF  9.978595e-01
## 2             RR    None         ERF  9.994762e-01
## 3             RR    None         ERF  9.946188e-01
## 4            ERS    None         ERF -1.306489e-03
## 5            ERS    None         ERF  6.000075e-05
## 6            ERS    None         ERF  1.413162e-03
## 7  Relative Hill    None         ERF -1.700291e-01
## 8           Step    None         ERF  1.000000e+02
## 9            ERS      BW         ERF  9.659261e+00
## 10           ERS      BW         ERF  2.270000e-02
## 11           TWI      BW         ERF  4.400000e+00
## 12           CSF      BW         ERF  9.967577e-06
## 13           TDI      BW         ERF  2.889000e-01
## 14            RR    None   Threshold  0.000000e+00
## 15            RR    None   Threshold  0.000000e+00
## 16            RR    None   Threshold  0.000000e+00
## 17           ERS    None   Threshold  0.000000e+00
## 18           ERS    None   Threshold  0.000000e+00
## 19           ERS    None   Threshold  0.000000e+00
## 20 Relative Hill    None   Threshold  4.700000e+01
## 21          Step    None   Threshold  1.000000e+01
## 22           ERS      BW   Threshold  0.000000e+00
## 23           ERS      BW   Threshold  0.000000e+00
## 24           TWI      BW   Threshold  0.000000e+00
## 25           CSF      BW   Threshold  0.000000e+00
## 26           TDI      BW   Threshold  0.000000e+00
oprint(summary(expo_dir,"mean"))
##     Gender     Age Exposure_agent        mean
## 1   Female   0 - 4            ALA  21.4337074
## 2     Male   0 - 4            ALA  28.9985453
## 3   Female   5 - 9            ALA  39.6865892
## 4     Male   5 - 9            ALA  53.6936207
## 5   Female 10 - 14            ALA  60.3859645
## 6     Male 10 - 14            ALA  81.6986578
## 7   Female 15 - 19            ALA  60.3859645
## 8     Male 15 - 19            ALA  81.6986578
## 9   Female 20 - 24            ALA  99.7525088
## 10    Male 20 - 24            ALA 134.9592767
## 11  Female 25 - 29            ALA  99.7525088
## 12    Male 25 - 29            ALA 134.9592767
## 13  Female 30 - 34            ALA  99.7525088
## 14    Male 30 - 34            ALA 134.9592767
## 15  Female 35 - 39            ALA  99.7525088
## 16    Male 35 - 39            ALA 134.9592767
## 17  Female 40 - 44            ALA  99.7525088
## 18    Male 40 - 44            ALA 134.9592767
## 19  Female 45 - 49            ALA  99.7525088
## 20    Male 45 - 49            ALA 134.9592767
## 21  Female 50 - 54            ALA  99.7525088
## 22    Male 50 - 54            ALA 134.9592767
## 23  Female 55 - 59            ALA  99.7525088
## 24    Male 55 - 59            ALA 134.9592767
## 25  Female 60 - 64            ALA  99.7525088
## 26    Male 60 - 64            ALA 134.9592767
## 27  Female 65 - 69            ALA 138.5860138
## 28    Male 65 - 69            ALA 187.4987246
## 29  Female 70 - 74            ALA 138.5860138
## 30    Male 70 - 74            ALA 187.4987246
## 31  Female 75 - 79            ALA 138.5860138
## 32    Male 75 - 79            ALA 187.4987246
## 33  Female 80 - 84            ALA 138.5860138
## 34    Male 80 - 84            ALA 187.4987246
## 35  Female     85+            ALA 138.5860138
## 36    Male     85+            ALA 187.4987246
## 37  Female   0 - 4            DHA  55.5313401
## 38    Male   0 - 4            DHA  76.5297273
## 39  Female   5 - 9            DHA  91.8105802
## 40    Male   5 - 9            DHA 125.6134050
## 41  Female 10 - 14            DHA 103.9070801
## 42    Male 10 - 14            DHA 141.9792579
## 43  Female 15 - 19            DHA 103.9070801
## 44    Male 15 - 19            DHA 141.9792579
## 45  Female 20 - 24            DHA 150.9747434
## 46    Male 20 - 24            DHA 205.6590377
## 47  Female 25 - 29            DHA 150.9747434
## 48    Male 25 - 29            DHA 205.6590377
## 49  Female 30 - 34            DHA 150.9747434
## 50    Male 30 - 34            DHA 205.6590377
## 51  Female 35 - 39            DHA 150.9747434
## 52    Male 35 - 39            DHA 205.6590377
## 53  Female 40 - 44            DHA 150.9747434
## 54    Male 40 - 44            DHA 205.6590377
## 55  Female 45 - 49            DHA 150.9747434
## 56    Male 45 - 49            DHA 205.6590377
## 57  Female 50 - 54            DHA 150.9747434
## 58    Male 50 - 54            DHA 205.6590377
## 59  Female 55 - 59            DHA 150.9747434
## 60    Male 55 - 59            DHA 205.6590377
## 61  Female 60 - 64            DHA 150.9747434
## 62    Male 60 - 64            DHA 205.6590377
## 63  Female 65 - 69            DHA 194.8920588
## 64    Male 65 - 69            DHA 265.0765820
## 65  Female 70 - 74            DHA 194.8920588
## 66    Male 70 - 74            DHA 265.0765820
## 67  Female 75 - 79            DHA 194.8920588
## 68    Male 75 - 79            DHA 265.0765820
## 69  Female 80 - 84            DHA 194.8920588
## 70    Male 80 - 84            DHA 265.0765820
## 71  Female     85+            DHA 194.8920588
## 72    Male     85+            DHA 265.0765820
## 73  Female   0 - 4            EPA  29.7100701
## 74    Male   0 - 4            EPA  43.4721516
## 75  Female   5 - 9            EPA  29.7100701
## 76    Male   5 - 9            EPA  43.4721516
## 77  Female 10 - 14            EPA  29.7100701
## 78    Male 10 - 14            EPA  43.4721516
## 79  Female 15 - 19            EPA  29.7100701
## 80    Male 15 - 19            EPA  43.4721516
## 81  Female 20 - 24            EPA  29.7100701
## 82    Male 20 - 24            EPA  43.4721516
## 83  Female 25 - 29            EPA  29.7100701
## 84    Male 25 - 29            EPA  43.4721516
## 85  Female 30 - 34            EPA  29.7100701
## 86    Male 30 - 34            EPA  43.4721516
## 87  Female 35 - 39            EPA  29.7100701
## 88    Male 35 - 39            EPA  43.4721516
## 89  Female 40 - 44            EPA  29.7100701
## 90    Male 40 - 44            EPA  43.4721516
## 91  Female 45 - 49            EPA  29.7100701
## 92    Male 45 - 49            EPA  43.4721516
## 93  Female 50 - 54            EPA  29.7100701
## 94    Male 50 - 54            EPA  43.4721516
## 95  Female 55 - 59            EPA  29.7100701
## 96    Male 55 - 59            EPA  43.4721516
## 97  Female 60 - 64            EPA  29.7100701
## 98    Male 60 - 64            EPA  43.4721516
## 99  Female 65 - 69            EPA  29.7100701
## 100   Male 65 - 69            EPA  43.4721516
## 101 Female 70 - 74            EPA  29.7100701
## 102   Male 70 - 74            EPA  43.4721516
## 103 Female 75 - 79            EPA  29.7100701
## 104   Male 75 - 79            EPA  43.4721516
## 105 Female 80 - 84            EPA  29.7100701
## 106   Male 80 - 84            EPA  43.4721516
## 107 Female     85+            EPA  29.7100701
## 108   Male     85+            EPA  43.4721516
## 109 Female   0 - 4           Fish   5.7621802
## 110   Male   0 - 4           Fish   7.7958909
## 111 Female   5 - 9           Fish  17.5698850
## 112   Male   5 - 9           Fish  23.7710209
## 113 Female 10 - 14           Fish  17.6141458
## 114   Male 10 - 14           Fish  23.8309032
## 115 Female 15 - 19           Fish  17.6141458
## 116   Male 15 - 19           Fish  23.8309032
## 117 Female 20 - 24           Fish  27.6701889
## 118   Male 20 - 24           Fish  37.4361379
## 119 Female 25 - 29           Fish  27.6701889
## 120   Male 25 - 29           Fish  37.4361379
## 121 Female 30 - 34           Fish  27.6701889
## 122   Male 30 - 34           Fish  37.4361379
## 123 Female 35 - 39           Fish  27.6701889
## 124   Male 35 - 39           Fish  37.4361379
## 125 Female 40 - 44           Fish  27.6701889
## 126   Male 40 - 44           Fish  37.4361379
## 127 Female 45 - 49           Fish  27.6701889
## 128   Male 45 - 49           Fish  37.4361379
## 129 Female 50 - 54           Fish  27.6701889
## 130   Male 50 - 54           Fish  37.4361379
## 131 Female 55 - 59           Fish  27.6701889
## 132   Male 55 - 59           Fish  37.4361379
## 133 Female 60 - 64           Fish  27.6701889
## 134   Male 60 - 64           Fish  37.4361379
## 135 Female 65 - 69           Fish  35.7717063
## 136   Male 65 - 69           Fish  48.3970145
## 137 Female 70 - 74           Fish  35.7717063
## 138   Male 70 - 74           Fish  48.3970145
## 139 Female 75 - 79           Fish  35.7717063
## 140   Male 75 - 79           Fish  48.3970145
## 141 Female 80 - 84           Fish  35.7717063
## 142   Male 80 - 84           Fish  48.3970145
## 143 Female     85+           Fish  35.7717063
## 144   Male     85+           Fish  48.3970145
## 145 Female   0 - 4           MeHg   0.5971765
## 146   Male   0 - 4           MeHg   0.8079447
## 147 Female   5 - 9           MeHg   1.3270854
## 148   Male   5 - 9           MeHg   1.7954684
## 149 Female 10 - 14           MeHg   1.6315120
## 150   Male 10 - 14           MeHg   2.2073398
## 151 Female 15 - 19           MeHg   1.6315120
## 152   Male 15 - 19           MeHg   2.2073398
## 153 Female 20 - 24           MeHg   2.8284141
## 154   Male 20 - 24           MeHg   3.8266780
## 155 Female 25 - 29           MeHg   2.8284141
## 156   Male 25 - 29           MeHg   3.8266780
## 157 Female 30 - 34           MeHg   2.8284141
## 158   Male 30 - 34           MeHg   3.8266780
## 159 Female 35 - 39           MeHg   2.8284141
## 160   Male 35 - 39           MeHg   3.8266780
## 161 Female 40 - 44           MeHg   2.8284141
## 162   Male 40 - 44           MeHg   3.8266780
## 163 Female 45 - 49           MeHg   2.8284141
## 164   Male 45 - 49           MeHg   3.8266780
## 165 Female 50 - 54           MeHg   2.8284141
## 166   Male 50 - 54           MeHg   3.8266780
## 167 Female 55 - 59           MeHg   2.8284141
## 168   Male 55 - 59           MeHg   3.8266780
## 169 Female 60 - 64           MeHg   2.8284141
## 170   Male 60 - 64           MeHg   3.8266780
## 171 Female 65 - 69           MeHg   3.7425454
## 172   Male 65 - 69           MeHg   5.0634437
## 173 Female 70 - 74           MeHg   3.7425454
## 174   Male 70 - 74           MeHg   5.0634437
## 175 Female 75 - 79           MeHg   3.7425454
## 176   Male 75 - 79           MeHg   5.0634437
## 177 Female 80 - 84           MeHg   3.7425454
## 178   Male 80 - 84           MeHg   5.0634437
## 179 Female     85+           MeHg   3.7425454
## 180   Male     85+           MeHg   5.0634437
## 181 Female   0 - 4         Omega3 141.4674764
## 182   Male   0 - 4         Omega3 192.7261932
## 183 Female   5 - 9         Omega3 248.8196536
## 184   Male   5 - 9         Omega3 337.9673743
## 185 Female 10 - 14         Omega3 295.6121154
## 186   Male 10 - 14         Omega3 401.2748226
## 187 Female 15 - 19         Omega3 295.6121154
## 188   Male 15 - 19         Omega3 401.2748226
## 189 Female 20 - 24         Omega3 443.9103092
## 190   Male 20 - 24         Omega3 601.9135554
## 191 Female 25 - 29         Omega3 443.9103092
## 192   Male 25 - 29         Omega3 601.9135554
## 193 Female 30 - 34         Omega3 443.9103092
## 194   Male 30 - 34         Omega3 601.9135554
## 195 Female 35 - 39         Omega3 443.9103092
## 196   Male 35 - 39         Omega3 601.9135554
## 197 Female 40 - 44         Omega3 443.9103092
## 198   Male 40 - 44         Omega3 601.9135554
## 199 Female 45 - 49         Omega3 443.9103092
## 200   Male 45 - 49         Omega3 601.9135554
## 201 Female 50 - 54         Omega3 443.9103092
## 202   Male 50 - 54         Omega3 601.9135554
## 203 Female 55 - 59         Omega3 443.9103092
## 204   Male 55 - 59         Omega3 601.9135554
## 205 Female 60 - 64         Omega3 443.9103092
## 206   Male 60 - 64         Omega3 601.9135554
## 207 Female 65 - 69         Omega3 599.1340842
## 208   Male 65 - 69         Omega3 811.9221921
## 209 Female 70 - 74         Omega3 599.1340842
## 210   Male 70 - 74         Omega3 811.9221921
## 211 Female 75 - 79         Omega3 599.1340842
## 212   Male 75 - 79         Omega3 811.9221921
## 213 Female 80 - 84         Omega3 599.1340842
## 214   Male 80 - 84         Omega3 811.9221921
## 215 Female     85+         Omega3 599.1340842
## 216   Male     85+         Omega3 811.9221921
## 217 Female   0 - 4           PFAS  47.2495764
## 218   Male   0 - 4           PFAS  63.9258974
## 219 Female   5 - 9           PFAS 149.4611117
## 220   Male   5 - 9           PFAS 202.2120924
## 221 Female 10 - 14           PFAS 151.7132151
## 222   Male 10 - 14           PFAS 205.2590557
## 223 Female 15 - 19           PFAS 151.7132151
## 224   Male 15 - 19           PFAS 205.2590557
## 225 Female 20 - 24           PFAS 230.7378180
## 226   Male 20 - 24           PFAS 312.1746949
## 227 Female 25 - 29           PFAS 230.7378180
## 228   Male 25 - 29           PFAS 312.1746949
## 229 Female 30 - 34           PFAS 230.7378180
## 230   Male 30 - 34           PFAS 312.1746949
## 231 Female 35 - 39           PFAS 230.7378180
## 232   Male 35 - 39           PFAS 312.1746949
## 233 Female 40 - 44           PFAS 230.7378180
## 234   Male 40 - 44           PFAS 312.1746949
## 235 Female 45 - 49           PFAS 230.7378180
## 236   Male 45 - 49           PFAS 312.1746949
## 237 Female 50 - 54           PFAS 230.7378180
## 238   Male 50 - 54           PFAS 312.1746949
## 239 Female 55 - 59           PFAS 230.7378180
## 240   Male 55 - 59           PFAS 312.1746949
## 241 Female 60 - 64           PFAS 230.7378180
## 242   Male 60 - 64           PFAS 312.1746949
## 243 Female 65 - 69           PFAS 308.7619762
## 244   Male 65 - 69           PFAS 417.7367914
## 245 Female 70 - 74           PFAS 308.7619762
## 246   Male 70 - 74           PFAS 417.7367914
## 247 Female 75 - 79           PFAS 308.7619762
## 248   Male 75 - 79           PFAS 417.7367914
## 249 Female 80 - 84           PFAS 308.7619762
## 250   Male 80 - 84           PFAS 417.7367914
## 251 Female     85+           PFAS 308.7619762
## 252   Male     85+           PFAS 417.7367914
## 253 Female   0 - 4            TEQ  54.2533243
## 254   Male   0 - 4            TEQ  73.4015564
## 255 Female   5 - 9            TEQ  97.8208278
## 256   Male   5 - 9            TEQ 132.3458259
## 257 Female 10 - 14            TEQ 160.1890499
## 258   Male 10 - 14            TEQ 216.7263616
## 259 Female 15 - 19            TEQ 160.1890499
## 260   Male 15 - 19            TEQ 216.7263616
## 261 Female 20 - 24            TEQ 301.8922163
## 262   Male 20 - 24            TEQ 408.4424103
## 263 Female 25 - 29            TEQ 301.8922163
## 264   Male 25 - 29            TEQ 408.4424103
## 265 Female 30 - 34            TEQ 301.8922163
## 266   Male 30 - 34            TEQ 408.4424103
## 267 Female 35 - 39            TEQ 301.8922163
## 268   Male 35 - 39            TEQ 408.4424103
## 269 Female 40 - 44            TEQ 301.8922163
## 270   Male 40 - 44            TEQ 408.4424103
## 271 Female 45 - 49            TEQ 301.8922163
## 272   Male 45 - 49            TEQ 408.4424103
## 273 Female 50 - 54            TEQ 301.8922163
## 274   Male 50 - 54            TEQ 408.4424103
## 275 Female 55 - 59            TEQ 301.8922163
## 276   Male 55 - 59            TEQ 408.4424103
## 277 Female 60 - 64            TEQ 301.8922163
## 278   Male 60 - 64            TEQ 408.4424103
## 279 Female 65 - 69            TEQ 541.2734024
## 280   Male 65 - 69            TEQ 732.3110739
## 281 Female 70 - 74            TEQ 541.2734024
## 282   Male 70 - 74            TEQ 732.3110739
## 283 Female 75 - 79            TEQ 541.2734024
## 284   Male 75 - 79            TEQ 732.3110739
## 285 Female 80 - 84            TEQ 541.2734024
## 286   Male 80 - 84            TEQ 732.3110739
## 287 Female     85+            TEQ 541.2734024
## 288   Male     85+            TEQ 732.3110739
## 289 Female   0 - 4      Vitamin D   3.3655848
## 290   Male   0 - 4      Vitamin D   5.1938819
## 291 Female   5 - 9      Vitamin D   4.5671267
## 292   Male   5 - 9      Vitamin D   6.8194974
## 293 Female 10 - 14      Vitamin D   4.4677666
## 294   Male 10 - 14      Vitamin D   6.6850691
## 295 Female 15 - 19      Vitamin D   4.4677666
## 296   Male 15 - 19      Vitamin D   6.6850691
## 297 Female 20 - 24      Vitamin D   5.5261030
## 298   Male 20 - 24      Vitamin D   8.1169360
## 299 Female 25 - 29      Vitamin D   5.5261030
## 300   Male 25 - 29      Vitamin D   8.1169360
## 301 Female 30 - 34      Vitamin D   5.5261030
## 302   Male 30 - 34      Vitamin D   8.1169360
## 303 Female 35 - 39      Vitamin D   5.5261030
## 304   Male 35 - 39      Vitamin D   8.1169360
## 305 Female 40 - 44      Vitamin D   5.5261030
## 306   Male 40 - 44      Vitamin D   8.1169360
## 307 Female 45 - 49      Vitamin D   5.5261030
## 308   Male 45 - 49      Vitamin D   8.1169360
## 309 Female 50 - 54      Vitamin D   5.5261030
## 310   Male 50 - 54      Vitamin D   8.1169360
## 311 Female 55 - 59      Vitamin D   5.5261030
## 312   Male 55 - 59      Vitamin D   8.1169360
## 313 Female 60 - 64      Vitamin D   5.5261030
## 314   Male 60 - 64      Vitamin D   8.1169360
## 315 Female 65 - 69      Vitamin D   6.2731073
## 316   Male 65 - 69      Vitamin D   9.1275888
## 317 Female 70 - 74      Vitamin D   6.2731073
## 318   Male 70 - 74      Vitamin D   9.1275888
## 319 Female 75 - 79      Vitamin D   6.2731073
## 320   Male 75 - 79      Vitamin D   9.1275888
## 321 Female 80 - 84      Vitamin D   6.2731073
## 322   Male 80 - 84      Vitamin D   9.1275888
## 323 Female     85+      Vitamin D   6.2731073
## 324   Male     85+      Vitamin D   9.1275888
oprint(summary(expo_dir,marginals="Exposure_agent"))
##   Exposure_agent       mean          sd          min    Q0.025     median
## 1            ALA 115.855874  261.262993 0.000000e+00 0.0000000   0.000000
## 2            DHA 176.411249  258.069683 1.758942e-04 3.4334285  64.772856
## 3            EPA  36.591111   22.326449 5.640929e-02 1.8609909  35.280707
## 4           Fish  31.794156   54.761899 0.000000e+00 0.0000000   0.000000
## 5           MeHg   3.225865    6.881612 0.000000e+00 0.0000000   0.000000
## 6         Omega3 521.734743  851.257430 1.272163e-01 6.6898725 134.685017
## 7           PFAS 269.319435  643.381952 0.000000e+00 0.0000000   0.000000
## 8            TEQ 385.349622 1095.053634 0.000000e+00 0.0000000   0.000000
## 9      Vitamin D   6.723405    6.145554 1.841267e-05 0.3059591   5.139873
##       Q0.975         max
## 1  955.26903  3844.68000
## 2  972.25417  3275.96222
## 3   81.06874    84.99556
## 4  186.87500   483.00000
## 5   24.17034    91.19623
## 6 3193.24906 11161.25027
## 7 1891.95351 11270.57710
## 8 2977.29742 18737.56553
## 9   23.06482    58.56516
oprint(summary(expo_indir,"mean"))
##    Gender     Age Exposure_agent Exposure      mean
## 1  Female   0 - 4           PFAS To child  2.798839
## 2    Male   0 - 4           PFAS To child  3.786665
## 3  Female   5 - 9           PFAS To child  8.813907
## 4    Male   5 - 9           PFAS To child 11.924698
## 5  Female 10 - 14           PFAS To child  8.970467
## 6    Male 10 - 14           PFAS To child 12.136514
## 7  Female 15 - 19           PFAS To child  8.970467
## 8    Male 15 - 19           PFAS To child 12.136514
## 9  Female 20 - 24           PFAS To child 13.461386
## 10   Male 20 - 24           PFAS To child 18.212464
## 11 Female 25 - 29           PFAS To child 13.461386
## 12   Male 25 - 29           PFAS To child 18.212464
## 13 Female 30 - 34           PFAS To child 13.461386
## 14   Male 30 - 34           PFAS To child 18.212464
## 15 Female 35 - 39           PFAS To child 13.461386
## 16   Male 35 - 39           PFAS To child 18.212464
## 17 Female 40 - 44           PFAS To child 13.461386
## 18   Male 40 - 44           PFAS To child 18.212464
## 19 Female 45 - 49           PFAS To child 13.461386
## 20   Male 45 - 49           PFAS To child 18.212464
## 21 Female 50 - 54           PFAS To child 13.461386
## 22   Male 50 - 54           PFAS To child 18.212464
## 23 Female 55 - 59           PFAS To child 13.461386
## 24   Male 55 - 59           PFAS To child 18.212464
## 25 Female 60 - 64           PFAS To child 13.461386
## 26   Male 60 - 64           PFAS To child 18.212464
## 27 Female 65 - 69           PFAS To child 18.293769
## 28   Male 65 - 69           PFAS To child 24.750394
## 29 Female 70 - 74           PFAS To child 18.293769
## 30   Male 70 - 74           PFAS To child 24.750394
## 31 Female 75 - 79           PFAS To child 18.293769
## 32   Male 75 - 79           PFAS To child 24.750394
## 33 Female 80 - 84           PFAS To child 18.293769
## 34   Male 80 - 84           PFAS To child 24.750394
## 35 Female     85+           PFAS To child 18.293769
## 36   Male     85+           PFAS To child 24.750394
## 37 Female   0 - 4            TEQ To child  7.235977
## 38   Male   0 - 4            TEQ To child  9.789851
## 39 Female   5 - 9            TEQ To child 13.026674
## 40   Male   5 - 9            TEQ To child 17.624323
## 41 Female 10 - 14            TEQ To child 20.566721
## 42   Male 10 - 14            TEQ To child 27.825564
## 43 Female 15 - 19            TEQ To child 20.566721
## 44   Male 15 - 19            TEQ To child 27.825564
## 45 Female 20 - 24            TEQ To child 40.395969
## 46   Male 20 - 24            TEQ To child 54.653369
## 47 Female 25 - 29            TEQ To child 40.395969
## 48   Male 25 - 29            TEQ To child 54.653369
## 49 Female 30 - 34            TEQ To child 40.395969
## 50   Male 30 - 34            TEQ To child 54.653369
## 51 Female 35 - 39            TEQ To child 40.395969
## 52   Male 35 - 39            TEQ To child 54.653369
## 53 Female 40 - 44            TEQ To child 40.395969
## 54   Male 40 - 44            TEQ To child 54.653369
## 55 Female 45 - 49            TEQ To child 40.395969
## 56   Male 45 - 49            TEQ To child 54.653369
## 57 Female 50 - 54            TEQ To child 40.395969
## 58   Male 50 - 54            TEQ To child 54.653369
## 59 Female 55 - 59            TEQ To child 40.395969
## 60   Male 55 - 59            TEQ To child 54.653369
## 61 Female 60 - 64            TEQ To child 40.395969
## 62   Male 60 - 64            TEQ To child 54.653369
## 63 Female 65 - 69            TEQ To child 72.655053
## 64   Male 65 - 69            TEQ To child 98.298013
## 65 Female 70 - 74            TEQ To child 72.655053
## 66   Male 70 - 74            TEQ To child 98.298013
## 67 Female 75 - 79            TEQ To child 72.655053
## 68   Male 75 - 79            TEQ To child 98.298013
## 69 Female 80 - 84            TEQ To child 72.655053
## 70   Male 80 - 84            TEQ To child 98.298013
## 71 Female     85+            TEQ To child 72.655053
## 72   Male     85+            TEQ To child 98.298013
oprint(summary(exposure,"mean"))
##     Gender     Age Exposure_agent Exposure        mean
## 1   Female   0 - 4            ALA   Direct  21.4337074
## 2     Male   0 - 4            ALA   Direct  28.9985453
## 3   Female   5 - 9            ALA   Direct  39.6865892
## 4     Male   5 - 9            ALA   Direct  53.6936207
## 5   Female 10 - 14            ALA   Direct  60.3859645
## 6     Male 10 - 14            ALA   Direct  81.6986578
## 7   Female 15 - 19            ALA   Direct  60.3859645
## 8     Male 15 - 19            ALA   Direct  81.6986578
## 9   Female 20 - 24            ALA   Direct  99.7525088
## 10    Male 20 - 24            ALA   Direct 134.9592767
## 11  Female 25 - 29            ALA   Direct  99.7525088
## 12    Male 25 - 29            ALA   Direct 134.9592767
## 13  Female 30 - 34            ALA   Direct  99.7525088
## 14    Male 30 - 34            ALA   Direct 134.9592767
## 15  Female 35 - 39            ALA   Direct  99.7525088
## 16    Male 35 - 39            ALA   Direct 134.9592767
## 17  Female 40 - 44            ALA   Direct  99.7525088
## 18    Male 40 - 44            ALA   Direct 134.9592767
## 19  Female 45 - 49            ALA   Direct  99.7525088
## 20    Male 45 - 49            ALA   Direct 134.9592767
## 21  Female 50 - 54            ALA   Direct  99.7525088
## 22    Male 50 - 54            ALA   Direct 134.9592767
## 23  Female 55 - 59            ALA   Direct  99.7525088
## 24    Male 55 - 59            ALA   Direct 134.9592767
## 25  Female 60 - 64            ALA   Direct  99.7525088
## 26    Male 60 - 64            ALA   Direct 134.9592767
## 27  Female 65 - 69            ALA   Direct 138.5860138
## 28    Male 65 - 69            ALA   Direct 187.4987246
## 29  Female 70 - 74            ALA   Direct 138.5860138
## 30    Male 70 - 74            ALA   Direct 187.4987246
## 31  Female 75 - 79            ALA   Direct 138.5860138
## 32    Male 75 - 79            ALA   Direct 187.4987246
## 33  Female 80 - 84            ALA   Direct 138.5860138
## 34    Male 80 - 84            ALA   Direct 187.4987246
## 35  Female     85+            ALA   Direct 138.5860138
## 36    Male     85+            ALA   Direct 187.4987246
## 37  Female   0 - 4            DHA   Direct  55.5313401
## 38    Male   0 - 4            DHA   Direct  76.5297273
## 39  Female   5 - 9            DHA   Direct  91.8105802
## 40    Male   5 - 9            DHA   Direct 125.6134050
## 41  Female 10 - 14            DHA   Direct 103.9070801
## 42    Male 10 - 14            DHA   Direct 141.9792579
## 43  Female 15 - 19            DHA   Direct 103.9070801
## 44    Male 15 - 19            DHA   Direct 141.9792579
## 45  Female 20 - 24            DHA   Direct 150.9747434
## 46    Male 20 - 24            DHA   Direct 205.6590377
## 47  Female 25 - 29            DHA   Direct 150.9747434
## 48    Male 25 - 29            DHA   Direct 205.6590377
## 49  Female 30 - 34            DHA   Direct 150.9747434
## 50    Male 30 - 34            DHA   Direct 205.6590377
## 51  Female 35 - 39            DHA   Direct 150.9747434
## 52    Male 35 - 39            DHA   Direct 205.6590377
## 53  Female 40 - 44            DHA   Direct 150.9747434
## 54    Male 40 - 44            DHA   Direct 205.6590377
## 55  Female 45 - 49            DHA   Direct 150.9747434
## 56    Male 45 - 49            DHA   Direct 205.6590377
## 57  Female 50 - 54            DHA   Direct 150.9747434
## 58    Male 50 - 54            DHA   Direct 205.6590377
## 59  Female 55 - 59            DHA   Direct 150.9747434
## 60    Male 55 - 59            DHA   Direct 205.6590377
## 61  Female 60 - 64            DHA   Direct 150.9747434
## 62    Male 60 - 64            DHA   Direct 205.6590377
## 63  Female 65 - 69            DHA   Direct 194.8920588
## 64    Male 65 - 69            DHA   Direct 265.0765820
## 65  Female 70 - 74            DHA   Direct 194.8920588
## 66    Male 70 - 74            DHA   Direct 265.0765820
## 67  Female 75 - 79            DHA   Direct 194.8920588
## 68    Male 75 - 79            DHA   Direct 265.0765820
## 69  Female 80 - 84            DHA   Direct 194.8920588
## 70    Male 80 - 84            DHA   Direct 265.0765820
## 71  Female     85+            DHA   Direct 194.8920588
## 72    Male     85+            DHA   Direct 265.0765820
## 73  Female   0 - 4            EPA   Direct  29.7100701
## 74    Male   0 - 4            EPA   Direct  43.4721516
## 75  Female   5 - 9            EPA   Direct  29.7100701
## 76    Male   5 - 9            EPA   Direct  43.4721516
## 77  Female 10 - 14            EPA   Direct  29.7100701
## 78    Male 10 - 14            EPA   Direct  43.4721516
## 79  Female 15 - 19            EPA   Direct  29.7100701
## 80    Male 15 - 19            EPA   Direct  43.4721516
## 81  Female 20 - 24            EPA   Direct  29.7100701
## 82    Male 20 - 24            EPA   Direct  43.4721516
## 83  Female 25 - 29            EPA   Direct  29.7100701
## 84    Male 25 - 29            EPA   Direct  43.4721516
## 85  Female 30 - 34            EPA   Direct  29.7100701
## 86    Male 30 - 34            EPA   Direct  43.4721516
## 87  Female 35 - 39            EPA   Direct  29.7100701
## 88    Male 35 - 39            EPA   Direct  43.4721516
## 89  Female 40 - 44            EPA   Direct  29.7100701
## 90    Male 40 - 44            EPA   Direct  43.4721516
## 91  Female 45 - 49            EPA   Direct  29.7100701
## 92    Male 45 - 49            EPA   Direct  43.4721516
## 93  Female 50 - 54            EPA   Direct  29.7100701
## 94    Male 50 - 54            EPA   Direct  43.4721516
## 95  Female 55 - 59            EPA   Direct  29.7100701
## 96    Male 55 - 59            EPA   Direct  43.4721516
## 97  Female 60 - 64            EPA   Direct  29.7100701
## 98    Male 60 - 64            EPA   Direct  43.4721516
## 99  Female 65 - 69            EPA   Direct  29.7100701
## 100   Male 65 - 69            EPA   Direct  43.4721516
## 101 Female 70 - 74            EPA   Direct  29.7100701
## 102   Male 70 - 74            EPA   Direct  43.4721516
## 103 Female 75 - 79            EPA   Direct  29.7100701
## 104   Male 75 - 79            EPA   Direct  43.4721516
## 105 Female 80 - 84            EPA   Direct  29.7100701
## 106   Male 80 - 84            EPA   Direct  43.4721516
## 107 Female     85+            EPA   Direct  29.7100701
## 108   Male     85+            EPA   Direct  43.4721516
## 109 Female   0 - 4           Fish   Direct   5.7621802
## 110   Male   0 - 4           Fish   Direct   7.7958909
## 111 Female   5 - 9           Fish   Direct  17.5698850
## 112   Male   5 - 9           Fish   Direct  23.7710209
## 113 Female 10 - 14           Fish   Direct  17.6141458
## 114   Male 10 - 14           Fish   Direct  23.8309032
## 115 Female 15 - 19           Fish   Direct  17.6141458
## 116   Male 15 - 19           Fish   Direct  23.8309032
## 117 Female 20 - 24           Fish   Direct  27.6701889
## 118   Male 20 - 24           Fish   Direct  37.4361379
## 119 Female 25 - 29           Fish   Direct  27.6701889
## 120   Male 25 - 29           Fish   Direct  37.4361379
## 121 Female 30 - 34           Fish   Direct  27.6701889
## 122   Male 30 - 34           Fish   Direct  37.4361379
## 123 Female 35 - 39           Fish   Direct  27.6701889
## 124   Male 35 - 39           Fish   Direct  37.4361379
## 125 Female 40 - 44           Fish   Direct  27.6701889
## 126   Male 40 - 44           Fish   Direct  37.4361379
## 127 Female 45 - 49           Fish   Direct  27.6701889
## 128   Male 45 - 49           Fish   Direct  37.4361379
## 129 Female 50 - 54           Fish   Direct  27.6701889
## 130   Male 50 - 54           Fish   Direct  37.4361379
## 131 Female 55 - 59           Fish   Direct  27.6701889
## 132   Male 55 - 59           Fish   Direct  37.4361379
## 133 Female 60 - 64           Fish   Direct  27.6701889
## 134   Male 60 - 64           Fish   Direct  37.4361379
## 135 Female 65 - 69           Fish   Direct  35.7717063
## 136   Male 65 - 69           Fish   Direct  48.3970145
## 137 Female 70 - 74           Fish   Direct  35.7717063
## 138   Male 70 - 74           Fish   Direct  48.3970145
## 139 Female 75 - 79           Fish   Direct  35.7717063
## 140   Male 75 - 79           Fish   Direct  48.3970145
## 141 Female 80 - 84           Fish   Direct  35.7717063
## 142   Male 80 - 84           Fish   Direct  48.3970145
## 143 Female     85+           Fish   Direct  35.7717063
## 144   Male     85+           Fish   Direct  48.3970145
## 145 Female   0 - 4           MeHg   Direct   0.5971765
## 146   Male   0 - 4           MeHg   Direct   0.8079447
## 147 Female   5 - 9           MeHg   Direct   1.3270854
## 148   Male   5 - 9           MeHg   Direct   1.7954684
## 149 Female 10 - 14           MeHg   Direct   1.6315120
## 150   Male 10 - 14           MeHg   Direct   2.2073398
## 151 Female 15 - 19           MeHg   Direct   1.6315120
## 152   Male 15 - 19           MeHg   Direct   2.2073398
## 153 Female 20 - 24           MeHg   Direct   2.8284141
## 154   Male 20 - 24           MeHg   Direct   3.8266780
## 155 Female 25 - 29           MeHg   Direct   2.8284141
## 156   Male 25 - 29           MeHg   Direct   3.8266780
## 157 Female 30 - 34           MeHg   Direct   2.8284141
## 158   Male 30 - 34           MeHg   Direct   3.8266780
## 159 Female 35 - 39           MeHg   Direct   2.8284141
## 160   Male 35 - 39           MeHg   Direct   3.8266780
## 161 Female 40 - 44           MeHg   Direct   2.8284141
## 162   Male 40 - 44           MeHg   Direct   3.8266780
## 163 Female 45 - 49           MeHg   Direct   2.8284141
## 164   Male 45 - 49           MeHg   Direct   3.8266780
## 165 Female 50 - 54           MeHg   Direct   2.8284141
## 166   Male 50 - 54           MeHg   Direct   3.8266780
## 167 Female 55 - 59           MeHg   Direct   2.8284141
## 168   Male 55 - 59           MeHg   Direct   3.8266780
## 169 Female 60 - 64           MeHg   Direct   2.8284141
## 170   Male 60 - 64           MeHg   Direct   3.8266780
## 171 Female 65 - 69           MeHg   Direct   3.7425454
## 172   Male 65 - 69           MeHg   Direct   5.0634437
## 173 Female 70 - 74           MeHg   Direct   3.7425454
## 174   Male 70 - 74           MeHg   Direct   5.0634437
## 175 Female 75 - 79           MeHg   Direct   3.7425454
## 176   Male 75 - 79           MeHg   Direct   5.0634437
## 177 Female 80 - 84           MeHg   Direct   3.7425454
## 178   Male 80 - 84           MeHg   Direct   5.0634437
## 179 Female     85+           MeHg   Direct   3.7425454
## 180   Male     85+           MeHg   Direct   5.0634437
## 181 Female   0 - 4         Omega3   Direct 141.4674764
## 182   Male   0 - 4         Omega3   Direct 192.7261932
## 183 Female   5 - 9         Omega3   Direct 248.8196536
## 184   Male   5 - 9         Omega3   Direct 337.9673743
## 185 Female 10 - 14         Omega3   Direct 295.6121154
## 186   Male 10 - 14         Omega3   Direct 401.2748226
## 187 Female 15 - 19         Omega3   Direct 295.6121154
## 188   Male 15 - 19         Omega3   Direct 401.2748226
## 189 Female 20 - 24         Omega3   Direct 443.9103092
## 190   Male 20 - 24         Omega3   Direct 601.9135554
## 191 Female 25 - 29         Omega3   Direct 443.9103092
## 192   Male 25 - 29         Omega3   Direct 601.9135554
## 193 Female 30 - 34         Omega3   Direct 443.9103092
## 194   Male 30 - 34         Omega3   Direct 601.9135554
## 195 Female 35 - 39         Omega3   Direct 443.9103092
## 196   Male 35 - 39         Omega3   Direct 601.9135554
## 197 Female 40 - 44         Omega3   Direct 443.9103092
## 198   Male 40 - 44         Omega3   Direct 601.9135554
## 199 Female 45 - 49         Omega3   Direct 443.9103092
## 200   Male 45 - 49         Omega3   Direct 601.9135554
## 201 Female 50 - 54         Omega3   Direct 443.9103092
## 202   Male 50 - 54         Omega3   Direct 601.9135554
## 203 Female 55 - 59         Omega3   Direct 443.9103092
## 204   Male 55 - 59         Omega3   Direct 601.9135554
## 205 Female 60 - 64         Omega3   Direct 443.9103092
## 206   Male 60 - 64         Omega3   Direct 601.9135554
## 207 Female 65 - 69         Omega3   Direct 599.1340842
## 208   Male 65 - 69         Omega3   Direct 811.9221921
## 209 Female 70 - 74         Omega3   Direct 599.1340842
## 210   Male 70 - 74         Omega3   Direct 811.9221921
## 211 Female 75 - 79         Omega3   Direct 599.1340842
## 212   Male 75 - 79         Omega3   Direct 811.9221921
## 213 Female 80 - 84         Omega3   Direct 599.1340842
## 214   Male 80 - 84         Omega3   Direct 811.9221921
## 215 Female     85+         Omega3   Direct 599.1340842
## 216   Male     85+         Omega3   Direct 811.9221921
## 217 Female   0 - 4           PFAS   Direct  47.2495764
## 218   Male   0 - 4           PFAS   Direct  63.9258974
## 219 Female   5 - 9           PFAS   Direct 149.4611117
## 220   Male   5 - 9           PFAS   Direct 202.2120924
## 221 Female 10 - 14           PFAS   Direct 151.7132151
## 222   Male 10 - 14           PFAS   Direct 205.2590557
## 223 Female 15 - 19           PFAS   Direct 151.7132151
## 224   Male 15 - 19           PFAS   Direct 205.2590557
## 225 Female 20 - 24           PFAS   Direct 230.7378180
## 226   Male 20 - 24           PFAS   Direct 312.1746949
## 227 Female 25 - 29           PFAS   Direct 230.7378180
## 228   Male 25 - 29           PFAS   Direct 312.1746949
## 229 Female 30 - 34           PFAS   Direct 230.7378180
## 230   Male 30 - 34           PFAS   Direct 312.1746949
## 231 Female 35 - 39           PFAS   Direct 230.7378180
## 232   Male 35 - 39           PFAS   Direct 312.1746949
## 233 Female 40 - 44           PFAS   Direct 230.7378180
## 234   Male 40 - 44           PFAS   Direct 312.1746949
## 235 Female 45 - 49           PFAS   Direct 230.7378180
## 236   Male 45 - 49           PFAS   Direct 312.1746949
## 237 Female 50 - 54           PFAS   Direct 230.7378180
## 238   Male 50 - 54           PFAS   Direct 312.1746949
## 239 Female 55 - 59           PFAS   Direct 230.7378180
## 240   Male 55 - 59           PFAS   Direct 312.1746949
## 241 Female 60 - 64           PFAS   Direct 230.7378180
## 242   Male 60 - 64           PFAS   Direct 312.1746949
## 243 Female 65 - 69           PFAS   Direct 308.7619762
## 244   Male 65 - 69           PFAS   Direct 417.7367914
## 245 Female 70 - 74           PFAS   Direct 308.7619762
## 246   Male 70 - 74           PFAS   Direct 417.7367914
## 247 Female 75 - 79           PFAS   Direct 308.7619762
## 248   Male 75 - 79           PFAS   Direct 417.7367914
## 249 Female 80 - 84           PFAS   Direct 308.7619762
## 250   Male 80 - 84           PFAS   Direct 417.7367914
## 251 Female     85+           PFAS   Direct 308.7619762
## 252   Male     85+           PFAS   Direct 417.7367914
## 253 Female   0 - 4            TEQ   Direct  54.2533243
## 254   Male   0 - 4            TEQ   Direct  73.4015564
## 255 Female   5 - 9            TEQ   Direct  97.8208278
## 256   Male   5 - 9            TEQ   Direct 132.3458259
## 257 Female 10 - 14            TEQ   Direct 160.1890499
## 258   Male 10 - 14            TEQ   Direct 216.7263616
## 259 Female 15 - 19            TEQ   Direct 160.1890499
## 260   Male 15 - 19            TEQ   Direct 216.7263616
## 261 Female 20 - 24            TEQ   Direct 301.8922163
## 262   Male 20 - 24            TEQ   Direct 408.4424103
## 263 Female 25 - 29            TEQ   Direct 301.8922163
## 264   Male 25 - 29            TEQ   Direct 408.4424103
## 265 Female 30 - 34            TEQ   Direct 301.8922163
## 266   Male 30 - 34            TEQ   Direct 408.4424103
## 267 Female 35 - 39            TEQ   Direct 301.8922163
## 268   Male 35 - 39            TEQ   Direct 408.4424103
## 269 Female 40 - 44            TEQ   Direct 301.8922163
## 270   Male 40 - 44            TEQ   Direct 408.4424103
## 271 Female 45 - 49            TEQ   Direct 301.8922163
## 272   Male 45 - 49            TEQ   Direct 408.4424103
## 273 Female 50 - 54            TEQ   Direct 301.8922163
## 274   Male 50 - 54            TEQ   Direct 408.4424103
## 275 Female 55 - 59            TEQ   Direct 301.8922163
## 276   Male 55 - 59            TEQ   Direct 408.4424103
## 277 Female 60 - 64            TEQ   Direct 301.8922163
## 278   Male 60 - 64            TEQ   Direct 408.4424103
## 279 Female 65 - 69            TEQ   Direct 541.2734024
## 280   Male 65 - 69            TEQ   Direct 732.3110739
## 281 Female 70 - 74            TEQ   Direct 541.2734024
## 282   Male 70 - 74            TEQ   Direct 732.3110739
## 283 Female 75 - 79            TEQ   Direct 541.2734024
## 284   Male 75 - 79            TEQ   Direct 732.3110739
## 285 Female 80 - 84            TEQ   Direct 541.2734024
## 286   Male 80 - 84            TEQ   Direct 732.3110739
## 287 Female     85+            TEQ   Direct 541.2734024
## 288   Male     85+            TEQ   Direct 732.3110739
## 289 Female   0 - 4      Vitamin D   Direct   3.3655848
## 290   Male   0 - 4      Vitamin D   Direct   5.1938819
## 291 Female   5 - 9      Vitamin D   Direct   4.5671267
## 292   Male   5 - 9      Vitamin D   Direct   6.8194974
## 293 Female 10 - 14      Vitamin D   Direct   4.4677666
## 294   Male 10 - 14      Vitamin D   Direct   6.6850691
## 295 Female 15 - 19      Vitamin D   Direct   4.4677666
## 296   Male 15 - 19      Vitamin D   Direct   6.6850691
## 297 Female 20 - 24      Vitamin D   Direct   5.5261030
## 298   Male 20 - 24      Vitamin D   Direct   8.1169360
## 299 Female 25 - 29      Vitamin D   Direct   5.5261030
## 300   Male 25 - 29      Vitamin D   Direct   8.1169360
## 301 Female 30 - 34      Vitamin D   Direct   5.5261030
## 302   Male 30 - 34      Vitamin D   Direct   8.1169360
## 303 Female 35 - 39      Vitamin D   Direct   5.5261030
## 304   Male 35 - 39      Vitamin D   Direct   8.1169360
## 305 Female 40 - 44      Vitamin D   Direct   5.5261030
## 306   Male 40 - 44      Vitamin D   Direct   8.1169360
## 307 Female 45 - 49      Vitamin D   Direct   5.5261030
## 308   Male 45 - 49      Vitamin D   Direct   8.1169360
## 309 Female 50 - 54      Vitamin D   Direct   5.5261030
## 310   Male 50 - 54      Vitamin D   Direct   8.1169360
## 311 Female 55 - 59      Vitamin D   Direct   5.5261030
## 312   Male 55 - 59      Vitamin D   Direct   8.1169360
## 313 Female 60 - 64      Vitamin D   Direct   5.5261030
## 314   Male 60 - 64      Vitamin D   Direct   8.1169360
## 315 Female 65 - 69      Vitamin D   Direct   6.2731073
## 316   Male 65 - 69      Vitamin D   Direct   9.1275888
## 317 Female 70 - 74      Vitamin D   Direct   6.2731073
## 318   Male 70 - 74      Vitamin D   Direct   9.1275888
## 319 Female 75 - 79      Vitamin D   Direct   6.2731073
## 320   Male 75 - 79      Vitamin D   Direct   9.1275888
## 321 Female 80 - 84      Vitamin D   Direct   6.2731073
## 322   Male 80 - 84      Vitamin D   Direct   9.1275888
## 323 Female     85+      Vitamin D   Direct   6.2731073
## 324   Male     85+      Vitamin D   Direct   9.1275888
## 325 Female   0 - 4           PFAS To child  13.4613861
## 326   Male   0 - 4           PFAS To child  13.4613861
## 327 Female   0 - 4            TEQ To child  40.3959687
## 328   Male   0 - 4            TEQ To child  40.3959687
#oprint(summary(fish_proportion,"mean"))
oprint(summary(incidence,"mean"))
##                    Response   Age Adjust  mean
## 1         Immunosuppression 0 - 4    BAU 4.810
## 2 Loss in child's IQ points 0 - 4    BAU 5.960
## 3       Sperm concentration 0 - 4    BAU 0.070
## 4   Yes or no dental defect 0 - 4    BAU 0.224
## 5         Immunosuppression 5 - 9    BAU 3.950
oprint(summary(PAF,"mean"))
##     Exposure_agent                                          Response Gender
## 1              DHA                         Loss in child's IQ points Female
## 2             MeHg                         Loss in child's IQ points Female
## 3              TEQ                           Cancer morbidity yearly Female
## 4              TEQ                               Sperm concentration Female
## 5              TEQ                           Yes or no dental defect Female
## 6             PFAS                                 Immunosuppression Female
## 7              TEQ Dioxin recommendation tolerable daily intake 2018 Female
## 8        Vitamin D                          Vitamin D recommendation Female
## 9             PFAS                                          PFAS TWI Female
## 10            Fish                               All-cause mortality Female
## 11          Omega3                                     Breast cancer Female
## 12            Fish                                        Depression Female
## 13          Omega3                                    CHD2 mortality Female
## 14             DHA                         Loss in child's IQ points   Male
## 15            MeHg                         Loss in child's IQ points   Male
## 16             TEQ                           Cancer morbidity yearly   Male
## 17             TEQ                               Sperm concentration   Male
## 18             TEQ                           Yes or no dental defect   Male
## 19            PFAS                                 Immunosuppression   Male
## 20             TEQ Dioxin recommendation tolerable daily intake 2018   Male
## 21       Vitamin D                          Vitamin D recommendation   Male
## 22            PFAS                                          PFAS TWI   Male
## 23            Fish                               All-cause mortality   Male
## 24          Omega3                                     Breast cancer   Male
## 25            Fish                                        Depression   Male
## 26          Omega3                                    CHD2 mortality   Male
## 27             TEQ                           Cancer morbidity yearly Female
## 28            PFAS                                 Immunosuppression Female
## 29             TEQ Dioxin recommendation tolerable daily intake 2018 Female
## 30       Vitamin D                          Vitamin D recommendation Female
## 31            PFAS                                          PFAS TWI Female
## 32            Fish                               All-cause mortality Female
## 33          Omega3                                     Breast cancer Female
## 34            Fish                                        Depression Female
## 35          Omega3                                    CHD2 mortality Female
## 36             TEQ                           Cancer morbidity yearly   Male
## 37            PFAS                                 Immunosuppression   Male
## 38             TEQ Dioxin recommendation tolerable daily intake 2018   Male
## 39       Vitamin D                          Vitamin D recommendation   Male
## 40            PFAS                                          PFAS TWI   Male
## 41            Fish                               All-cause mortality   Male
## 42          Omega3                                     Breast cancer   Male
## 43            Fish                                        Depression   Male
## 44          Omega3                                    CHD2 mortality   Male
## 45             TEQ                           Cancer morbidity yearly Female
## 46             TEQ Dioxin recommendation tolerable daily intake 2018 Female
## 47       Vitamin D                          Vitamin D recommendation Female
## 48            PFAS                                          PFAS TWI Female
## 49            Fish                               All-cause mortality Female
## 50          Omega3                                     Breast cancer Female
## 51            Fish                                        Depression Female
## 52          Omega3                                    CHD2 mortality Female
## 53             TEQ                           Cancer morbidity yearly   Male
## 54             TEQ Dioxin recommendation tolerable daily intake 2018   Male
## 55       Vitamin D                          Vitamin D recommendation   Male
## 56            PFAS                                          PFAS TWI   Male
## 57            Fish                               All-cause mortality   Male
## 58          Omega3                                     Breast cancer   Male
## 59            Fish                                        Depression   Male
## 60          Omega3                                    CHD2 mortality   Male
## 61             TEQ                           Cancer morbidity yearly Female
## 62             TEQ Dioxin recommendation tolerable daily intake 2018 Female
## 63       Vitamin D                          Vitamin D recommendation Female
## 64            PFAS                                          PFAS TWI Female
## 65            Fish                               All-cause mortality Female
## 66          Omega3                                     Breast cancer Female
## 67            Fish                                        Depression Female
## 68          Omega3                                    CHD2 mortality Female
## 69             TEQ                           Cancer morbidity yearly   Male
## 70             TEQ Dioxin recommendation tolerable daily intake 2018   Male
## 71       Vitamin D                          Vitamin D recommendation   Male
## 72            PFAS                                          PFAS TWI   Male
## 73            Fish                               All-cause mortality   Male
## 74          Omega3                                     Breast cancer   Male
## 75            Fish                                        Depression   Male
## 76          Omega3                                    CHD2 mortality   Male
## 77             TEQ                           Cancer morbidity yearly Female
## 78             TEQ Dioxin recommendation tolerable daily intake 2018 Female
## 79       Vitamin D                          Vitamin D recommendation Female
## 80            PFAS                                          PFAS TWI Female
## 81            Fish                               All-cause mortality Female
## 82          Omega3                                     Breast cancer Female
## 83            Fish                                        Depression Female
## 84          Omega3                                    CHD2 mortality Female
## 85             TEQ                           Cancer morbidity yearly   Male
## 86             TEQ Dioxin recommendation tolerable daily intake 2018   Male
## 87       Vitamin D                          Vitamin D recommendation   Male
## 88            PFAS                                          PFAS TWI   Male
## 89            Fish                               All-cause mortality   Male
## 90          Omega3                                     Breast cancer   Male
## 91            Fish                                        Depression   Male
## 92          Omega3                                    CHD2 mortality   Male
## 93             TEQ                           Cancer morbidity yearly Female
## 94             TEQ Dioxin recommendation tolerable daily intake 2018 Female
## 95       Vitamin D                          Vitamin D recommendation Female
## 96            PFAS                                          PFAS TWI Female
## 97            Fish                               All-cause mortality Female
## 98          Omega3                                     Breast cancer Female
## 99            Fish                                        Depression Female
## 100         Omega3                                    CHD2 mortality Female
## 101            TEQ                           Cancer morbidity yearly   Male
## 102            TEQ Dioxin recommendation tolerable daily intake 2018   Male
## 103      Vitamin D                          Vitamin D recommendation   Male
## 104           PFAS                                          PFAS TWI   Male
## 105           Fish                               All-cause mortality   Male
## 106         Omega3                                     Breast cancer   Male
## 107           Fish                                        Depression   Male
## 108         Omega3                                    CHD2 mortality   Male
## 109            TEQ                           Cancer morbidity yearly Female
## 110            TEQ Dioxin recommendation tolerable daily intake 2018 Female
## 111      Vitamin D                          Vitamin D recommendation Female
## 112           PFAS                                          PFAS TWI Female
## 113           Fish                               All-cause mortality Female
## 114         Omega3                                     Breast cancer Female
## 115           Fish                                        Depression Female
## 116         Omega3                                    CHD2 mortality Female
## 117            TEQ                           Cancer morbidity yearly   Male
## 118            TEQ Dioxin recommendation tolerable daily intake 2018   Male
## 119      Vitamin D                          Vitamin D recommendation   Male
## 120           PFAS                                          PFAS TWI   Male
## 121           Fish                               All-cause mortality   Male
## 122         Omega3                                     Breast cancer   Male
## 123           Fish                                        Depression   Male
## 124         Omega3                                    CHD2 mortality   Male
## 125            TEQ                           Cancer morbidity yearly Female
## 126            TEQ Dioxin recommendation tolerable daily intake 2018 Female
## 127      Vitamin D                          Vitamin D recommendation Female
## 128           PFAS                                          PFAS TWI Female
## 129           Fish                               All-cause mortality Female
## 130         Omega3                                     Breast cancer Female
## 131           Fish                                        Depression Female
## 132         Omega3                                    CHD2 mortality Female
## 133            TEQ                           Cancer morbidity yearly   Male
## 134            TEQ Dioxin recommendation tolerable daily intake 2018   Male
## 135      Vitamin D                          Vitamin D recommendation   Male
## 136           PFAS                                          PFAS TWI   Male
## 137           Fish                               All-cause mortality   Male
## 138         Omega3                                     Breast cancer   Male
## 139           Fish                                        Depression   Male
## 140         Omega3                                    CHD2 mortality   Male
## 141            TEQ                           Cancer morbidity yearly Female
## 142            TEQ Dioxin recommendation tolerable daily intake 2018 Female
## 143      Vitamin D                          Vitamin D recommendation Female
## 144           PFAS                                          PFAS TWI Female
## 145           Fish                               All-cause mortality Female
## 146         Omega3                                     Breast cancer Female
## 147           Fish                                        Depression Female
## 148         Omega3                                    CHD2 mortality Female
## 149            TEQ                           Cancer morbidity yearly   Male
## 150            TEQ Dioxin recommendation tolerable daily intake 2018   Male
## 151      Vitamin D                          Vitamin D recommendation   Male
## 152           PFAS                                          PFAS TWI   Male
## 153           Fish                               All-cause mortality   Male
## 154         Omega3                                     Breast cancer   Male
## 155           Fish                                        Depression   Male
## 156         Omega3                                    CHD2 mortality   Male
## 157            TEQ                           Cancer morbidity yearly Female
## 158            TEQ Dioxin recommendation tolerable daily intake 2018 Female
## 159      Vitamin D                          Vitamin D recommendation Female
## 160           PFAS                                          PFAS TWI Female
## 161           Fish                               All-cause mortality Female
## 162         Omega3                                     Breast cancer Female
## 163           Fish                                        Depression Female
## 164         Omega3                                    CHD2 mortality Female
## 165            TEQ                           Cancer morbidity yearly   Male
## 166            TEQ Dioxin recommendation tolerable daily intake 2018   Male
## 167      Vitamin D                          Vitamin D recommendation   Male
## 168           PFAS                                          PFAS TWI   Male
## 169           Fish                               All-cause mortality   Male
## 170         Omega3                                     Breast cancer   Male
## 171           Fish                                        Depression   Male
## 172         Omega3                                    CHD2 mortality   Male
## 173            TEQ                           Cancer morbidity yearly Female
## 174            TEQ Dioxin recommendation tolerable daily intake 2018 Female
## 175      Vitamin D                          Vitamin D recommendation Female
## 176           PFAS                                          PFAS TWI Female
## 177           Fish                               All-cause mortality Female
## 178         Omega3                                     Breast cancer Female
## 179           Fish                                        Depression Female
## 180         Omega3                                    CHD2 mortality Female
## 181            TEQ                           Cancer morbidity yearly   Male
## 182            TEQ Dioxin recommendation tolerable daily intake 2018   Male
## 183      Vitamin D                          Vitamin D recommendation   Male
## 184           PFAS                                          PFAS TWI   Male
## 185           Fish                               All-cause mortality   Male
## 186         Omega3                                     Breast cancer   Male
## 187           Fish                                        Depression   Male
## 188         Omega3                                    CHD2 mortality   Male
## 189            TEQ                           Cancer morbidity yearly Female
## 190            TEQ Dioxin recommendation tolerable daily intake 2018 Female
## 191      Vitamin D                          Vitamin D recommendation Female
## 192           PFAS                                          PFAS TWI Female
## 193           Fish                               All-cause mortality Female
## 194         Omega3                                     Breast cancer Female
## 195           Fish                                        Depression Female
## 196         Omega3                                    CHD2 mortality Female
## 197            TEQ                           Cancer morbidity yearly   Male
## 198            TEQ Dioxin recommendation tolerable daily intake 2018   Male
## 199      Vitamin D                          Vitamin D recommendation   Male
## 200           PFAS                                          PFAS TWI   Male
## 201           Fish                               All-cause mortality   Male
## 202         Omega3                                     Breast cancer   Male
## 203           Fish                                        Depression   Male
## 204         Omega3                                    CHD2 mortality   Male
## 205            TEQ                           Cancer morbidity yearly Female
## 206            TEQ Dioxin recommendation tolerable daily intake 2018 Female
## 207      Vitamin D                          Vitamin D recommendation Female
## 208           PFAS                                          PFAS TWI Female
## 209           Fish                               All-cause mortality Female
## 210         Omega3                                     Breast cancer Female
## 211           Fish                                        Depression Female
## 212         Omega3                                    CHD2 mortality Female
## 213            TEQ                           Cancer morbidity yearly   Male
## 214            TEQ Dioxin recommendation tolerable daily intake 2018   Male
## 215      Vitamin D                          Vitamin D recommendation   Male
## 216           PFAS                                          PFAS TWI   Male
## 217           Fish                               All-cause mortality   Male
## 218         Omega3                                     Breast cancer   Male
## 219           Fish                                        Depression   Male
## 220         Omega3                                    CHD2 mortality   Male
## 221            TEQ                           Cancer morbidity yearly Female
## 222            TEQ Dioxin recommendation tolerable daily intake 2018 Female
## 223      Vitamin D                          Vitamin D recommendation Female
## 224           PFAS                                          PFAS TWI Female
## 225           Fish                               All-cause mortality Female
## 226         Omega3                                     Breast cancer Female
## 227           Fish                                        Depression Female
## 228         Omega3                                    CHD2 mortality Female
## 229            TEQ                           Cancer morbidity yearly   Male
## 230            TEQ Dioxin recommendation tolerable daily intake 2018   Male
## 231      Vitamin D                          Vitamin D recommendation   Male
## 232           PFAS                                          PFAS TWI   Male
## 233           Fish                               All-cause mortality   Male
## 234         Omega3                                     Breast cancer   Male
## 235           Fish                                        Depression   Male
## 236         Omega3                                    CHD2 mortality   Male
## 237            TEQ                           Cancer morbidity yearly Female
## 238            TEQ Dioxin recommendation tolerable daily intake 2018 Female
## 239      Vitamin D                          Vitamin D recommendation Female
## 240           PFAS                                          PFAS TWI Female
## 241           Fish                               All-cause mortality Female
## 242         Omega3                                     Breast cancer Female
## 243           Fish                                        Depression Female
## 244         Omega3                                    CHD2 mortality Female
## 245            TEQ                           Cancer morbidity yearly   Male
## 246            TEQ Dioxin recommendation tolerable daily intake 2018   Male
## 247      Vitamin D                          Vitamin D recommendation   Male
## 248           PFAS                                          PFAS TWI   Male
## 249           Fish                               All-cause mortality   Male
## 250         Omega3                                     Breast cancer   Male
## 251           Fish                                        Depression   Male
## 252         Omega3                                    CHD2 mortality   Male
## 253            TEQ                           Cancer morbidity yearly Female
## 254            TEQ Dioxin recommendation tolerable daily intake 2018 Female
## 255      Vitamin D                          Vitamin D recommendation Female
## 256           PFAS                                          PFAS TWI Female
## 257           Fish                               All-cause mortality Female
## 258         Omega3                                     Breast cancer Female
## 259           Fish                                        Depression Female
## 260         Omega3                                    CHD2 mortality Female
## 261            TEQ                           Cancer morbidity yearly   Male
## 262            TEQ Dioxin recommendation tolerable daily intake 2018   Male
## 263      Vitamin D                          Vitamin D recommendation   Male
## 264           PFAS                                          PFAS TWI   Male
## 265           Fish                               All-cause mortality   Male
## 266         Omega3                                     Breast cancer   Male
## 267           Fish                                        Depression   Male
## 268         Omega3                                    CHD2 mortality   Male
## 269            TEQ                           Cancer morbidity yearly Female
## 270            TEQ Dioxin recommendation tolerable daily intake 2018 Female
## 271      Vitamin D                          Vitamin D recommendation Female
## 272           PFAS                                          PFAS TWI Female
## 273           Fish                               All-cause mortality Female
## 274         Omega3                                     Breast cancer Female
## 275           Fish                                        Depression Female
## 276         Omega3                                    CHD2 mortality Female
## 277            TEQ                           Cancer morbidity yearly   Male
## 278            TEQ Dioxin recommendation tolerable daily intake 2018   Male
## 279      Vitamin D                          Vitamin D recommendation   Male
## 280           PFAS                                          PFAS TWI   Male
## 281           Fish                               All-cause mortality   Male
## 282         Omega3                                     Breast cancer   Male
## 283           Fish                                        Depression   Male
## 284         Omega3                                    CHD2 mortality   Male
## 285            TEQ                           Cancer morbidity yearly Female
## 286            TEQ Dioxin recommendation tolerable daily intake 2018 Female
## 287      Vitamin D                          Vitamin D recommendation Female
## 288           PFAS                                          PFAS TWI Female
## 289           Fish                               All-cause mortality Female
## 290         Omega3                                     Breast cancer Female
## 291           Fish                                        Depression Female
## 292         Omega3                                    CHD2 mortality Female
## 293            TEQ                           Cancer morbidity yearly   Male
## 294            TEQ Dioxin recommendation tolerable daily intake 2018   Male
## 295      Vitamin D                          Vitamin D recommendation   Male
## 296           PFAS                                          PFAS TWI   Male
## 297           Fish                               All-cause mortality   Male
## 298         Omega3                                     Breast cancer   Male
## 299           Fish                                        Depression   Male
## 300         Omega3                                    CHD2 mortality   Male
##         Age Adjust         mean
## 1     0 - 4    BAU -0.012266039
## 2     0 - 4    BAU  0.013555654
## 3     0 - 4    BAU  0.002022785
## 4     0 - 4    BAU  0.080238634
## 5     0 - 4    BAU  0.599425454
## 6     0 - 4    BAU  0.004093076
## 7     0 - 4    BAU  1.478125000
## 8     0 - 4    BAU  4.536363636
## 9     0 - 4    BAU  0.023000000
## 10    0 - 4    BAU -0.012037344
## 11    0 - 4    BAU -0.067563144
## 12    0 - 4    BAU -0.029025880
## 13    0 - 4    BAU -0.108185673
## 14    0 - 4    BAU -0.016953394
## 15    0 - 4    BAU  0.018340003
## 16    0 - 4    BAU  0.002415030
## 17    0 - 4    BAU  0.095548935
## 18    0 - 4    BAU  0.720707302
## 19    0 - 4    BAU  0.005217378
## 20    0 - 4    BAU  1.478125000
## 21    0 - 4    BAU  4.386363636
## 22    0 - 4    BAU  0.041000000
## 23    0 - 4    BAU -0.016153270
## 24    0 - 4    BAU -0.089580840
## 25    0 - 4    BAU -0.038463280
## 26    0 - 4    BAU -0.117243257
## 27    5 - 9    BAU  0.002006699
## 28    5 - 9    BAU  0.012270406
## 29    5 - 9    BAU  0.731250000
## 30    5 - 9    BAU  4.259090909
## 31    5 - 9    BAU  0.153000000
## 32    5 - 9    BAU -0.035385812
## 33    5 - 9    BAU -0.111739372
## 34    5 - 9    BAU -0.082525344
## 35    5 - 9    BAU -0.119095150
## 36    5 - 9    BAU  0.002714945
## 37    5 - 9    BAU  0.016601137
## 38    5 - 9    BAU  0.734375000
## 39    5 - 9    BAU  3.709090909
## 40    5 - 9    BAU  0.199000000
## 41    5 - 9    BAU -0.046910387
## 42    5 - 9    BAU -0.143846858
## 43    5 - 9    BAU -0.106356630
## 44    5 - 9    BAU -0.127066194
## 45  10 - 14    BAU  0.003539928
## 46  10 - 14    BAU  0.931250000
## 47  10 - 14    BAU  4.418181818
## 48  10 - 14    BAU  0.139000000
## 49  10 - 14    BAU -0.036171863
## 50  10 - 14    BAU -0.128667793
## 51  10 - 14    BAU -0.084317128
## 52  10 - 14    BAU -0.123494264
## 53  10 - 14    BAU  0.004789314
## 54  10 - 14    BAU  0.931250000
## 55  10 - 14    BAU  3.777272727
## 56  10 - 14    BAU  0.215000000
## 57  10 - 14    BAU -0.048131254
## 58  10 - 14    BAU -0.163487103
## 59  10 - 14    BAU -0.109766315
## 60  10 - 14    BAU -0.130289980
## 61  15 - 19    BAU  0.003539928
## 62  15 - 19    BAU  0.931250000
## 63  15 - 19    BAU  4.418181818
## 64  15 - 19    BAU  0.139000000
## 65  15 - 19    BAU -0.036171863
## 66  15 - 19    BAU -0.128667793
## 67  15 - 19    BAU -0.084317128
## 68  15 - 19    BAU -0.123494264
## 69  15 - 19    BAU  0.004789314
## 70  15 - 19    BAU  0.931250000
## 71  15 - 19    BAU  3.777272727
## 72  15 - 19    BAU  0.215000000
## 73  15 - 19    BAU -0.048131254
## 74  15 - 19    BAU -0.163487103
## 75  15 - 19    BAU -0.109766315
## 76  15 - 19    BAU -0.130289980
## 77  20 - 24    BAU  0.006863037
## 78  20 - 24    BAU  0.737500000
## 79  20 - 24    BAU  3.881818182
## 80  20 - 24    BAU  0.230000000
## 81  20 - 24    BAU -0.053393022
## 82  20 - 24    BAU -0.156347282
## 83  20 - 24    BAU -0.114854981
## 84  20 - 24    BAU -0.116898061
## 85  20 - 24    BAU  0.009285286
## 86  20 - 24    BAU  0.737500000
## 87  20 - 24    BAU  3.363636364
## 88  20 - 24    BAU  0.261000000
## 89  20 - 24    BAU -0.069612671
## 90  20 - 24    BAU -0.192266258
## 91  20 - 24    BAU -0.143629678
## 92  20 - 24    BAU -0.124708188
## 93  25 - 29    BAU  0.006863037
## 94  25 - 29    BAU  0.737500000
## 95  25 - 29    BAU  3.881818182
## 96  25 - 29    BAU  0.230000000
## 97  25 - 29    BAU -0.053393022
## 98  25 - 29    BAU -0.156347282
## 99  25 - 29    BAU -0.114854981
## 100 25 - 29    BAU -0.116898061
## 101 25 - 29    BAU  0.009285286
## 102 25 - 29    BAU  0.737500000
## 103 25 - 29    BAU  3.363636364
## 104 25 - 29    BAU  0.261000000
## 105 25 - 29    BAU -0.069612671
## 106 25 - 29    BAU -0.192266258
## 107 25 - 29    BAU -0.143629678
## 108 25 - 29    BAU -0.124708188
## 109 30 - 34    BAU  0.006863037
## 110 30 - 34    BAU  0.737500000
## 111 30 - 34    BAU  3.881818182
## 112 30 - 34    BAU  0.230000000
## 113 30 - 34    BAU -0.053393022
## 114 30 - 34    BAU -0.156347282
## 115 30 - 34    BAU -0.114854981
## 116 30 - 34    BAU -0.116898061
## 117 30 - 34    BAU  0.009285286
## 118 30 - 34    BAU  0.737500000
## 119 30 - 34    BAU  3.363636364
## 120 30 - 34    BAU  0.261000000
## 121 30 - 34    BAU -0.069612671
## 122 30 - 34    BAU -0.192266258
## 123 30 - 34    BAU -0.143629678
## 124 30 - 34    BAU -0.124708188
## 125 35 - 39    BAU  0.006863037
## 126 35 - 39    BAU  0.737500000
## 127 35 - 39    BAU  3.881818182
## 128 35 - 39    BAU  0.230000000
## 129 35 - 39    BAU -0.053393022
## 130 35 - 39    BAU -0.156347282
## 131 35 - 39    BAU -0.114854981
## 132 35 - 39    BAU -0.116898061
## 133 35 - 39    BAU  0.009285286
## 134 35 - 39    BAU  0.737500000
## 135 35 - 39    BAU  3.363636364
## 136 35 - 39    BAU  0.261000000
## 137 35 - 39    BAU -0.069612671
## 138 35 - 39    BAU -0.192266258
## 139 35 - 39    BAU -0.143629678
## 140 35 - 39    BAU -0.124708188
## 141 40 - 44    BAU  0.006863037
## 142 40 - 44    BAU  0.737500000
## 143 40 - 44    BAU  3.881818182
## 144 40 - 44    BAU  0.230000000
## 145 40 - 44    BAU -0.053393022
## 146 40 - 44    BAU -0.156347282
## 147 40 - 44    BAU -0.114854981
## 148 40 - 44    BAU -0.116898061
## 149 40 - 44    BAU  0.009285286
## 150 40 - 44    BAU  0.737500000
## 151 40 - 44    BAU  3.363636364
## 152 40 - 44    BAU  0.261000000
## 153 40 - 44    BAU -0.069612671
## 154 40 - 44    BAU -0.192266258
## 155 40 - 44    BAU -0.143629678
## 156 40 - 44    BAU -0.124708188
## 157 45 - 49    BAU  0.006863037
## 158 45 - 49    BAU  0.737500000
## 159 45 - 49    BAU  3.881818182
## 160 45 - 49    BAU  0.230000000
## 161 45 - 49    BAU -0.053393022
## 162 45 - 49    BAU -0.156347282
## 163 45 - 49    BAU -0.114854981
## 164 45 - 49    BAU -0.116898061
## 165 45 - 49    BAU  0.009285286
## 166 45 - 49    BAU  0.737500000
## 167 45 - 49    BAU  3.363636364
## 168 45 - 49    BAU  0.261000000
## 169 45 - 49    BAU -0.069612671
## 170 45 - 49    BAU -0.192266258
## 171 45 - 49    BAU -0.143629678
## 172 45 - 49    BAU -0.124708188
## 173 50 - 54    BAU  0.006863037
## 174 50 - 54    BAU  0.737500000
## 175 50 - 54    BAU  3.881818182
## 176 50 - 54    BAU  0.230000000
## 177 50 - 54    BAU -0.053393022
## 178 50 - 54    BAU -0.156347282
## 179 50 - 54    BAU -0.114854981
## 180 50 - 54    BAU -0.116898061
## 181 50 - 54    BAU  0.009285286
## 182 50 - 54    BAU  0.737500000
## 183 50 - 54    BAU  3.363636364
## 184 50 - 54    BAU  0.261000000
## 185 50 - 54    BAU -0.069612671
## 186 50 - 54    BAU -0.192266258
## 187 50 - 54    BAU -0.143629678
## 188 50 - 54    BAU -0.124708188
## 189 55 - 59    BAU  0.006863037
## 190 55 - 59    BAU  0.737500000
## 191 55 - 59    BAU  3.881818182
## 192 55 - 59    BAU  0.230000000
## 193 55 - 59    BAU -0.053393022
## 194 55 - 59    BAU -0.156347282
## 195 55 - 59    BAU -0.114854981
## 196 55 - 59    BAU -0.116898061
## 197 55 - 59    BAU  0.009285286
## 198 55 - 59    BAU  0.737500000
## 199 55 - 59    BAU  3.363636364
## 200 55 - 59    BAU  0.261000000
## 201 55 - 59    BAU -0.069612671
## 202 55 - 59    BAU -0.192266258
## 203 55 - 59    BAU -0.143629678
## 204 55 - 59    BAU -0.124708188
## 205 60 - 64    BAU  0.006863037
## 206 60 - 64    BAU  0.737500000
## 207 60 - 64    BAU  3.881818182
## 208 60 - 64    BAU  0.230000000
## 209 60 - 64    BAU -0.053393022
## 210 60 - 64    BAU -0.156347282
## 211 60 - 64    BAU -0.114854981
## 212 60 - 64    BAU -0.116898061
## 213 60 - 64    BAU  0.009285286
## 214 60 - 64    BAU  0.737500000
## 215 60 - 64    BAU  3.363636364
## 216 60 - 64    BAU  0.261000000
## 217 60 - 64    BAU -0.069612671
## 218 60 - 64    BAU -0.192266258
## 219 60 - 64    BAU -0.143629678
## 220 60 - 64    BAU -0.124708188
## 221 65 - 69    BAU  0.012182501
## 222 65 - 69    BAU  0.978125000
## 223 65 - 69    BAU  3.718181818
## 224 65 - 69    BAU  0.257000000
## 225 65 - 69    BAU -0.066671202
## 226 65 - 69    BAU -0.192635847
## 227 65 - 69    BAU -0.142398049
## 228 65 - 69    BAU -0.120106805
## 229 65 - 69    BAU  0.016482207
## 230 65 - 69    BAU  0.978125000
## 231 65 - 69    BAU  3.259090909
## 232 65 - 69    BAU  0.298000000
## 233 65 - 69    BAU -0.086254649
## 234 65 - 69    BAU -0.232492989
## 235 65 - 69    BAU -0.174942328
## 236 65 - 69    BAU -0.126324226
## 237 70 - 74    BAU  0.012182501
## 238 70 - 74    BAU  0.978125000
## 239 70 - 74    BAU  3.718181818
## 240 70 - 74    BAU  0.257000000
## 241 70 - 74    BAU -0.066671202
## 242 70 - 74    BAU -0.192635847
## 243 70 - 74    BAU -0.142398049
## 244 70 - 74    BAU -0.120106805
## 245 70 - 74    BAU  0.016482207
## 246 70 - 74    BAU  0.978125000
## 247 70 - 74    BAU  3.259090909
## 248 70 - 74    BAU  0.298000000
## 249 70 - 74    BAU -0.086254649
## 250 70 - 74    BAU -0.232492989
## 251 70 - 74    BAU -0.174942328
## 252 70 - 74    BAU -0.126324226
## 253 75 - 79    BAU  0.012182501
## 254 75 - 79    BAU  0.978125000
## 255 75 - 79    BAU  3.718181818
## 256 75 - 79    BAU  0.257000000
## 257 75 - 79    BAU -0.066671202
## 258 75 - 79    BAU -0.192635847
## 259 75 - 79    BAU -0.142398049
## 260 75 - 79    BAU -0.120106805
## 261 75 - 79    BAU  0.016482207
## 262 75 - 79    BAU  0.978125000
## 263 75 - 79    BAU  3.259090909
## 264 75 - 79    BAU  0.298000000
## 265 75 - 79    BAU -0.086254649
## 266 75 - 79    BAU -0.232492989
## 267 75 - 79    BAU -0.174942328
## 268 75 - 79    BAU -0.126324226
## 269 80 - 84    BAU  0.012182501
## 270 80 - 84    BAU  0.978125000
## 271 80 - 84    BAU  3.718181818
## 272 80 - 84    BAU  0.257000000
## 273 80 - 84    BAU -0.066671202
## 274 80 - 84    BAU -0.192635847
## 275 80 - 84    BAU -0.142398049
## 276 80 - 84    BAU -0.120106805
## 277 80 - 84    BAU  0.016482207
## 278 80 - 84    BAU  0.978125000
## 279 80 - 84    BAU  3.259090909
## 280 80 - 84    BAU  0.298000000
## 281 80 - 84    BAU -0.086254649
## 282 80 - 84    BAU -0.232492989
## 283 80 - 84    BAU -0.174942328
## 284 80 - 84    BAU -0.126324226
## 285     85+    BAU  0.012182501
## 286     85+    BAU  0.978125000
## 287     85+    BAU  3.718181818
## 288     85+    BAU  0.257000000
## 289     85+    BAU -0.066671202
## 290     85+    BAU -0.192635847
## 291     85+    BAU -0.142398049
## 292     85+    BAU -0.120106805
## 293     85+    BAU  0.016482207
## 294     85+    BAU  0.978125000
## 295     85+    BAU  3.259090909
## 296     85+    BAU  0.298000000
## 297     85+    BAU -0.086254649
## 298     85+    BAU -0.232492989
## 299     85+    BAU -0.174942328
## 300     85+    BAU -0.126324226
oprint(summary(population,"mean"))
##    Gender       Age    mean
## 1  Female     0 - 4  125040
## 2    Male     0 - 4  130884
## 3  Female     5 - 9  149633
## 4    Male     5 - 9  156654
## 5  Female   10 - 14  151113
## 6    Male   10 - 14  157712
## 7  Female   15 - 19  144441
## 8    Male   15 - 19  152230
## 9  Female   20 - 24  152265
## 10   Male   20 - 24  161679
## 11 Female   25 - 29  172593
## 12   Male   25 - 29  183092
## 13 Female   30 - 34  169653
## 14   Male   30 - 34  181115
## 15 Female   35 - 39  174660
## 16   Male   35 - 39  186122
## 17 Female   40 - 44  168547
## 18   Male   40 - 44  177928
## 19 Female   45 - 49  154391
## 20   Male   45 - 49  159982
## 21 Female   50 - 54  176612
## 22   Male   50 - 54  179182
## 23 Female   55 - 59  185152
## 24   Male   55 - 59  183719
## 25 Female   60 - 64  183336
## 26   Male   60 - 64  176283
## 27 Female   65 - 69  185685
## 28   Male   65 - 69  171275
## 29 Female   70 - 74  186034
## 30   Male   70 - 74  163697
## 31 Female   75 - 79  118190
## 32   Male   75 - 79   93987
## 33 Female   80 - 84   96256
## 34   Male   80 - 84   65140
## 35 Female       85+  103429
## 36   Male       85+   47581
## 37 Female Undefined 2797030
## 38   Male Undefined 2728262
oprint(summary(RR,"mean"))
##     Exposure_agent            Response   ER_function Scaling Gender     Age
## 1             Fish All-cause mortality            RR    None Female   0 - 4
## 2           Omega3       Breast cancer            RR    None Female   0 - 4
## 3             Fish          Depression            RR    None Female   0 - 4
## 4           Omega3      CHD2 mortality Relative Hill    None Female   0 - 4
## 5             Fish All-cause mortality            RR    None   Male   0 - 4
## 6           Omega3       Breast cancer            RR    None   Male   0 - 4
## 7             Fish          Depression            RR    None   Male   0 - 4
## 8           Omega3      CHD2 mortality Relative Hill    None   Male   0 - 4
## 9             Fish All-cause mortality            RR    None Female   5 - 9
## 10          Omega3       Breast cancer            RR    None Female   5 - 9
## 11            Fish          Depression            RR    None Female   5 - 9
## 12          Omega3      CHD2 mortality Relative Hill    None Female   5 - 9
## 13            Fish All-cause mortality            RR    None   Male   5 - 9
## 14          Omega3       Breast cancer            RR    None   Male   5 - 9
## 15            Fish          Depression            RR    None   Male   5 - 9
## 16          Omega3      CHD2 mortality Relative Hill    None   Male   5 - 9
## 17            Fish All-cause mortality            RR    None Female 10 - 14
## 18          Omega3       Breast cancer            RR    None Female 10 - 14
## 19            Fish          Depression            RR    None Female 10 - 14
## 20          Omega3      CHD2 mortality Relative Hill    None Female 10 - 14
## 21            Fish All-cause mortality            RR    None   Male 10 - 14
## 22          Omega3       Breast cancer            RR    None   Male 10 - 14
## 23            Fish          Depression            RR    None   Male 10 - 14
## 24          Omega3      CHD2 mortality Relative Hill    None   Male 10 - 14
## 25            Fish All-cause mortality            RR    None Female 15 - 19
## 26          Omega3       Breast cancer            RR    None Female 15 - 19
## 27            Fish          Depression            RR    None Female 15 - 19
## 28          Omega3      CHD2 mortality Relative Hill    None Female 15 - 19
## 29            Fish All-cause mortality            RR    None   Male 15 - 19
## 30          Omega3       Breast cancer            RR    None   Male 15 - 19
## 31            Fish          Depression            RR    None   Male 15 - 19
## 32          Omega3      CHD2 mortality Relative Hill    None   Male 15 - 19
## 33            Fish All-cause mortality            RR    None Female 20 - 24
## 34          Omega3       Breast cancer            RR    None Female 20 - 24
## 35            Fish          Depression            RR    None Female 20 - 24
## 36          Omega3      CHD2 mortality Relative Hill    None Female 20 - 24
## 37            Fish All-cause mortality            RR    None   Male 20 - 24
## 38          Omega3       Breast cancer            RR    None   Male 20 - 24
## 39            Fish          Depression            RR    None   Male 20 - 24
## 40          Omega3      CHD2 mortality Relative Hill    None   Male 20 - 24
## 41            Fish All-cause mortality            RR    None Female 25 - 29
## 42          Omega3       Breast cancer            RR    None Female 25 - 29
## 43            Fish          Depression            RR    None Female 25 - 29
## 44          Omega3      CHD2 mortality Relative Hill    None Female 25 - 29
## 45            Fish All-cause mortality            RR    None   Male 25 - 29
## 46          Omega3       Breast cancer            RR    None   Male 25 - 29
## 47            Fish          Depression            RR    None   Male 25 - 29
## 48          Omega3      CHD2 mortality Relative Hill    None   Male 25 - 29
## 49            Fish All-cause mortality            RR    None Female 30 - 34
## 50          Omega3       Breast cancer            RR    None Female 30 - 34
## 51            Fish          Depression            RR    None Female 30 - 34
## 52          Omega3      CHD2 mortality Relative Hill    None Female 30 - 34
## 53            Fish All-cause mortality            RR    None   Male 30 - 34
## 54          Omega3       Breast cancer            RR    None   Male 30 - 34
## 55            Fish          Depression            RR    None   Male 30 - 34
## 56          Omega3      CHD2 mortality Relative Hill    None   Male 30 - 34
## 57            Fish All-cause mortality            RR    None Female 35 - 39
## 58          Omega3       Breast cancer            RR    None Female 35 - 39
## 59            Fish          Depression            RR    None Female 35 - 39
## 60          Omega3      CHD2 mortality Relative Hill    None Female 35 - 39
## 61            Fish All-cause mortality            RR    None   Male 35 - 39
## 62          Omega3       Breast cancer            RR    None   Male 35 - 39
## 63            Fish          Depression            RR    None   Male 35 - 39
## 64          Omega3      CHD2 mortality Relative Hill    None   Male 35 - 39
## 65            Fish All-cause mortality            RR    None Female 40 - 44
## 66          Omega3       Breast cancer            RR    None Female 40 - 44
## 67            Fish          Depression            RR    None Female 40 - 44
## 68          Omega3      CHD2 mortality Relative Hill    None Female 40 - 44
## 69            Fish All-cause mortality            RR    None   Male 40 - 44
## 70          Omega3       Breast cancer            RR    None   Male 40 - 44
## 71            Fish          Depression            RR    None   Male 40 - 44
## 72          Omega3      CHD2 mortality Relative Hill    None   Male 40 - 44
## 73            Fish All-cause mortality            RR    None Female 45 - 49
## 74          Omega3       Breast cancer            RR    None Female 45 - 49
## 75            Fish          Depression            RR    None Female 45 - 49
## 76          Omega3      CHD2 mortality Relative Hill    None Female 45 - 49
## 77            Fish All-cause mortality            RR    None   Male 45 - 49
## 78          Omega3       Breast cancer            RR    None   Male 45 - 49
## 79            Fish          Depression            RR    None   Male 45 - 49
## 80          Omega3      CHD2 mortality Relative Hill    None   Male 45 - 49
## 81            Fish All-cause mortality            RR    None Female 50 - 54
## 82          Omega3       Breast cancer            RR    None Female 50 - 54
## 83            Fish          Depression            RR    None Female 50 - 54
## 84          Omega3      CHD2 mortality Relative Hill    None Female 50 - 54
## 85            Fish All-cause mortality            RR    None   Male 50 - 54
## 86          Omega3       Breast cancer            RR    None   Male 50 - 54
## 87            Fish          Depression            RR    None   Male 50 - 54
## 88          Omega3      CHD2 mortality Relative Hill    None   Male 50 - 54
## 89            Fish All-cause mortality            RR    None Female 55 - 59
## 90          Omega3       Breast cancer            RR    None Female 55 - 59
## 91            Fish          Depression            RR    None Female 55 - 59
## 92          Omega3      CHD2 mortality Relative Hill    None Female 55 - 59
## 93            Fish All-cause mortality            RR    None   Male 55 - 59
## 94          Omega3       Breast cancer            RR    None   Male 55 - 59
## 95            Fish          Depression            RR    None   Male 55 - 59
## 96          Omega3      CHD2 mortality Relative Hill    None   Male 55 - 59
## 97            Fish All-cause mortality            RR    None Female 60 - 64
## 98          Omega3       Breast cancer            RR    None Female 60 - 64
## 99            Fish          Depression            RR    None Female 60 - 64
## 100         Omega3      CHD2 mortality Relative Hill    None Female 60 - 64
## 101           Fish All-cause mortality            RR    None   Male 60 - 64
## 102         Omega3       Breast cancer            RR    None   Male 60 - 64
## 103           Fish          Depression            RR    None   Male 60 - 64
## 104         Omega3      CHD2 mortality Relative Hill    None   Male 60 - 64
## 105           Fish All-cause mortality            RR    None Female 65 - 69
## 106         Omega3       Breast cancer            RR    None Female 65 - 69
## 107           Fish          Depression            RR    None Female 65 - 69
## 108         Omega3      CHD2 mortality Relative Hill    None Female 65 - 69
## 109           Fish All-cause mortality            RR    None   Male 65 - 69
## 110         Omega3       Breast cancer            RR    None   Male 65 - 69
## 111           Fish          Depression            RR    None   Male 65 - 69
## 112         Omega3      CHD2 mortality Relative Hill    None   Male 65 - 69
## 113           Fish All-cause mortality            RR    None Female 70 - 74
## 114         Omega3       Breast cancer            RR    None Female 70 - 74
## 115           Fish          Depression            RR    None Female 70 - 74
## 116         Omega3      CHD2 mortality Relative Hill    None Female 70 - 74
## 117           Fish All-cause mortality            RR    None   Male 70 - 74
## 118         Omega3       Breast cancer            RR    None   Male 70 - 74
## 119           Fish          Depression            RR    None   Male 70 - 74
## 120         Omega3      CHD2 mortality Relative Hill    None   Male 70 - 74
## 121           Fish All-cause mortality            RR    None Female 75 - 79
## 122         Omega3       Breast cancer            RR    None Female 75 - 79
## 123           Fish          Depression            RR    None Female 75 - 79
## 124         Omega3      CHD2 mortality Relative Hill    None Female 75 - 79
## 125           Fish All-cause mortality            RR    None   Male 75 - 79
## 126         Omega3       Breast cancer            RR    None   Male 75 - 79
## 127           Fish          Depression            RR    None   Male 75 - 79
## 128         Omega3      CHD2 mortality Relative Hill    None   Male 75 - 79
## 129           Fish All-cause mortality            RR    None Female 80 - 84
## 130         Omega3       Breast cancer            RR    None Female 80 - 84
## 131           Fish          Depression            RR    None Female 80 - 84
## 132         Omega3      CHD2 mortality Relative Hill    None Female 80 - 84
## 133           Fish All-cause mortality            RR    None   Male 80 - 84
## 134         Omega3       Breast cancer            RR    None   Male 80 - 84
## 135           Fish          Depression            RR    None   Male 80 - 84
## 136         Omega3      CHD2 mortality Relative Hill    None   Male 80 - 84
## 137           Fish All-cause mortality            RR    None Female     85+
## 138         Omega3       Breast cancer            RR    None Female     85+
## 139           Fish          Depression            RR    None Female     85+
## 140         Omega3      CHD2 mortality Relative Hill    None Female     85+
## 141           Fish All-cause mortality            RR    None   Male     85+
## 142         Omega3       Breast cancer            RR    None   Male     85+
## 143           Fish          Depression            RR    None   Male     85+
## 144         Omega3      CHD2 mortality Relative Hill    None   Male     85+
##     Exposure      mean
## 1     Direct 0.9879627
## 2     Direct 0.9324418
## 3     Direct 0.9709743
## 4     Direct 0.8918143
## 5     Direct 0.9838467
## 6     Direct 0.9104266
## 7     Direct 0.9615370
## 8     Direct 0.8827567
## 9     Direct 0.9646142
## 10    Direct 0.8884120
## 11    Direct 0.9174749
## 12    Direct 0.8809048
## 13    Direct 0.9530896
## 14    Direct 0.8563940
## 15    Direct 0.8936439
## 16    Direct 0.8729338
## 17    Direct 0.9638281
## 18    Direct 0.8714502
## 19    Direct 0.9156829
## 20    Direct 0.8765057
## 21    Direct 0.9518687
## 22    Direct 0.8367223
## 23    Direct 0.8902337
## 24    Direct 0.8697100
## 25    Direct 0.9638281
## 26    Direct 0.8714502
## 27    Direct 0.9156829
## 28    Direct 0.8765057
## 29    Direct 0.9518687
## 30    Direct 0.8367223
## 31    Direct 0.8902337
## 32    Direct 0.8697100
## 33    Direct 0.9466070
## 34    Direct 0.8437245
## 35    Direct 0.8851684
## 36    Direct 0.8831019
## 37    Direct 0.9303873
## 38    Direct 0.8078619
## 39    Direct 0.8564131
## 40    Direct 0.8752918
## 41    Direct 0.9466070
## 42    Direct 0.8437245
## 43    Direct 0.8851684
## 44    Direct 0.8831019
## 45    Direct 0.9303873
## 46    Direct 0.8078619
## 47    Direct 0.8564131
## 48    Direct 0.8752918
## 49    Direct 0.9466070
## 50    Direct 0.8437245
## 51    Direct 0.8851684
## 52    Direct 0.8831019
## 53    Direct 0.9303873
## 54    Direct 0.8078619
## 55    Direct 0.8564131
## 56    Direct 0.8752918
## 57    Direct 0.9466070
## 58    Direct 0.8437245
## 59    Direct 0.8851684
## 60    Direct 0.8831019
## 61    Direct 0.9303873
## 62    Direct 0.8078619
## 63    Direct 0.8564131
## 64    Direct 0.8752918
## 65    Direct 0.9466070
## 66    Direct 0.8437245
## 67    Direct 0.8851684
## 68    Direct 0.8831019
## 69    Direct 0.9303873
## 70    Direct 0.8078619
## 71    Direct 0.8564131
## 72    Direct 0.8752918
## 73    Direct 0.9466070
## 74    Direct 0.8437245
## 75    Direct 0.8851684
## 76    Direct 0.8831019
## 77    Direct 0.9303873
## 78    Direct 0.8078619
## 79    Direct 0.8564131
## 80    Direct 0.8752918
## 81    Direct 0.9466070
## 82    Direct 0.8437245
## 83    Direct 0.8851684
## 84    Direct 0.8831019
## 85    Direct 0.9303873
## 86    Direct 0.8078619
## 87    Direct 0.8564131
## 88    Direct 0.8752918
## 89    Direct 0.9466070
## 90    Direct 0.8437245
## 91    Direct 0.8851684
## 92    Direct 0.8831019
## 93    Direct 0.9303873
## 94    Direct 0.8078619
## 95    Direct 0.8564131
## 96    Direct 0.8752918
## 97    Direct 0.9466070
## 98    Direct 0.8437245
## 99    Direct 0.8851684
## 100   Direct 0.8831019
## 101   Direct 0.9303873
## 102   Direct 0.8078619
## 103   Direct 0.8564131
## 104   Direct 0.8752918
## 105   Direct 0.9333288
## 106   Direct 0.8075121
## 107   Direct 0.8576020
## 108   Direct 0.8798932
## 109   Direct 0.9137454
## 110   Direct 0.7677808
## 111   Direct 0.8250577
## 112   Direct 0.8736758
## 113   Direct 0.9333288
## 114   Direct 0.8075121
## 115   Direct 0.8576020
## 116   Direct 0.8798932
## 117   Direct 0.9137454
## 118   Direct 0.7677808
## 119   Direct 0.8250577
## 120   Direct 0.8736758
## 121   Direct 0.9333288
## 122   Direct 0.8075121
## 123   Direct 0.8576020
## 124   Direct 0.8798932
## 125   Direct 0.9137454
## 126   Direct 0.7677808
## 127   Direct 0.8250577
## 128   Direct 0.8736758
## 129   Direct 0.9333288
## 130   Direct 0.8075121
## 131   Direct 0.8576020
## 132   Direct 0.8798932
## 133   Direct 0.9137454
## 134   Direct 0.7677808
## 135   Direct 0.8250577
## 136   Direct 0.8736758
## 137   Direct 0.9333288
## 138   Direct 0.8075121
## 139   Direct 0.8576020
## 140   Direct 0.8798932
## 141   Direct 0.9137454
## 142   Direct 0.7677808
## 143   Direct 0.8250577
## 144   Direct 0.8736758
###################
# Graphs
ggplot(trim(amount), aes(x=Group, weight=amountResult, fill=Fish))+geom_bar()+
  facet_wrap(~Gender)+
  labs(
    title="Kalansyönti Suomessa ikäryhmittäin",
    y="Syönti (g/d)"
  )

ggsave("Kalansyönti Suomessa ikäryhmittäin.svg")
## Saving 7 x 5 in image
plot_ly(trim(EvalOutput(total_amount)), x=~Kala, y=~total_amountResult, color=~Kala, type="bar") %>%
  layout(yaxis=list(title="Kalan kokonaiskulutus Suomessa (milj kg /a)"), barmode="stack")
## Warning: `arrange_()` is deprecated as of dplyr 0.7.0.
## Please use `arrange()` instead.
## See vignette('programming') for more help
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_warnings()` to see where this warning was generated.
plot_ly(trim(conc_vit), x=~Nutrient, y=~conc_vitResult, color=~Kala, type="scatter", mode="markers") %>%
  layout(yaxis=list(title="Concentrations of nutrients (mg or ug /g)"))
## Warning in RColorBrewer::brewer.pal(N, "Set2"): n too large, allowed maximum for palette Set2 is 8
## Returning the palette you asked for with that many colors
## Warning in RColorBrewer::brewer.pal(N, "Set2"): n too large, allowed maximum for palette Set2 is 8
## Returning the palette you asked for with that many colors
ggplot(conc_mehg@output, aes(x=conc_mehgResult, colour=Area))+stat_ecdf()+
  scale_x_log10()+facet_wrap(~Fish, scales="free_x")

ggplot(trim(conc_pfas), aes(x=Fish, y=conc_pfasResult))+geom_point()+ # Inputed data for missing species.
  coord_flip()

ggplot(conc_pfas@output, aes(x=conc_pfasResult, color=Compound, linetype=Area))+stat_ecdf()+
  scale_x_log10()+
  stat_ecdf(data=conc_eukalat@output, aes(x=conc_eukalatResult))+
  scale_linetype_manual(values=c("dotted","solid","twodash"))+
  labs(
    title="PFAS concentration in fishes in Finland",
    x="PFAS concentration (ng/g fresh weight)",
    y="Cumulative probability"
  )
## Warning in self$trans$transform(x): NaNs produced
## Warning: Transformation introduced infinite values in continuous x-axis
## Warning: Removed 27 rows containing non-finite values (stat_ecdf).

# The code may produce some negative values, which are removed from the graph
ggsave("PFAS-pitoisuus kalassa Suomessa.svg")
## Saving 7 x 5 in image
## Warning in self$trans$transform(x): NaNs produced
## Warning: Transformation introduced infinite values in continuous x-axis
## Warning: Removed 27 rows containing non-finite values (stat_ecdf).
ggplot(conc@output, aes(x=concResult, colour=Fish))+stat_ecdf()+
  facet_wrap(~Compound, scales="free_x")

#ggplot(oapply(expo_dir, cols=c("Iter"),FUN=mean)@output,
#       aes(x=Gender, weight=expo_dirResult,fill=Fish))+geom_bar()+
#  facet_wrap(~Compound, scales="free_y")+
#  labs(title="Eri yhdisteiden saanti kalasta")
#ggsave("Yhdisteiden saanti kalasta Suomessa.svg")

plot_ly(trim(oapply(exposure,NULL,mean,"Gender")), x=~Age, y=~exposureResult, color=~Exposure_agent, text=~Exposure_agent, type="bar") %>%
  layout(yaxis=list(title="Exposure to nutrients (g or ug /d)"))
## Warning in RColorBrewer::brewer.pal(N, "Set2"): n too large, allowed maximum for palette Set2 is 8
## Returning the palette you asked for with that many colors
## Warning in RColorBrewer::brewer.pal(N, "Set2"): n too large, allowed maximum for palette Set2 is 8
## Returning the palette you asked for with that many colors
ggplot(trim(oapply(exposure,NULL,mean,"Gender")), aes(x=Age, weight=exposureResult))+geom_bar()+
  facet_wrap(~Exposure_agent, scales="free_y")+
  labs(
    title="Exposure to compounds",
    y="(omega: mg/d; vit D: ug/d, PFAS: ng/d)"
  )

cat("Kalaperäisiä tautitaakkoja Suomessa\n")
## Kalaperäisiä tautitaakkoja Suomessa
if(openv$N>1) {
  tmp <- summary(oapply(BoDattr,NULL,sum,c("Age","Gender","Response")))
  tmp <- data.frame(
    Altiste = tmp$Exposure_agent,
#    Vaikutus = tmp$Response,
    Keskiarvo = as.character(signif(tmp$mean,2)),
    "95 luottamusväli" = paste0(signif(tmp$Q0.025,2)," - ", signif(tmp$Q0.975,2)),
    Keskihajonta = signif(tmp$sd,2)
  )#[rev(match(lev, tmp$Exposure_agent)),]

  oprint(tmp)
  
  tmp <- summary(oapply(BoDattr,NULL,sum,c("Age","Gender","Exposure_agent")))
  tmp <- data.frame(
    Terveysvaikutus = tmp$Response,
    Keskiarvo = signif(tmp$mean,2),
    "95 luottamusväli" = paste0(signif(tmp$Q0.025,2)," - ", signif(tmp$Q0.975,2)),
    Keskihajonta = signif(tmp$sd,2)
  )
  oprint(tmp)
}
##     Altiste Keskiarvo X95.luottamusväli Keskihajonta
## 1       DHA      -500       -1900 - -77          450
## 2      Fish    -68000     -260000 - -30        76000
## 3      MeHg       540          0 - 4800         1400
## 4    Omega3    -46000   -83000 - -16000        18000
## 5      PFAS      3000         0 - 24000         6100
## 6       TEQ      7600         0 - 38000        12000
## 7 Vitamin D     10000       440 - 48000        13000
##                                      Terveysvaikutus Keskiarvo
## 1                                All-cause mortality    -62000
## 2                                      Breast cancer     -4200
## 3                            Cancer morbidity yearly      3400
## 4                                     CHD2 mortality    -42000
## 5                                         Depression     -5800
## 6  Dioxin recommendation tolerable daily intake 2018      3000
## 7                                  Immunosuppression        11
## 8                          Loss in child's IQ points        47
## 9                                           PFAS TWI      3000
## 10                               Sperm concentration       730
## 11                          Vitamin D recommendation     10000
## 12                           Yes or no dental defect       450
##    X95.luottamusväli Keskihajonta
## 1      -250000 - -25        72000
## 2        -14000 - 14         4100
## 3          0 - 23000         8600
## 4    -73000 - -16000        15000
## 5         -22000 - 0         6300
## 6          0 - 21000         5900
## 7             0 - 59           17
## 8        -800 - 3200         1100
## 9          0 - 24000         6100
## 10        -11 - 6200         1800
## 11       440 - 48000        13000
## 12          0 - 3700         1000
ggplot(trim(BoDattr), aes(x=Exposure_agent, weight=BoDattrResult, fill=Response))+geom_bar()

ggplot(trim(BoDattr), aes(x=Response, weight=BoDattrResult, fill=Exposure_agent))+geom_bar()

ggplot(trim(BoDattr), aes(x=Age, weight=BoDattrResult, fill=Response))+geom_bar(position="stack")+
  scale_x_discrete(breaks = levels(ages)[seq(1,length(ages),3)])

ggplot(trim(BoDattr), aes(x=Age, weight=BoDattrResult, fill=Response))+geom_bar(position="stack")+
  scale_x_discrete(breaks = levels(ages)[seq(1,length(ages),3)])+
  coord_cartesian(ylim=c(-2000,2000))

plot_ly(trim(BoDattr), x=~Age, y=~BoDattrResult, color=~Response, text=~paste(Age, Exposure_agent, sep=": "), type="bar") %>%
  layout(yaxis=list(title="Disease burden (DALY /a); CHD2=coronary heart disease"), barmode="stack")
## Warning in RColorBrewer::brewer.pal(N, "Set2"): n too large, allowed maximum for palette Set2 is 8
## Returning the palette you asked for with that many colors

## Warning in RColorBrewer::brewer.pal(N, "Set2"): n too large, allowed maximum for palette Set2 is 8
## Returning the palette you asked for with that many colors
plot_ly(trim(BoDattr), x=~Exposure_agent, y=~BoDattrResult, color=~Response, text=~paste(Age, Exposure_agent, sep=": "), type="bar") %>%
  layout(yaxis=list(title="Disease burden (DALY /a); CHD2=coronary heart disease"), barmode="stack")
## Warning in RColorBrewer::brewer.pal(N, "Set2"): n too large, allowed maximum for palette Set2 is 8
## Returning the palette you asked for with that many colors

## Warning in RColorBrewer::brewer.pal(N, "Set2"): n too large, allowed maximum for palette Set2 is 8
## Returning the palette you asked for with that many colors
################ Insight network
#gr <- scrape(type="assessment")
#objects.latest("Op_en3861", "makeGraph") # [[Insight network]]
#gr <- makeGraph(gr)
#export_graph(gr, "ruori.svg")
#render_graph(gr) # Does not work: Error in generate_dot(graph) : object 'attribute' not found
##################### Diagnostics
objects.latest("Op_en6007", code_name="diagnostics")
## Loading objects:
##   showind
##   binoptest
##   showLoctable
##   ovashapetest
showLoctable()
showind()
## Hg is not an ovariable.
## subgrouping is not an ovariable.
## mc2d is not an ovariable.
## sumExposcen is not an ovariable.
## mc2d is not an ovariable.
## mc2dparam is not an ovariable.
## TEFversion is not an ovariable.
# Run all above